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1.  INTRODUCTION 


The  progress  achieved  in  the  design  and  fabrication  of  a  GaAlAs 
semiconductor  laser  capable  of  high  average  power  levels  (0.1  to  1.0  W)  and 
low  (~1  mrad)  beam  divergence  is  reported.  The  progran  began  in  June  1977 
with  a  three-month  design  phase.  The  design  included  (1)  the  development  of  a 
laser  array  concept,  (2)  an  analysis  of  laser  injection  locking,  (3)  a  thermal 
analysis  of  the  suggested  array,  (4)  the  development  of  an  integrated  optical 
beam  combining  network,  and  (3)  the  development  of  an  epitaxial  structure 
design. 

The  results  of  the  initial  analyses  phase  suggested  the  use  of  several 
relatively  low-powered  GaAs  lasers  coherently  locked  together  in  an  array  to 
obtain  the  desired  optical  output,  and  the  use  of  waveguide  corrugations  to 
perform  optical  beam  forming.  A  detailed  description  of  the  design  study  is 
included  in  the  first  year  final  report.1  The  design  concept  is  illustrated 
in  Fig.  1-1.  The  lasing  cavity  is  defined  by  the  two  distributed  Bragg 
reflectors  (DBR)  at  right  angles  to  each  other.  The  distributed  Bragg  deflec¬ 
tor  (DBD)  cascade  aids  in  locking  the  laser  array  and  expands  the  beam  in  the 
vertical  direction.  The  distributed  output  coupler  provides  beam  expansion  in 
the  horizontal  direction,  and  the  resultant  laser  output  aperture  has  an  area 
of  several  millimeters  compared  to  a  few  microns  for  conventional  Fabry-Perot 
semiconductor  lasers  (Fig.  1-2).  It  is  the  large  increase  in  output  aperture 
of  the  high-brightness  laser  that  will  result  in  low  beam  divergence. 

The  high-brightness  laser  project  has  been  an  ambitious  undertaking. 
Multilayer  GaAlAs  material  had  to  be  grown  to  micron  and  subnicron  dimensions 
over  centimeter  square  areas,  evaluated,  and  fabricated  into  complex  devices. 

The  high-brightness  laser  shown  in  Figs.  1-1  and  1-2  is  actually  many 
devices  integrated  onto  a  single  crystal:  DBR  lasers,  a  DBD  cascade,  and  an 


Hl.M.  Stoll,  High  Brightness  GaAs  Laser  Source,  Contract  No.  77G-101820-000, 
The  Aerospace  Corporation  (30  June  1978). 
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output  coupler.  At  Aerospace,  we  have  succeeded  in  fabricating  the  individual 
components.  The  fabrication  techniques  and  experimental  results,  lor  the  indi¬ 
vidual  devices  are  discussed  in  Section  4,  as  are  some  results  on  coupling  of 
stripe  geometry  lasers.  We  have  not,  however,  advanced  to  the  point  of 
integrating  the  components  together. 

It  is  our  belief  that  the  successful  fabrication  of  a  high-brightness 
laser  will  profoundly  impact  optical  systems  of  the  future.  For  this  reason, 
Aerospace  is  investing  corporate  funds  to  support  the  basic  research  in  the 
materials  growth  area  required  for  sustained  progress.  In  parallel  with  the 
materials  effort,  further  refining  and  optimization  of  the  fabrication 
techniques  required  to  produce  the  key  elements  of  the  high-brightness  laser 
is  under  way. 
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2.  SUPPORTING  MATERIALS  TECHNOLOGY 


2.1  LIQUID  PHASE  EPITAXY  SYSTEM 

The  Aerospace  liquid  phase  epitaxy  (LPE)  system  is  shown  schematically  in 
Fig.  2-1.  Photographs  of  the  system  are  shown  in  Fig.  2-2.  The  important 
features  are  a  temperature-controlled  furnace,  a  high-purity  gas-flow  system, 
and  a  graphite  boat  and  slider. 

The  furnace  is  a  three-zone  Thermco  minibrute  with  a  Eurotherm  programmer 
having  a  quoted  stability  of  ±  0.25°C.  In  fact,  we  measured  a  temperature 
stability  of  ±  0.1°C  over  typical  growth  times.  After  careful  adjustment,  we 
obtained  the  temperature  profile  shown  in  Fig.  2-3.  The  temperature  was  the 
sane  within  0.33°C  over  the  full  9  in.  of  the  boat  and  slider.  The  programmer 
dropped  the  temperature  at  rates  typically  l°C/min  but  was  capable  of  rates  as 
slow  as  0.2°C/min. 

A  high-purity  H2  stream  at  a  rate  of  10  cm'Vmin  (5  standard  ft^/hr)  was 
produced  from  commercial  hydrogen  by  diffusing  through  a  palladium  hydrogen 
purifier,  resulting  in  —0.1  ppn  impurities.  The  gas  flow  system  also  included 
N2  flushing  capabilities,  and  dust,  oil,  and  preliminary  oxygen-  removing 
filters.  Figure  2-4  is  a  photograph  of  the  graphite  boat  used  for  DH  laser 
growth.  This  is  machined  out  of  high-density  POCO  graphite  (DP-7)  and  placed 
in  the  center  of  the  furnace.  An  unusual  feature  of  this  particular  boat  is 
evident,  i.e.,  the  two  sapphire  plates  used  as  wiper  blades  to  remove  melt. 
As  part  of  this  contract,  new  and  more  sophisticated  growth  techniques  were 
developed.  The  growth  procedure  used  now  will  be  described  in  detail  in 
Appendix  A. 

The  typical  double  heterostructure  laser  structure  is  shown  in  Fig.  2-5, 
including  an  active  GaAs  layer  of  thickness  t,  surrounded  by  Gaj_xAlxAs 
several  microns  thick  with  an  AlAs  concentration  of  x  =  0.47.  The  structure 
also  includes  a  cap  of  GaAs,  which  occurred  as  a  consequence  of  melt  remaining 
on  the  surface  after  growth  and  causing  GaAs  deposition  as  the  sample  cooled. 
A  SEM  photograph  of  a  cleaved  and  stained  cross  section  of  a  typical  double 
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heterostructure  is  shown  in  Fig.  2-5.  The  stain  used  for  this  photograph  was 
Monsanto  etch  (1:1:4;  HF^NC^l^O)  for  10  sec. 

After  the  double  heterostructures  were  grown,  ohmic  contacts  were  fabri¬ 
cated.  First,  500  A  of  chromium  and  then  200  A  of  gold  were  evaporated  on  the 
p-side,  with  the  substrate  at  300°C.  Then,  the  samples  were  thinned  to  7  mil 
thick,  and  2000  A  of  AuGe  was  evaporated  on  the  p-side,  with  the  sample  at 
180°C.  These  contacts  were  then  alloyed  by  heating  to  300°C  for  3  min  while 
still  in  the  vacuum  system.  As  part  of  this  contract,  new  ohmic  contact  tech¬ 
niques  were  developed,  which  will  be  described  later  in  this  report.  After 
application  of  ohmic  contacts,  the  wafer  was  cleaved  into  bars  16  nils  wide 
and  then  scribed  into  dice  7  mils  wide.  The  dice  then  were  mounted  on  a  TO-5 
header  with  indium  solder,  and  a  gold  wire  lead  was  attached  with  a  therno- 
compression  bonder. 

Laser  thresholds  were  measured  by  monitoring  the  light  output  through  a 
monochromator.  The  results  are  summarized  in  Table  2-1.  The  results  demon¬ 
strated  that  we  could  easily  and  reproducihly  achieve  cleaved  double¬ 
heterostructure  GaAs  lasers  with  state-of-the-art  properties,  and  work  was 
begun  on  development  of  a  high-brightness  double  heterostructure  (DH)  laser 
using  distributed  Bragg  reflector  (DBR)  structures,  a  beam  expander,  and  an 
output  collimator. 

First,  a  preliminary  design  study  was  carried  out.  The  results  of  this 
study  indicated  that  the  integrated  optics  laser  source  could  most  conve¬ 
niently  be  fabricated  by  means  of  chemical  etch  techniques.  The  structure 
decided  upon  is  shown  in  Fig.  2-6,  and  is  called  an  asymmetric  large  optical 
cavity  (ALOC). 

The  high-brightness  integrated-optics  laser  source  requires  an  active 
laser  layer  coupled  to  a  passive  waveguide  layer  (the  region  in  which  beam 
expansion  takes  place).  The  passive  waveguide  layer  must  have  a  composition 
different  from  the  laser  layer  in  order  to  minimize  absorption  loss  and  in¬ 
efficient  operation.  The  need  for  different  layer  compositions  requires  the 
use  of  different  melts  during  the  liquid-phase  epitaxy  (LFE)  epigrowth  pro¬ 
cess.  In  the  original  proposal,  it  was  suggested  that  this  could  be  achieved 
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Table  2-1.  Experimental  Results  for  Double  Heterostructure  Lasers 


Wafer 

X 

t,  pm 

Average 

Measured 

Threshold, 

kA/cnr 

Average 

Measured 

Threshold, 

kA/cnr 

Wave¬ 

length 

X,A 

L7 

0.47 

0.2 

2.4  ±  0.5 

1.9 

8750  ±  50 

L8 

0.47 

0.6 

3.8  ±  0.8 

2.0 

8825  ±  25 

L10 

0.47 

0.1 

2.6  ±  0.9 

2.0 

8600  ±  20 

LI  1 

0.44 

0.46 

4.4  ±  1.5 

2.9 

8800  ±  70 

GaAs 


Fig.  2-6.  Fptlayer  Structure,  Integrating  Laser,  and  Passive  Waveguide  on 
Same  Substrate.  Planar  epitaxy  with  chemical  etch  and  ion 
milling  used  for  measurements. 


by  nonplanar  epitaxy.  During  the  design  study,  however,  it  was  discovered 
that  the  ALOC  structure  allows  for  planar  layers  to  be  grown  and  then  the 
active-passive  transition  to  be  fabricated  later  by  chemically  etching  through 
a  photoresist  mask.  In  the  active  region,  the  guided  light  is  not  confined 
only  to  the  laser  layer  but  extends  throughout  the  waveguide  layer.  When  the 
absorbing  GaAs  layer  in  the  passive  region  is  removed,  the  light  propagates  in 
the  waveguide  layer.  A  detailed  analysis  of  the  design  of  the  DBR  structure 
is  included  as  Appendix  B. 

The  result  of  this  analysis  was  the  calculation  that  the  layer  structure 
shown  in  Fig.  2-7  should  have  a  threshold  less  than  2  kamp/cm^  and  an  effi¬ 
ciency  87%  that  of  typical  cleaved  lasers.  An  additional  conclusion  of  the 
analysis  is  that  there  is  a  five-layer  symmetric  large  optical  cavity  (SLOC) 
structure,  which  has  an  even  lower  threshold  and  higher  efficiency.  However, 
because  of  the  greater  difficulty  in  fabrication,  we  decided  to  make  our  first 
high-br ightness  lasers  using  the  asymmetric  large  optical  cavity  (ALOC)  struc¬ 
ture.  Later  research  efforts  to  increase  efficiency  of  the  devices  can  be 
directed  toward  development  of  the  SLOC. 

The  decision  was  made  to  use  the  ALOC  structure,  and,  accordingly, 
crystal  growth  was  begun  in  order  to  fabricate  these  epilayers.  Six  wafers 
were  grown,  and  cleaved  lasers  were  fabricated.  The  structure  is  shown  in 
Fig.  2-7.  The  results  are  summarized  in  Table  2-2.  Lasers  L15  and  L16  pro¬ 
duced  the  expected  results,  that  is,  the  thresholds  for  the  ALOC  were  lower 
than  any  of  the  DH  lasers.  In  fact.  Fig.  2-8  demonstrates  this  dramatically 
by  showing  the  threshold  plotted  on  the  same  scale  as  the  DH  lasers.  The 
results  of  these  LI  5  and  1.16  ALOC  lasers  gave  us  confidence  in  the  basic 
design.  There  were  problems,  however,  with  samples  L17  through  L19,  indica¬ 
ting  that  there  was  something  wrong  with  the  thin  active  regions.  Later  the 
cause  was  discovered.  It  was  attributed  to  the  use  of  a  too  high  aluminum 
concentration  in  the  layer  above  the  active  region,  which  caused  strain  in 
very  thin  active  regions. 

In  parallel  with  these  crystal  growth  runs,  we  began  to  study  selective 
chemical  etching  of  the  epilayers.  For  this  purpose,  we  grew  layers  L20 
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Fig.  2-7.  ALOC  Laser  Layer  Structure 


Table  2-2.  ALOC  Laser 


um  J= 


cra  kA/cn 


2.8  ±  1.1 
1.6  ±  0.4 
No  lasing 
4.0 

No  lasing 


8760  ±  1000 
8750  ±  30 


5] 


THRESHOLD,  kamp/cm' 


ACTIVE  REGION  THICKNESS, /xm 


Fig,  2-8.  DH  Laser 
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through  L24,  which  were  not  double  heterostructures,  but  were  layers  of  high 
AlAs  concentration  grown  on  layers  of  low  AlAs  concentration. 

In  all  these  runs,  we  found  that  the  layers  we  grew  had  microscopic  pits, 
which  did  not  seem  to  interfere  with  cleaved  laser  operation,  but  which  were 
unacceptable  for  the  high-brightness  laser.  It  became  clear  that  we  required 
more  sophisticated  crystal  growth  techniques.  We  attributed  these  pits  to 
aluminum  oxide,  and  decided  to  develop  a  new  boat  that  would  permit  us  to  add 
the  aluminum  dopant  only  after  an  initial  prebakeout  period.  We  followed  a 
procedure  outlined  by  one  of  the  research  scientists  at  RCA.*  This  boat  was 
designed  and  fabricated  and  is  shown  in  Fig.  2-9. 

A  sliding  one-piece  lid  was  made  for  the  boat,  which  would  drop  into  the 
melt  wells  when  positioned  properly.  This  lid  is  shown  at  the  bottom  of  Fig. 
2-9(a).  The  boat  is  loaded  before  growth,  with  the  aluminum  dopants  sitting 
on  the  edges  between  melt  wells.  The  lid  rests  on  top  of  the  melt  wells, 
positioned  so  as  to  be  able  to  push  the  aluminum  into  the  melts  when  re¬ 
quired.  After  a  preliminary  bakeout  the  lid  is  pushed  over  the  melts,  the 
aluminum  falLs  in,  and  the  lid  drops  down  onto  the  melt.  The  melts  are  then 
taken  to  growth  temperature. 

Experiments  began  with  the  newly  designed  boat.  Crystals  Nl  through  N4 
were  simple  double  heterostructure  lasers,  and  01  through  04  were  grown  as 
low-threshold  symmetric  large  optical  cavities.  It  appeared  that  most  of  the 
pits  were  eliminated,  but  there  appeared  to  be  excess  GaAs  on  the  surface, 
which  made  it  difficult  to  use  for  the  chemical  etching. 

We  decided  that  the  best  way  to  provide  a  good  quality  thin  cap  was  to 
grow  a  surface  layer  of  Ga^  7A1q  gAs  on  top  of  the  cap.  It  would  not  matter 
if  this  layer  were  uneven  since  the  very  first  step  in  the  processing  would  be 
tiio  removal  of  this  layer  hv  HC1. 

Wafers  M  through  A 1 2  were  grown  as  ALOC  lasers  with  this  AlAs  surface 
layer.  These  layers  were  verv  flat  and  were  of  the  desired  thickness.  A 


M.  Ktteuherg,  RCA  Corporation,  private  communication. 
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typical  SUM  of  a  cleaved  and  stained  wafer  is  shown  in  Fig.  2-10.  Processing 
was  begun  on  all  these  wafers  before  making  the  active  or  passive  lasers.  The 
layer  structures  grown  are  summarized  in  Table  2-3. 

We  also  decided  that  one  crystal  growth  furnace  was  insufficient  for  the 
demands  of  the  program,  so  a  second  identical  furnace  was  introduced  and 
plumbed. 

Next,  wafers  M3  through  Alb  were  grown.  These  were  ALOC  structures 
essentially  identical  to  A1  through  A12  and  were  intended  for  DBR  lasers.  The 
layer  surfaces  were  smooth,  and  the  thicknesses  were  as  expected.  During  this 
period  of  growth,  it  was  known  that  the  current-voltage  (I-V)  characteristics 
of  lasers  made  for  ALOC  material  were  poor,  but  the  reason  for  this  was  not 
understood.  After  the  growth  of  Alb,  however,  it  was  discovered  that  they 
were  poor  because  the  surface  AlAs  layer  that  was  to  have  been  etched  away  was 
not  in  fact  fully  removed.  Because  we  had  anticipated  its  removal,  this  layer 
had  not  been  electrically  doped.  The  incomplete  removal  resulted,  therefore 
in  an  i-laver  instead  of  an  ohmic  contact.  With  this  fact  now  fully  under¬ 
stood,  the  ALOC  heterostructures  were  grown  with  p-doped  surface  layers. 
Layers  \17  through  A23  were  grown  with  this  configuration. 

As  a  result  of  later  discoveries,  that  is,  the  AB  etch  removed  layers 
with  lower  AlAs  concent  rat  Ion,  the  AlAs  concentration  was  reduced  for  A20 
through  \ 2 2 .  When  it  was  discovered  that  thin  active  and  waveguide  regions 
required  even  lower  AlAs  concentrations  to  lase,  the  last  sample  in  this 
series  was  grown  witli  considerably  lower  AlAs  concentrations.  A  sunmary  of 
the  samples  grown  and  their  structure  is  provided  in  Table  2-4. 

As  a  result  of  passive  waveguide  tests,  we  decided  to  grow  a  buffer  layer 
of  GaAs  on  the  substrate  before  the  first  GaAlAs  layer  was  grown.  There  is 
evidence  in  the  literature  that  this  would  reduce  the  dislocations  and  improve 
the  quaiitv  of  the  passive  waveguides,  a  necessity  For  efficient  high¬ 
brightness  lasers.  Therefore,  crystals  Cl  through  Cb  were  grown.  The  layer 
structure  for  the  Olivers  is  given  in  Table  2-3.  Because  the  buffer  layer 
required  the  use  of  one  of  the  melt  holes  of  the  boat  designed  for  the  ALOC 
structure,  the  surface  1  ivr  could  not  he  grown.  As  a  result,  the  cap  layer 
was  verv  thick,  arid  I  hose  wafer's  were  not  suitable  for  DBF  lasers. 
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Table  2-3.  Layer  Structure  of  Wafers 


n-Isolation  Waveguide  Active  p-Isolation  Cap 


Wafer 

X 

t,  pm 

X 

t,  pm 

t ,  pm 

X 

t,  pm 

t,  pm 

A1 

0.3 

3.6 

0.29 

0.43 

0.11 

0.6 

8.6 

1 

A2 

0.3 

12 

0.29 

1 

— 

0.6 

8 

1 

A3 

0.3 

13 

0.29 

2 

— 

0.6 

6 

3 

A4 

0.2 

8 

0.18 

1 

— 

0.6 

9 

1 

A5a 

A6 

0.32 

8b 

0.29 

3b 

0.46b 

0.6 

4.6 

2 

A7 

0.32 

3.9 

0.29 

0.7 

— 

0.6 

3.3 

3 

A8 

0.32 

2 

0.29 

1.6 

0.46 

0.6 

2.3 

3 

A9 

0.32 

3 

0.29 

1.2 

0.3 

0.6 

5 

3 

A10C 

All 

0.32 

3.8 

0.25 

_ 

0.2 

0.6 

5.1 

1.8 

A12 

0.4 

3.7 

0.9 

1.2 

0.6 

0.6 

7 

1.8 

aBad  growth,  need  to  readjust  furnace  profile. 
bGrown  thick  to  assist  in  identifying  layers. 
cBad  growth,  mistake. 
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Table  2-4.  Layer  Structure  for  ALOC  Wafers 


Wafer 

n-Isolation 
x  t ,  pm 

Waveguide 
x  t ,  pm 

Active 
t,  pm 

p-Isolation 
x  t ,  pm 

Cap 
t,  pm 

A13a 

Al4a 

Al  5 

0.6 

0.27 

0.3 

1.1 

0.15 

0.62 

2.8 

3.5 

A16b 

Al  7 

0.4 

3.14 

0.29 

1.1 

0.27 

0.62 

3.5 

0.8 

A18 

0.4 

3.5 

0.29 

1.1 

0.27 

0.62 

3.8 

0.4 

A19 

0.4 

2.3 

0.29 

0.96 

0.5 

0.62 

2.5 

2.4 

A20 

0.3 

3.5 

0.2 

— 

0.5 

0.5 

4 

0.5 

A21 

0.3 

3.6 

0.2 

0.4 

— 

0.5 

4.8 

0.6 

A22 

0.3 

2.6 

0.2 

0.6 

0.15 

0.5 

2.8 

1.6 

A23 

0.3 

1.2 

0.1 

0.6 

0.15 

0.3 

3.2 

0.5 

^Aluminum  did  not  fall  in. 
DOxygen  leak. 
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Table  2-5.  DBR  Structure  E-Series 


A  new  boat  was  therefore  designed  with  an  additional  melt  hole  for  the 
buffer  layer.  This  boat  was  also  designed  to  make  it  possible  to  grow  samples 
twice  the  size  of  the  previous  boat.  The  new  boat  design  did  increase  the 
waveguide  performance  and  lower  the  threshold. 

While  this  boat  was  being  built,  a  new  E-series  of  DBR  crystals  was  grown 
in  a  six-well  small  boat.  Previous  runs.  Cl  through  C6,  had  indicated  that 
growth  at  800°C  was  not  as  good  as  at  850°C  because  of  composition  phase 
changes  at  800°C.  These  E-series  runs  produced  fairly  good  DBR  samples,  but 
they  had  a  waveguide  thickness  greater  than  0.6  pm,  which  was  later  decided  to 
be  the  optimum  desired  thickness.  Also,  the  active  layer  tended  to  be  too 
thin  causing  thresholds  greater  than  2  kamps/cm  . 

When  the  new  seven-hole  boat  was  ready,  a  new  temperature  profile  had  to 
be  made  of  the  furnace  containing  the  new  boat.  A  profile  was  made  with  a 
gradient  1  °C  colder  at  the  final  melt  hole  than  the  first  hole  and  a  neutral 
top-to-bottom  gradient.  We  did  not  realize  until  the  complete  series  of  crys¬ 
tals  (Table  2-6)  were  grown  that  the  oven  profile  caused  etchback  of  the  just 
grown  epilayer.  The  proper  profile  should  be  flat  or  increase  in  temperature 
about  1 °C  from  the  first  hole  to  the  last  hole  to  prevent  etchback. 

Because  the  epilayers  grown  in  the  first  five  crystals  of  the  H-series 
were  uneven,  a  boat  modification  was  made  to  permit  excess  melt  to  flow  into 
channels  along  corners  of  the  melt  well.  We  thought  this  modification  would 
help  the  layers  grow  flat  and  parallel,  since  the  channels  would  equalize  the 
heights  of  all  the  melts  and  result  in  a  uniform  melt  pressure  from  well  to 
well.  This  modification  did  not  produce  the  hoped  for  improvement  in  layer 
uniformity.  The  typical  layer  variation  for  the  H-series  crystals  is  shown  in 
Fig.  2-11.  The  corresponding  surface  roughness  of  the  same  crystal  is  shown 
in  Fig.  2-12.  Typical  layer  uniformity  and  surface  smoothness  for  LPE-grown 
crystals  are  illustrated  in  Figs.  2-13  and  2-14.  At  present,  we  are  experi¬ 
menting  with  the  old  six-well  boat  and  the  latest  seven-well  boat  with  flat 
and  increasing  temperature  profiles.  The  preliminary  results  indicate  much 
smoother  layer  structure. 

The  H-series  crystals  were  grown  from  a  new  ingot  of  substrate  materials 
that  had  an  etch-pit  density  of  1500/in  .  In  previous  runs  of  A-series 
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Table  2-6.  DBR  Structure  H-Series 
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through  E-series  crystals  substrates  with  an  etch-pit  density  of  1000/in. ^  or 
less  were  used.  The  difference  in  etch  pit  density  did  not  seem  to  adversely 
affect  the  resulting  +  laser  thresholds.  Before  we  discovered  that  the 
furnace  profile  used  for  the  H-series  was  incorrect,  we  checked  the 
orientation  of  the  H-series  substrates  by  x-ray  diffraction.  We  found  that 
the  substrates  were  oriented  about  1  deg  off  the  100  axis,  which  is  very  close 
to  the  optimum  orientation  of  0.8  deg  off  the  100  axis  for  suppression  of  ter- 
race  formation  for  growth  at  850°C.  However,  to  achieve  the  high  optical 
and  electrical  properties  over  the  large  areas  required  to  make  a  high¬ 
brightness  laser,  we  completed  the  design  and  are  presently  assembling  a  metal 
organic  chemical  vapor  deposition  (MOCVD)  system.  The  expected  capabilities 
of  the  MOCVO  system  are  described  in  Section  2.2. 

With  proper  furnace  temperature  profiles  and  control,  a  pure  hydrogen 
atmosphere,  a  well-designed  boat,  and  high-quality  substrates  at  the  proper 
orientation,  we  hope  to  solve  the  problem  of  producing  crystals  with  high- 
quality  optical  and  electrical  properties  for  low-threshold  DRR  lasers. 

We  had  two  crystals  grown  in  the  large  optical  cavity  configuration  by 
Laser  Diode  Laboratories  (LDL).  The  measured  parameters  of  the  LDL  crystals 
are  shown  in  Table  2-7.  As  shown  in  Figs.  2-13  and  2—14,  the  layers  are  even, 
and  the  surface  is  relatively  smooth.  The  optical  quality,  measured  by  wave- 
guiding  through  the  crystal,  was  very  good.  The  Fabry-Perot  (FP)  lasers 
fabricated  from  the  LOT.  material  generally  were  a  factor  of  2  or  more  higher 
than  the  H-series  thresholds.  The  DBR  lasers  that  were  fabricated  during  this 
contract  period  were  made  from  sample  LDL-102. 

2.2  METAL  ORGANIC  CHEMICAL  VAPOR  DEPOSITION  SYSTEM 

The  metal  organic  chemical  vapor  deposition  (MOCVD)  system  being  designed 
and  built  at  Aerospace  consists  of  a  number  of  subsystems.  The  systems 
structure  is  shown  in  block  diagram  form  in  Fig.  2-15.  A  hydrogen  carrier  gas 
is  bubbled  through  the  metal  organic  liquids  trimethyl  gallium,  trinethyl 

^D.L.  Rode,  R.W.  Wagner,  and  N.F,.  Schunaker,  Appl.  Phys.  Lett.  30,  75  (1977). 
’O.L.  Rode,  Phys.  Status  Solidi  J)2,  425  (1975). 


aluminum,  and  the  dopants  diethyl  zinc  and  diethyl  telluride.  These  liquids 
are  all  kept  at  constant  temperatures  in  constant  temperature  baths.  The 
various  streams  of  carrier  gas  now  saturated  with  the  different  metal  organic 
components  are  mixed  with  arsine  gas  and  finally  diluted  with  hydrogen.  The 
mixture  enters  the  reaction  chamber  and  reaction  occurs  upon  a  substrate 
heated  on  a  rotating  pedestal  to  approximately  700°C.  The  exhaust  gas  is  then 
decomposed  in  a  furnace  heated  above  300°C.  The  gases  exiting  this  scrubber 
are  mainly  H7,  traces  of  arsine  and  metal  organics,  and  vapors  of  the  un- 
rcacted  metals.  These  trace  components  are  eliminated  with  a  small  in-line 
carbon  scrubber.  The  H7  is  then  burned,  and  the  exhaust  gas  is  sampled  for 
arsine  by  an  arsine  detector.  The  dashed  lines  in  Fig.  2-15  indicate  that 
parts  of  the  system  shown  in  the  figure  are  in  a  fume  hood.  The  parts  of  the 
system  not  shown  that  are  outside  the  hood  are  the  vacuum  system  or  some  part 
of  the  gas  handling  system. 

Figure  2-16  is  a  detailed  presentation  of  the  gas-handling  system  used 
for  MOCVD.  The  important  points  to  be  noted  are:  (1)  The  system  is  entirely 
stainless  steel  and  connections  are  VCR  type  (two  opposing  ring-like  faces 
bite  into  a  stainless-steel  washer).  (2)  All  valves  have  welded  bellows,  and 
air-operated  valves  are  used  in  key  positions  for  safety  and  ease  of  opera¬ 
tion.  (3)  Flow  control  valves  can  either  be  manually  operated  or  remotely 
operated  using  a  computer.  The  system  has  been  designed  for  safety,  for  ver¬ 
satility,  for  system  integrity  needed  to  do  high-quality  semi conduc tor  growth, 
and  for  the  experimental  control  needed  for  reproducibility  and  composition 
control.  Future'  needs  will  hi'  taken  into  account  in  design  of  the  system 
(Table  2-8). 

It  is  our  belief  that  the  MOCVD  system  will  permit  the  growth  of  uniform 
layers  over  large  (0.5  by  0.5  cm)  areas  with  high  electrical  and  optical 
quality  because  of  the  excellent  growth  characteristics  of  MOCVD  material  as 
compared  to  liquid-phase  grown  material.  Table  2-9  provides  a  comparison  of 
the  growth  characteristics  of  the  two  techniques.  MOCVD  appears  superior  to 
1,1’F  growth  in  all  respects.  MOCVD  is  more  reproducible  than  l.!n'.  because  MOCVD 
growth  does  not  depend  as  critically  on  temperature  for  the  growth  rate  or  the 
prior  history  of  a  melt  for  layer  composition;  it  is  easier  to  grow  a  layer  to 
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-16.  Gas-Handling  System  for  MOCVD 


Table  2-8.  Comparison  of  Capabilities  of 
MOCVD  and  LPE  Systems 


Parameter 

MOCVD 

System 

LPE 

Different  Material 

Systems  (same  apparatus) 

Yes 

No 

Controlled  Vertical 

Layer  Composition 

Yes 

No 

Artificial  Crystal 

Lattices 

Yes 

Virtually 

impossible 

Minimum  Layer 

Thickness,  A 

50 

500 

Turnaround 

1  day 

1  week 

Table  2-9.  Comparison 
MOCVD  and 

of  Characteristics 
LPE  Systems 

of  Growth  by 

System 

Parameter 

MOCVD 

LPE 

Reproducibility 

Excellent3 

Poorb 

Layer  Composition 
Homogeneity 

Excellent3 

Fair  to  goodb 

Layer  Thickness 
Variations 

Minima la 

Intrinsic 

variations'5 

Layer  Composition 
Accuracy 

Excellent 

Fair  to  good 

Final  Surface 

Morphology 

Smooth 

Rough 

aNonequi librium 
^Unstable  equilibrium 
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a  specific  thickness  and  composition  and  repeat  the  growth  of  these  layers 
with  MOCVD.  LPE  layers  are  not  homogeneous  in  composition  from  top  to  bottom 
because  the  composition  of  the  melt  is  constantly  changing  during  deposi¬ 
tion.  This  is  not  a  problem  for  MOCVD  where  gas  is  constantly  entering  the 
chamber,  flowing,  by  the  substrate,  and  exiting  from  the  chamber.  In  this 

I 

case,  a  steady-state  condition  is  set  up  and  growth,  if  uninterrupted,  remains 
unchanged  with  time.  There  is  a  minimal  variation  in  thickness  across  the 
MOCVD  grown  layers  because  the  rate  of  growth  is  not  a  sensitive  function  of 
temperature,  and  any  temperature  fluctuations  across  the  heating  stage  do  not 
show  up  as  thickness  variations.  This  is  a  very  serious  problem  for  LPE  grown 
material  because  LPE  is  so  sensitive  to  temperature.  Desired  layer  composi¬ 
tion  can  be  obtained  with  MOCVD,  but  it  is  harder  to  obtain  with  LPE  growth. 
The  growth  of  GaAlAs  by  LPE  involves  the  addition  of  aluminum  to  the  melt  in 
situ.  It  is  not  necessarily  true  that  all  the  aluminum  will  dissolve  in  the 
melt  nor  that  the  solution  will  be  uniform.  The  resultant  composition  of 
the  layer  can  vary  greatly  from  the  desired  composition.  This  error  is  ampli¬ 
fied  by  the  fact  that  a  small  error  in  the  concentration  of  aluminum  in  the 
melt  causes  a  large  error  in  the  amount  of  aluminum  in  the  GaAlAs  layers.  The 
morphology  of  the  top  face  and  interfaces  between  the  layers  is  poorer  for  LPE 
grown  material  than  MOCVD  material.  The  top  face  has  poor  morphology  because 
there  is  no  way  to  remove  the  excess  liquid  on  the  sample  before  removing  the 
sample  from  the  growth  chamber.  The  interfaces  have  poor  morphology  because 
the  underlying  layers  are  subjected  to  a  nonequilibrium  environment,  and  there 
is  some  solution  and  dissolution  of  material  from  the  underlying  layer  when 
the  new  layer  is  initially  grown.  No  such  problem  occurs  with  MOCVD  because 
it  occurs  so  far  from  equilibrium  that  the  only  two  processes  that  can  occur 
are  deposition  or  no  deposition  of  a  new  layer;  evaporation  of  a  layer  does 
not  occur  at  normal  working  temperatures. 

MOCVD  is  a  more  versatile  experimental  tool  than  LPE.  With  MOCVD,  mate¬ 
rials  that  have  not  been  duplicated  with  LPE  are  synthesized.  One  can,  in 
principle,  make  artificial  lattices  of  any  complexity  with  individual  layer 
thickness  to  50  A.  With  electrical  control  of  the  flow  controllers,  it  is 
possible  to  control  the  vertical  composition  or  doping  of  a  layer.  None  of 
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these  controls  on  composition,  doping,  or  thickness  can  be  matched  with  the 
current  LPE  system. 


The  MOCVD  project  started  in  late  1978.  At  this  time,  we  have  all  of  the 
critically  needed  components  that  are  necessary  for  a  working  system.  The 
laboratory  that  is  to  house  the  experiment  is  near  completion,  and  we  are  in 
the  process  of  assembling  the  system.  The  system  will  take  another  three 
months  for  completion,  but  will  probably  not  run  until  the  early  part  of  the 
summer  because  of  the  additional  equipment  that  has  to  be  installed  to  meet 
safety  standards. 

The  experiment  that  we  are  setting  up  is  a  very  complex  and  ambitious 
project.  Instead  of  starting  with  a  simple  inflexible  system,  we  have 
designed  and  are  building  a  very  sophisticated  and  very  flexible  system  that 
we  will  be  able  to  modify  in  the  future  to  meet  our  needs. 
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3.  MATERIAL  CHARACTERIZATION 


3. 1  LAYER  THICKNESS  AND  COMPOSITION 

In  order  Co  predict  and  evaluate  both  device  and  crystal  growing 
performance,  the  epilayers  must  be  fully  characterized  with  respect  to  layer 
thicknesses  and  compositions.  Layer  thicknesses  are  determined  by  means  of 
both  optical  and  electron  microscopy.  Figure  3-1  is  an  SEM  photograph  of  a 
cleaved  edge  of  sanple  H5. 

Three  techniques  have  been  investigated  for  obtaining  layer  composi¬ 
tions:  (l)  photoluminescence;  (2)  electron  dispersive  analysis  by  x  rays 

(EDAX);  and  (3)  ion  microprobe  mass  analysis  (IMMA).  Layer  compositions  of  a 
k-deg-angle  lapped  specimen  were  obtained  with  the  use  of  all  of  these 
techniques. 

Because  it  yields  real-time  data  and  does  not  require  the  use  of  a 
correction  factor,  the  IMMA  technique  has  been  deemed  the  most  accurate  and 
useful  of  the  three.  The  photoluminescence  data,  while  not  inconsistent  with 
the  IMMA  data,  must  be  corrected  for  wavelength  shifts  because  of  dopants. 
This  effect  leads  to  an  underestimation  of  the  aluminum  content  of  a  given 
layer.  In  addition,  the  photoluminescence  data  must  be  taken  on  a  point-by¬ 
point  basis  and  is  not,  therefore,  real  tine.  The  EDAX  technique  suffers  from 
similar  drawbacks.  While  not  influenced  by  dopants,  the  I.  emission  line  of 

arsenic  interferes  with  both  the  K  aluminum  line  and  the  K  gallium  line. 

a  a 

Stripping  the  arsenic  spectrum  mathematically  from  the  composite  spectrum 
alleviates  this  problem  to  some  degree  but  still  yields  what  appears  to  be  a 
consistently  low  aluminum  concentration.  Finally,  as  in  the  case  of  the 
photolumi nescence  technique,  data  must  be  gathered  on  a  point-by-point  basis 
and  is  not  real  time. 

Composition  traces  obtained  by  means  of  the  IMMA  technique  (sample  H5) 
are  shown  in  Fig.  3-2.  This  technique  has  been  used  on  a  routine  basis  to 
determine  AlAs  concentrations. 
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EXTRA  LAYER  7.4  fim 
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Fig.  3-1.  SEM  Photograph  of  Cleaved  Edge  of  Sample  H5 
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Fig.  3-2.  IMMA  Scan  for  Sample  H5 
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3.2  WAVEGUIDE  EVALUATION 


Experiments  were  begun  on  the  guiding  of  1.15-pra  light  in  the  GaAlAs  and 
GaAs  waveguides  to  aid  in  the  development  of  the  beam  expander  and  to  evaluate 
the  material  quality  for  DBR  lasers.  In  the  early  samples,  it  was  discovered 
that  the  optical  quality  of  the  GaAlAs  layers  grown  by  LPE  was  of  poor  qual¬ 
ity.  In  Fig.  3-3,  near-field  photographs  of  light  guided  through  various  epl- 
layers  are  shown.  The  vertical  light  is  the  guided  light;  the  haze  left  of 
the  line  is  light  going  through  the  substrate;  the  haze  right  of  the  line  is 
light  going  through  the  air.  In  Fig.  3-3(a),  the  light  is  shown  guided 
through  a  piece  of  vapor  phase  epitaxial  GaAs,  which  was  the  same  material  as 
used  for  the  passive  studies  of  beam  expansion.  The  guided  light  forms  a 
smooth  line,  corresponding  to  good-quality  epitaxial  material. 

In  contrast,  passive  guiding  experiments  with  sample  A19  indicate  that 
there  is  a  tendency  for  the  light  to  break  up  as  it  travels  in  the  waveguide, 
causing  a  characteristic  beading  pattern  shown  in  Fig.  3-3(b). 

A  program  was  therefore  begun  to  improve  the  quality  of  the  epilayers 
grown  by  LPF,.  It  was  discovered  that  if  a  GaAs  buffer  epilayer  is  grown 
first,  that  is,  before  the  layers  of  GaAlAs,  the  quality  of  the  waveguides  is 
improved.  In  Fig.  3-3(c),  light  has  been  guided  through  sample  C2,  an  epi¬ 
layer  that  includes  a  buffer  GaAs  layer.  There  is  almost  no  beading,  and  the 
quality  of  the  optical  waveguide  is  almost  as  good  as  that  of  the  vapor-phase 
epitaxy  (VPE).  The  structure  of  sample  C2  is  shown  in  Fig.  3-4,  including  the 
GaAs  buffer  layer.  This  huffer  layer  was  included  in  all  later  DBR  laser 
structures  to  improve  their  efficiency. 

The  effectiveness  of  the  GaAs  buffer  layer  is  the  result  of  the  poor 
optical  quality  of  the  GaAlAs  layers  caused  by  nonunifornities  in  the  re¬ 
fractive  index  that  are  caused  by  the  nonuniform  growth  of  GaAlAs  around  dis¬ 
locations.  It  is  a  well-known  fact  that  an  epilayer  of  GaAs  has  fewer  dislo¬ 
cations  than  the  substrate  upon  which  it  is  grown.  By  contrast,  an  epilayer 
of  GaMAs  tends  to  increase  the  number  of  dislocations,  since  there  is  a 
slight  lattice  mismatch.  Therefore,  the  buffer  layer  is  required  to  remove  as 
many  dislocations  as  possible  before  the  wavoguiding  GaAlAs  layers  are  grown. 
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Fig.  3-3.  Passive  Waveguiding  Through  (a)  Vapor-Phase  Epitaxy  of  GaAs, 
(b)  ALOC  Material,  and  (c)  GaAlAs  LPE  with  Buffer  Layer 
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3. 3  PHOTOLUMINESCKNCE  MEASUREMENTS 


In  the  process  of  developing,  etch  techniques  for  the  selective  renoval  of 
epilavers  from  the  Al.hC  structure,  we  found  it  important  to  monitor  whether  or 
not  the  laver  had  indeed  been  removed.  Ue  developed  a  photoluninesc.ence  tech¬ 
nique  to  make  this  evaluation. 

\  diagram  of  the  experimental  apparatus  is  shown  in  Fig.  3-5.  The  peak 
oc  the  phot o l tin i noseenoe  spectrum  is  a  function  of  the  aluminum  concentration 
and  the  doping  level.  In  Fig.  3-b,  a  nap  of  the  photoluninescence  spectrum  of 
a  2-dog  angle  lapped  section  of  crystal  186  is  shown.  The  photolmninescence 
peaks  corresponding,  to  the  cap,  p-isolat i on,  active,  waveguide,  n-isolation, 
buffer,  and  substrate  lavers  are  indicated.  In  Fig.  3-7(a),  a  photolumi- 
losciMH'e  scan  of  sample  lSbE  before  processing  is  shown.  After  20  min  of  ion 
milling,,  the  scan  shown  in  Fig.  3 — 7  ( b ")  indicates  the  cap  layer  has  been  re- 
•  noved  tod  the  p-isolat ion  layer  is  exposed.  After  10  min  of  exposure  to  AB 
etch.  Fig.  i-7(c)  indicates  that  the  cap  layer  has  been  removed,  and  two 
peaks,  corresponding  to  the  ictive  and  waveguide  layer,  are  seen.  After  an 
additional  5  min  exposure  to  an  AS  etch,  the  active  layer  is  removed,  with 
'll:  'he  waveguide  layer  exposed,  Fig.  3-7(d). 

"rhis  photo] uni nescenco  technique  is  used  extensively  in  fabricating  D3R 
1  isers  anil  also  provides  additional  information  on  the  quality  of  the  crystal. 

1 . 4  -  AHHiCA  ''  I '  IN  (>F  H  R' )  A 1 1  —  \hK  \  F.VbiY-PFBOT  l.ASERS 

\s  a  test  of  the  electrical  and  lasing  qualities  of  the  grown  material, 
broad-area  Fuhry-I’erot  (FF)  lasers  were  rahrieated,  and  their  lasing  spectra 
and  threshold  ch a r  ie t er i s t i c s  were  obtained. 

Figure  3-b  is  a  flow  chart  showing  tile  sequence  in  which  the  material 
chn rae ter i  -  i  ng  steps  ar.-  performed.  The  results  of  the  characterization  steps 
for  the  K,  H,  and  12  V.  series  of  LUF  grow  >  material  are  summarized  in  Tables 
2-5,  2-6,  and  2-7,  respectively. 
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FiR.  3-5.  Experimental  Apparatus  for  Photoluminescence  Measurements 
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Fig.  3-6.  Photoluminescence  Map  of  Sample  LDL  186.  Dots  indicate  peak  of 

photoluminescence;  vertical  lines  indicate  the  relative  intensitv 
of  the  photoluninescence  within  the  epilaver. 
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RELATIVE  INTENSITY 


t'i;;.  1-7.  Photo turn  nose once  Scan  With  (a)  Cap  Layer  Intact,  (b)  Cap  Layer 
Removed  and  p- Isolation  Layer  exposed,  (c)  p-Tsolation  Layer 
Removed  and  Active  Layer  Exposed,  and  (d)  Active  Layer  Removed 
and  Waveguide  Layer  Exposed. 


3-10 


\ 


EPILAYER 

STRUCTURE 


Fig;.  3-8.  Flow  Chart  for  Characterizing  Material 
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4.  FABRICATION  TECHNIQUES  AND  EXPERIMENTAL  RESULTS 

4. 1  DBR  LASERS 

4.1.1  Fabrication 

The  development  of  the  fabrication  techniques  has  been  described  in  a 
series  of  monthly  progress  reports.  The  fabrication  processes  that  are 
described  here  are  a  result  of  refining  many  different  approaches.  A  de¬ 
tailed,  step-by-step  process  schedule  is  included  in  Appendix  C. 

Figure  4-1  is  a  cross  section  through  an  asymmetric,  large-optical 
cavity  DBR  laser.  In  the  DBR  laser,  feedback  is  provided  by  a  cleave  on  one 
end  and  a  distributed  Bragg  reflector  (DBR)  etched  into  the  GaAlAs  waveguide- 
Si02  interface  on  the  other  end  of  the  active  layer.  The  DBR  end  is  sawed  to 
eliminate  end  reflections. 

The  active-passive  laser  is  identical  to  the  DBR  laser  in  construction, 
except  that  the  DBR  grating  is  eliminated  and  feedback  is  provided  at  both 
ends  by  cleaves.  A  primary  reason  for  fabrication  of  the  active  passive  de¬ 
vice  is  to  evaluate  the  coupling  between  the  active-passive  interface. 

DBR  lasers  and  active-passive  lasers  are  fabricated  in  the  same 
manner.  Standard  photolithographic  procedures  are  used;  0.5-mm-wide  strips  of 
AZ1375  Shipley  photoresist  spaced  1  mm  apart  and  3  Pm  thick  are  developed  on 
an  asymmetric  large  optical  cavity  (LOC)  crystal,  Fig.  4-2(a).  The  sizes  of 
the  crystals  that  have  been  processed  range  from  3  by  4  mm  to  4  by  12  mm. 

The  photoresist  stripes  protect  the  underlying  GaAs  from  erosion  when 
the  sample  is  subjected  to  chemical  etchants  and  ion-beam  etching.  A 
combination  of  chemical  and  ion-beam  etching  is  used  to  remove  the  first  three 
layers  (cap,  p-isolation,  and  active)  between  the  photoresist  stripes.  After 
the  photoresist  is  removed,  mesa  structures  composed  of  the  original  six 
epilayers  remain  separated  by  valleys  of  three  epilayers  (waveguide,  n- 
isolation,  and  buffer). 
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Fig.  4-1.  Cross  Section  Through  Asymmetric  LOC  DBR  Laser 
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2.  Sample  Before  Processing.  (a)  Top  view.  (b)  Side  view  aft 
removal  of  three  livers  by  ion-bean  and  chemical  etching. 


If  DBR  lasers  are  being  fabricated,  a  holographically  formed  grating 
with  a  period  of  about  3600  A  is  ion-beam  etched  into  the  top  surface  of  the 
crystal.  The  grating  structure  over  the  waveguide  layer  provides  third-order 
Bragg  scattering  and  forms  the  DBR.  The  grating  that  is  formed  on  the  mesa 
region  has  no  effect  on  the  device  performance.  Since  the  procedure  for  fab¬ 
ricating  gratings  for  DBR  lasers  is  common  to  making  gratings  for  distributed 
Bragg  deflectors  (DBD)  and  output  couplers,  a  discussion  of  grating  fabrica¬ 
tion  is  given  in  the  next  section. 

Once  the  grating  is  in  place  (for  DBR  lasers),  a  3000-A-thick  layer  of 
SiO>  is  deposited  over  the  entire  surface  of  the  sample.  A  1 . 6-jjm-thick  layer 
of  Shipley  AZ1350J  photoresist  is  deposited  over  the  oxide,  and  a  series  of 
windows  (25  by  450  pm)  spaced  250  pm  apart  are  developed  through  the  photo¬ 
resist  down  to  the  oxide  layer  on  the  mesa  regions.  The  oxide  layer  exposed 

by  the  windows  is  etched  away  reactively  in  a  plasma  gas  of  CF^.  Windows 
exist  through  the  photoresist  and  oxide  layer  down  to  the  GaAs  cap  layer  on 
the  mesas.  A  layer  of  platinum  and  gold  contacts  can  now  be  electroplated  to 
the  CaAs  in  the  windows  through  the  photoresist  and  the  oxide  layer  (Fig. 
4-3).  The  photoresist  and  oxide,  layers  provide  a  mask  to  confine  the  plating 
to  the  conductive  GaAs  stripe  and  also  provide  sidewall  confinement  of  the 
plated  material.  The  back  side  of  the  sample  is  also  electroplated  to  form 
the  bottom  contact.  A  detailed  explanation  of  the  plating  procedure  is  given 
in  Appendix  D.  After  cleaning  the  sample,  the  crystal  is  cut  into  individual 

lasers.  The  bottom  side  of  the  laser  is  soldered  to  a  TO-18  transistor 

header,  and  a  gold  wire  is  ball  bonded  to  the  top  contact  and  to  a  header 
post.  A  complete  DBR  laser  is  shown  in  Fig.  4-4. 

Preliminary  tests  of  the  laser  diodes  include  checking  the  I-V 
characteristics  on  a  transistor  curve  tracer  and  observing  the  luminescence 
through  a  hand-held  infrared  viewer. 

4.1.2  Hxperimental  Results  for  DBR  Lasers 

The  mode  spectrum,  total  power,  light  current  characteristics  and  beam 
divergence  '•ere  measured  for  the  fabricated  DBR  lasers.  Figure  4-5  is  a  block 
diagram  of  the  laser  test  setup.  A  typical  spectrum  of  a  DBR  laser  is  shown 


Fig;.  4-5.  Experimental  Setup  for  Laser  Spectrum  Measurements 
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in  Fig.  4-6.  In  order  to  verify  that  the  laser  action  occurred  from  the 
grating  and  not  from  spurious  reflections,  several  tests  were  performed. 


It  was  observed  that  the  DBR  spectrum  was  a  single  mode.  Comparable 
spectra  for  Fabry-Perot  lasers  from  the  same  material  were  multimode  (Fig. 
4-7).  In  addition,  a  DBR  laser,  which  operated  at  room  temperature,  was 

cooled  to  liquid  nitrogen  temperature,  and  laser  action  ceased,  since  the 
reflectivity  of  the  grating  was  no  longer  at  the  frequency  of  the  fluo¬ 
rescence.  This  is  in  contrast  to  Fabry-Perot  lasers  whose  frequency  shifts 
with  temperature. 

In  some  of  the  DBR  lasers,  two  lines  occurred  (Fig.  4-8)  because  the 
waveguide  was  sufficiently  thick  that  not  only  the  lowest  even  mode  but  also 
the  lowest  odd  mode  could  propagate. 

An  example  of  a  two-mode  DBR  laser  is  shown  in  Fig.  4-9.  As  shown  in 
Fig.  4-10,  a  grating-wave  vector  K  (K  e  2-ti/a)  can  cause  a  DBR  laser  to  oscil¬ 
late  in  more  than  one  transverse  mode,  depending  upon  the  spacing  of  the 

waveguide  modes  to  the  width  AX  of  the  gain  profile  and  the  relative  strength 
of  the  intermode  and  intramode  coupling  coefficients.  A  graphical  solution 
for  the  DBR  lasing  spectra  that  accounts  for  the  refractive  index  variation 
with  wavelength  is  indicated  in  Fig.  4-11,  The  graphical  solution  compares 
favorably  with  the  40-A  spacing  between  modes  observed  experimentally. 

The  linewidths  just  above  threshold  were  typically  1  A.  DBR  laser  G8, 
(Fig.  4-6)  has  a  linewidth  of  1.1  A,  and  DBR  laser  G9  (Fig.  4-8)  has  a  width 
of  1 . 6  A.  The  apparent  linewidth  broadened  slightly  with  pumping  (Fig.  4-12) 
in  which  DBR  laser  G5  is  pumped  at  three  times  threshold.  The  line  had 

broadened  to  about  3  A,  but  it  was  observed  on  the  oscilloscope  that  the  fre¬ 

quency  changed  during  the  100-nsec  input  pulse.  The  broadening  is  believed  to 
be  a  chirp  caused  by  heating  in  the  junction.  Such  effects  have  been  observed 
in  Fabry-Perot  diode  lasers. 

Diode  G9  had  a  complicated  behavior,  since  the  two  frequencies  corre¬ 
sponding  to  two  guided  modes  occurred  at  high  pumpir.g  rates.  It  was  observed 
that  at  low  excitation  levels,  the  shorter  wavelength  was  the  strongest  peak, 
and  the  longer  wavelength  occurred  only  during  the  initial  portion  of  the 
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Fig.  4-9,  Two-Mode  DBR  Laser 
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4-10.  (o-B  01  a;; ram  of  Three  Scattering  Processes.  The  {’rating  wave 
vector  K  =  2tt/A  has  same  length  for  each  process.  Backward 
scat  ter  in;;  between  modes  1-1,  1-2,  and  2-2  occurs  at 
frequencies  <o  to  ( ,  and  to.,,  respectively. 
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Fig.  4-12.  DBR  Laser  Spectra.  (a)  Diode  G5  at  three  times  threshold 
current.  (b)  Diode  G9  just  above  threshold. 
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pulse.  Evidently,  heating  caused  the  resonance  condition  to  change.  Near  the 
end  of  the  pulse,  the  shorter  wavelength  had  a  lower  threshold  and  was  the 
only  frequency  above  threshold.  These  data  demonstrate  that,  even  when  a 
grating  is  used,  it  is  possible  to  get  several  wavelengths  oscillating,  and 
not  necessarily  simultaneously. 

The  complexity  of  device  G9  can  be  seen  from  the  light-current  curve 
(Fig.  4-13),  which  has  several  kinks.  The  DBR  laser  causes  only  one  or  two 
frequencies  to  oscillate,  but  does  not  necessarily  improve  light-current 
characteristics,  in  relation  to  Fabry-Perot  lasers. 

Since  kinks  seem  to  be  related  to  the  filamentary  structure  of  lasers,^ 
and  since  we  were  interested  in  beam  quality,  in  order  to  perform  beam  ex¬ 
pansion,  we  made  a  careful  study  of  the  diffraction  qualities  of  the  light 
emitted  from  the  DBR  lasers.  Depending  on  the  device  used,  we  were  able  to 
observe  diffraction-limited  operation  just  above  threshold.  This  is  dem¬ 
onstrated  by  studying  the  divergence  of  sample  G9,  right  above  threshold. 
Near-field  patterns  were  observed  both  from  the  cleave  and  the  grating  face 
(Fig.  4-14).  In  addition,  the  in-plane  beam  divergence  was  measured  coming 
out  of  the  cleave.  This  was  accomplished  by  monitoring  the  spot  size  as  a 
function  of  distance  from  the  cleave.  The  data  are  shown  in  Fig.  4-15. 

The  data  are  analyzed  by  comparing  the  measured  divergence  with  that 
expected  for  a  Gaussian  beam  of  the  width  determined  from  the  near  field 
intensity  profiles  (Fig.  4-14). 

The  angular  divergence  of  a  Gaussian  beam  is 


0  =  X/irw  (4-1 ) 

o  o 


where  0  is  the  half-angle,  and  w  is  the  half  width  of  the  beam  at  l/e^  in 
o  o 

intensity.  Redefining  the  quantities  in  terms  of  full  width  half  maximum, 


e 


\ 

0.  72tio) 


(4-2) 


N.  Chinone,  J.  Appl. 


Phys.  48,  3237  (1977). 
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Fig.  4-13.  Light-Current  Characteristics  for  Diode  G9 
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Fig.  4-14.  Near-Field  Profiles  for  DBR  Laser  with  Diode  G9  at  Threshold 
(a)  Cleaved  face.  (b)  Grating  face. 
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The  filament  measured  in  the  near  field  had  an  apparent  width  of 
6.4  Pm(correct ion  has  not  been  made  for  resolution  of  the  microscope  objective 
used  in  this  measurement).  The  far-field  beam  divergence  out  the  cleave  was 
0.072  rad,  which  corresponds  to  the  diffraction  expected  from  a  filament  of 
width  5.4  Pm.  This  agrees  well  with  the  measured  filament  width  and  supports 
the  assumption  that  this  laser  operated  in  a  single  spatial  mode. 

The  laser  output  from  the  cleaved  face  was  observed  to  have  a 
near-field  profile  that  was  a  single  filament  22  pm  wide.  This  profile  is  the 
beam  that  has  diffracted  from  the  5.4-pm  filament  in  the  laser  region  as  a 
result  of  a  1-mm  grating  (Fig.  4-16).  If  the  filament  is  assumed  to  extend 
the  length  of  the  active  region  of  the  laser  diode  and  begins  to  expand  when 
it  enters  the  grating,  the  diffraction  should  occur  for  a  distance  of  1  mm,  at 
an  angle  given  by  the  far-field  divergence  angle  divided  by  the  index  of 
refraction  of  the  medium  (0.02  rad).  This  divergence  would  cause  the  output 
of  the  grating  to  have  a  width  of  20  pm,  which  compares  favorably  with  the 
measured  filament  width  of  22  Pm. 

For  completeness,  out-of-plane  beam  divergence  was  measured  also.  It 
was  observed  that  the  frequency  corresponding  to  the  higher  order  waveguide 
mode  caused  a  double-humped  far-field  pattern  from  the  cleave.  The  divergence 
was  greater  than  0.48  rad  and  could  not  be  measured  since  it  was  limited  by 
the  aperture  of  the  lens  that  was  used.  The  other  frequency  was  a  lowest 
order  waveguide  mode  and  had  a  beam  divergence  of  0.3  rad.  This  supports  the 
assumption  discussed  earlier  that  the  two  frequencies  are  caused  by  the 
different  modes  occurring  In  the  laser  region. 

Other  diodes  did  not  necessarily  show  in-place  diffraction-limited 
behavior.  The  near-field  profiles  of  diode  G5  are  shown  in  Fig.  4-17.  Al¬ 
though  the  apparent  near-field  profile  was  15  Pm  wide,  the  beam  divergence  was 
0.165  rad,  corresponding  to  a  filament  size  of  2.4  pm.  In  other  words,  the 
diode  had  a  highly  filamentary  character,  which  caused  a  great  deal  of  beam 
divergence.  This  fact  was  made  clear  by  looking  at  the  profile  out  of  the 
grating.  The  output  had  several  filaments  and  an  overall  width  of  almost 
100  Pm. 
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Fig.  4-16.  Geometry  for  In-Plane  Divergence  of  DBR  Lase 
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**  4_1 7‘  Near-Field  Profiles  for  DBR  User  with  Diode  G5.  (a)  Cleaved 

face.  (b)  Crating  face. 


Additional  measurements  were  made  of  the  front-to-rear  power  ratio  of 
DBR  lasers  of  differing  lengths  to  determine  the  coupling  coefficient  of  the 
grating  and  to  check  for  its  reproducibility.  Two  devices  were  compared  in 
detail;  the  data  are  shown  in  Table  4-1. 

According  to  the  theory,  a  retroreflector  of  length  L  on  resonance  has 
the  following  behavior  for  transmission: 


and  reflection: 


When  KL  »  1 , 


T 


1 

cosh2KL 


R 


slnh2KL 

2 

cosh  KL 


(4-3) 


(4-4) 


T  =  4e_2KL,  R  =  1  -  T  (4-5) 

The  coupling  coefficient  was  determined  from  Eq.  (4-5)  with  the 
assumption  that  the  front-to-back  power  ratio  will  be  the  same  as  the  ratio  of 
the  transmission  of  the  grating  to  the  transmission  of  the  GaAs  cleave  on  the 
front  output  face  (85%).  The  results  give  an  average  coupling  coefficient  of 
K  =  0.0036  um 

The  theoretical  expression  for  the  coupling  coefficient  of  a  grating  is 
given  for  square  teeth  in  Appendix  B.  The  experimental  grating  looked  more 
like  triangular  teeth  than  square  teeth,  however.  The  comparison  between 
square  and  triangular  shapes  may  be  made  by  comparing  the  Fourier  transforms 
for  the  third-order  coefficients  for  square  and  triangular  shapes.  Using  the 
calculations  shown  in  Fig.  4-18,  we  see  that  the  polarization  induced  by  a 
triangular  grating  will  be  two-thirds  that  induced  by  the  step  index.  Thus, 
we  expect  the  grating  coupling  coefficient  to  be 
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Table  4-1. 

Data  Used  to  Determine 
Gratings  in  DBR  Lasers 

Coupling  Coefficient  of 

G9 

C5 

Power-Out  Grating,  mW 

1 

0.5 

Power-Out  Cleave,  mW 

30 

10 

Power  Ratio  (grating 

to  cleave),  T 

0.03 

0.05 

Ln  T/4 

-5.3 

-4.4 

Grating  Length,  pm 

1000 

464 

K 

0.0026 

0.0047 

Lg,  um 


192 


106 


-7 T 


0  IT 


jr  Jrl2  2a 

F  = J  y  (x)  cos  3  x  dx  =  a J  cos  3  x  dx  =  — 

-7T  -7T/2  3 

(a) 


Fig.  4-18.  Third  Harmonic  Calculation  for  Square  and  Triangular-Shaped 
Crating  Teeth 
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where  n  is  the  refractive  index  of  the  guide  and  where  we  have  taken  the 
result  of  Yariv,^  mu^t;iplied  by  two-thirds,  and  used  teff  instead  of  t  for 
guides  that  are  not  far  from  cutoff. 

From  SEM  photographs  of  the  grating,  we  estimate  a  =  A/3,  t  =  1.54  pm, 

and  A  =  0.36.  The  grating  was  used  in  third  order,  and  the  refractive  index 

of  GaAs  at  0.88  Pm  is  3.6.  These  numbers  give  the  theoretically  expected 

coupling  coefficient  of  K  =  0.0044  pm  ^ .  In  order  to  calculate  this  number, 

the  effective  guide  width  was  required.  This  was  determined  by  solving  the 

eigenvalue  equation  for  the  waveguide  to  determine  the  exponential  profile 

outside  of  the  guide  and  to  calculate  teff  =  t  +  1/y.  For  the  DBR  lasers, 

2  2 

the  V  parameter  for  the  waveguide  was  Ack^  d  =  146,  resulting  in  two  modes, 
given  by  the  following  eigenvalue  equations:^ 


tan  u 
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(4-7) 


The  solutions  are  kd  =  1.7  and  3.4,  which  yield  Yd  =  12,  with  an  effective 
thickness  of  1.62  Pm.  This  is  the  value  used  in  Eq.  (4-6)  to  determine  the 
theoretical  coupling  coefficient. 


The  excellent  agreement  between  theory  and  experiment  permits  us  to 
predict  behavior  for  other  waveguides.  Me  are  presently  in  the  process  of 
fabricating  devices  from  material  with  a  waveguide  thickness  of  0.6  Pm.  These 
samples  should  yield  devices  with  coupling  coefficients  as  much  as  eight  times 
larger. 


A.  Variv,  Quantum  Electronics,  Second  ed. ,  Wiley,  New  York,  1975. 
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4.2  GRATING  FABRICATION  BY  HOLOGRAPHIC  EXPOSURE 

Holographic  exposure  is  accomplished  by  exposing  the  substrate  to  two 
interfering  argon-ion  laser  beams.  The  argon  laser  beams  are  spatially  fil¬ 
tered,  expanded,  and  collimated  prior  to  being  split.  The  split  beams  are 
interfered  at  an  angle  6,  necessary  for  a  given  periodicity  of  grating.  The 
relationship  between  0,  the  grating  periodicity  and  laser  wavelength,  is 


A 


1 

2  sin  9 


(4-8) 


where  A  is  the  grating  periodicity,  X  is  the  laser  wavelength  (0.4570  Pm),  and 
0  is  the  angle  of  interference.  In  Rig.  4-19,  the  optical  system  for  holo¬ 
graphic  exposure  is  shown. 

In  order  to  achieve  precise  alignment  of  the  substrate  with  respect  to 
the  two  interfering  argon  laser  beams  (and  therefore  accurate  control  of  the 
grating  orientation),  an  HeNe  laser  is  used.  The  HeNe  laser,  aligned  level 
with  respect  to  the  argon  beam,  is  reflected  off  a  (110)  crystal  plane  of  the 
substrate  and  back  onto  the  front  laser  window,  ensuring  that  the  (100) 

crystal  plane  is  normal  to  the  interfering  argon  beams  and  thereby  accurately 
orienting  the  grating. 

The  substrate  is  mounted  on  a  precision  coaxial  rotational  stage 
accurate  to  0.01  deg  of  arc  so  that  once  the  substrate  has  been  oriented  to 
the  normal  it  can  be  rotated  to  the  required  angle  of  grating  orientation. 
The  substrate  is  then  exposed  (~30  sec  at  a  beam  intensity  of  10  mW/cm^)  and 
developed.  The  sample  is  now  ready  for  ion-milling  replication  of  the 
photoresist  grating  into  the  substrate.  A  rotary  ion-milling  technique  is 

used  with  a  beam  current  of  0.5  mA/cm2  and  an  acceleration  voltage  of  500  Vdc. 

A  grating  with  a  period  of  ~3600  A  fabricated  by  this  technique  is 
illustrated  n  Fig.  4-20.  Although  the  techniques  for  grating  fabrication 

were  well  developed  before  the  beginning  of  the  contract  once  etched  struc¬ 

tures  were  fabricated,  it  was  found  that  the  photoresist  deposition  techniques 
required  for  grating  exposure  were  not  sufficiently  smooth.  A  new  technique 
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Fig.  4-20.  Cross  Section  Through  Ion-Milled  DBR  Reflector. 
Grating  periodicity  is  3600  A. 
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was  therefore  developed  in  which  the  application  of  thin  films  of  photoresist 
for  grating  fabrication  was  accomplished  by  a  technique  other  than  centrifugal 
spin  coat  application.  The  spin-coat  technique  was  acceptable  for  a  flat 
sample  with  no  changes  in  surface  topography.  Once  surface  steps  were 
introduced  by  the  etch  techniques,  bordering  effects  occurred,  which  origi¬ 
nated  from  the  bottom  edge  of  the  mesa,  anywhere  from  5  pm  to  0.3  or  0.4  mm, 
depending  on  photoresist  viscosities  and  speed  of  spinning.  Studies  were 
therefore  initiated  to  develop  alternative  techniques  of  photoresist 
deposition. 

Several  techniques  have  been  used  to  apply  photoresist  films, ^  such  as 
dip  coating,  which  is  simply  dipping  the  substrate  into  the  photoresist  solu¬ 
tion  and  drawing  it  out  at  a  fixed  rate  for  a  selected  thickness.  By  means  of 
this  technique  layers  are  produced  whose  thicknesses  vary  from  5  to  50  pm. 
Another  technique  is  roller  coating,  which  involves  the  application  of 
photoresist  by  the  contact  of  a  photoresist  impregnated  roller  against  the 
substrate.  By  means  of  this  technique  layers  are  produced  whose  thicknesses 
vary  from  25  to  300  pm.  One  final  technique  is  a  spray  coat  application, 
which  is  simply  passing  the  substrate  through  a  heavy  spray  of  photoresist  at 
a  fixed  rate.  Typical  layer  thicknesses  for  the  technique  are  3  to  25  pm. 

All  of  the  aforementioned  techniques  result  in  somewhat  thick,  films, 
ranging  from  a  thickness  of  3  pm  for  the  spray  coating  technique  to  0.4  mm  for 
the  roller  technique.  The  film  thickness  necessary  for  holographic  exposures 
of  gratings  is  — 0 .15  pm,  a  factor  of  20  thinner  than  obtainable  with  the 
spray  coating  technique.  However,  since  the  spray  coating  technique  proved  to 
be  the  most  promising,  work  commenced  to  refine  and  modify  the  technique  to 
facilitate  films  as  thin  as  0.1  pm. 

The  problems  to  be  addressed  were  the  proper  atomization  of  the  photo¬ 
resist  and  the  ability  to  precisely  control  the  time  that  the  substrate  was 


^W. S.  Deforest,  Photoresists:  Materials  and  Processes,  McGraw-Hill,  New  York, 
197  5.  ~  . .  . 
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exposed  to  the  spray.  The  first  problem  was  solved  by  determining  the  optimum 
viscosity  of  the  photoresist;  the  selected  viscosity  in  simple  measures  was 
4:1  Shipley  Thinner  to  Shipley  to  AZ1350  B,  respectively. 


Two  other  important  factors  in  the  atomization  process  are  the  proper 
air  to  photoresist  mixture,  which  was  ~100:1,  and  the  supply  pressure,  which 
was  ~80  psi. 

The  second  problem  was  controlling  the  amount  of  deposition  time  for  a 
given  film  thickness.  The  standard  technique  is  the  movement  of  the  substrate 
in  front  of  the  spray.  This  technique  is  unacceptable  because  of  the  mechan- 
ics  involved  with  moving  the  substrate  at  a  high  rate  of  speed  in  order  to 
obtain  a  thin  film.  Two  techniques  were  evaluated  to  accomplish  a  short  de¬ 
position  time,  one  of  which  was  the  switching  of  the  air  pressure  used  to 
propel  the  photoresist.  This  technique  proved  to  be  too  slow  because  of 
residual  pressure  in  the  system  after  the  supply  pressure  had  been  switched 
off.  The  other  technique,  which  was  successful,  was  an  electromechanical 
shutter  placed  between  the  spray  and  the  substrate. 

Because  of  a  lack  of  such  a  device,  one  had  to  be  designed  and  built. 
The  shutter  had  to  be  capable  of  actuation  times  of  less  than  30  msec,  while 
being  subjected  to  a  flow  at  ~80  psi.  A  standard  electromechanical  shutter 
actuator  was  used  to  switch  a  relay  supplying  the  power  to  the  shutter. 

One  last  parameter  in  the  deposition  process  was  the  distance  between 
the  spray  source  and  substrate.  This  variable  was  solved  by  several  test 
depositions.  The  optimum  distance  for  the  most  uniform  coverage  in  this 
system  was  1 1  cm. 

This  photoresist  deposition  system  proved  successful  in  the  uniform 
application  of  thin  films  of  photoresist  on  large  nonflat  surfaces  for  the 
purpose  of  high-resolution  photolithography  and  holography. 

In  parallel  with  the  effort  to  fabricate  DBR  lasers,  we  developed  the 
techniques  for  multiple-grating  exposures  that  would  be  needed  for  fabrication 
of  the  high-brightness  laser. 
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Device  integration  on  a  common  substrate  generates  severe  placement  and 
orientation  tolerances  on  each  of  the  individual  gratings.  In  addition,  each 
grating  must  bear  a  different  periodicity,  thus  requiring  multiple  holographic 
exposure  and  development  steps.  As  a  result  of  these  constraints,  we  must 
have  the  ability  to  control  the  placement,  orientation,  periodicity,  and  the 
amount  of  area  to  be  covered  by  each  individual  device  grating.  The  three 
primary  devices  and  their  relative  area,  placement,  and  orientation  on  the 
common  substrate  are  shown  in  Fig.  4-21. 

A  very  important  process  in  the  multiple-placement  procedure  is 
multiple  masking  for  control  of  the  area  and  placement  of  the  0ratings.  This 
is  photolithographically  facilitated  with  the  use  of  a  novel  photoresist 
mask.  The  first  step  in  the  photolithographic  process  is  the  application  of 
the  primary  photoresist  layer  (~  0.05  Pm  thick)  whose  application  can  be 
accomplished  either  hy  centrifugal  spin  or  by  spray  deposition.  Following 
photoresist  application,  a  window  is  exposed  and  developed  in  the  photoresist 
photolithographically,  Fig.  4-22(a).  At  this  point,  the  photoresist  is 
hardened  at  210°C  for  1  hr.  The  photoresist  mask  is  now  insoluble  in  CH3OH 
and  C2HCI3,  facilitating  a  suitable  surface  for  the  holographic  processing. 

The  next  process  in  the  procedure  is  the  holographic  production  of 
photo-resist  gratings.  The  preholographic ,  surface  preparation  consists  of  an 
abrasive  application  of  CH3OH  followed  by  a  flowing  of  C2HCI3  over  the  sur¬ 
face,  immediately  followed  by  a  drying  flow  of  dry  N2.  The  holographic  photo¬ 
resist  is  then  applied,  using  either  the  centrifugal  spin  or  spray  deposition 
technique,  to  a  thickness  of  0.15  Urn,  Fig.  4-22(b).  Following  the  photoresist 
application  and  prior  to  exposure,  the  coated  sample  is  baked  at  70°C  for  25 
min.  The  photoresist  is  then  holographically  exposed,  Fig.  4-22(c),  and  the 
resultant  grating  replicated  into  the  GaAs  material  by  ion  beam  milling,  Fig. 
4-22(d).  Following  this  step,  the  residual  photoresist  from  the  first  ion 
mask  is  now  removed  in  an  02  plasma.  Fig.  4-22(e). 

The  second  phase  of  the  multiple-masking  procedure  involves  the 
addition  of  succeeding  gratings  on  the  same  substrate.  A  0. 5-Pm-thick  film  of 
photoresist  (for  the  protection  of  the  first  grating  during  the  ion  milling  of 
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Fig.  4  21.  Three  Primary  Devices  on  Common  Substrate 
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A  window  is  then 


the  second  grating)  is  first  deposited  onto  the  substrate, 
exposed  photolithographically  and  developed  in  the  photoresist  film,  Fig. 
4-23(a),  thereby  facilitating  the  application  of  a  0. 13-pm  film  of  photoresist 
on  the  substrate  surface  while  still  providing  an  ion-milling  mask  for  the 
protection  of  the  first  grating.  Fig.  4-23(b). 

The  substrate  is  now  ready  for  the  second  holographic  exposure.  The 
substrate  is  aligned,  exposed,  and  developed.  Fig.  s-23(c),  as  described 
earlier.  At  this  point,  the  substrate  is  again  rotarv  ion  milled  in  order  to 
replicate  the  photoresist  grating  into  the  surface.  Fig.  4-23(d).  Whatever 
residual  photoresist  remains  is  then  removed  in  an  o.,  plasma  (power  density  of 
~8  W/cm^),  Fig.  4-23(e).  What  remains  are  two  precisely  placed  diffraction 

gratings  with  different  periodicities  and  orientations  on  a  common  substrate. 

This  technique  can  be  repeated  any  number  of  times,  in  the  integration 
of  an  arbitrary  number  of  thin  film  optical  devices,  each  requiring 
diffraction  gratings  of  differing  periodicity,  orientation,  depth,  and 
geomet  rv . 

This  multiple-placement  technique  has  proven  to  be  useful  in  the  inte¬ 
gration  of  various  thin-fllm  optical  devices  as  a  result  of  the  resolution  and 
placement  capabilities  of  high- resolut ion  photolithography  (placement 
acc  racies  of  ~3  to  5  pm  and  edge  resolutions  of  ~1  pm).  In  addition,  grating 

orientat  ions  can  be  controlled  to  an  accuracy  of  0.16  deg  of  arc  and  grating 

-4 

periodicities  to  accuracies  of  three  parts  in  10  pm. 

4.3  LOCKING  OF  I  A3ER  STRIPE  ARRAYS 

The  concept  of  the  high-brightness  laser  included  obtaining  power 
levels  of  0.1  to  1.0  W  by  locking  several  DBR  lasers  together.  In  parallel 
with  DBR  laser  development,  some  initial  studies  were  made  on  locking  of 
Fabrv-Perot  stripe  geometry  lasers.  A  theoretical  study  of  coupled  lasers 
appears  in  Appendix  E.  The  geometry  of  the  first  laser  array  studied  is  shown 
in  Fig.  4-24. 

The  lasers  were  pulsed,  and  device  behavior  was  mapped  out  from 
near  field  profiles  obtained  as  a  function  of  frequency  for  different  input 


4-35 


V 


I 


1 


(a) 


\A.  '‘w  V  u/~V  \\  \r\.  "V'V 


(c) 


AWvV w 


Cl 


(d) 


Fig.  4-23.  Multiple  Grating  Fabrication 
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Optical  coupling  was  indicated  by  noting  the  development  of  a 
new  laser  frequency  in  one  stripe  due  to  the  presence  of  adjacent  stripes 
(Fig.  4-25).  At  0.38  V,  only  stripe  A  was  lasing,  at  frequency  V2,  and  the 
spectrum  of  stripe  B  showed  only  a  weak  component  at  frequency  V2  because  of 
optical  coupling  from  laser  A.  The  threshold  for  stripe  B  was  0.39  V,  at 
which  point  frequency  began  to  lase,  taking  over  and  becoming  dominant  as 
the  current  was  increased.  This  new  laser  frequency  was  optically  coupled 
to  stripe  A  and  caused  laser  action  at  frequency  in  stripe  A. 

This  interpretation  was  verified  by  observing  near-field  profiles  at 
the  two  frequencies  (Fig-  4-26).  Frequency  V2  occurred  strongly  only  in 
stripe  A,  while  at  frequency  v7 ,  all  three  stripes  were  lasing  at  high  current 
levels.  The  double  peak  at  frequency  in  stripe  A  was  because  of  depletion 
resulting  from  action  at  frequency  v,.  Suppression  of  frequency  \>2,  such  as 
with  use  of  DBR,  would  result  in  three  coupled  stripes  lasing  only  at 
frequency  Vj. 

A  model  for  optical  coupling  between  the  laser  stripes  considers  direc¬ 
tionally  coupled  single  mode  waveguides.  In  the  strong  coupling  limit,  the 
directional  coupling  coefficient  becomes 
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o 


exp(-Aew)ko  s 


1/2 
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Basov,  ^  who  showed  that  lasers  of  length  L  can  be  locked  as  long  as  they 
differ  in  wavelength  by  less  than  6X 

where 


6X 

X  <  r/2nL 


(4-10) 


and  r  is  the  exchange  of  energy  between  adjacent  lasers.  Since  the  experi¬ 
mentally  observed  optical  coupling  was  on  the  order  of  5%,  the  estimated 
locking  range  is  0.04  A.  If  the  parallel  DH  lasers  are  not  close  enough  in 
frequency,  they  may  be  tuned  with  separate  excitation  of  the  individual  laser 
stripe.  This  varies  the  local  temperature  or  refractive  index  the  small 
amount  necessary  to  lock  the  lasers. 

Additional  measurements  were  completed  to  demonstrate  locking  by  using 
a  high-resolution  spectrometer  and  comparing  individual  laser  mode  frequencies 
of  two  adjacent  laser  stripes.  The  device  geometry  was  the  same  as  that  shown 
in  Fig.  4-24.  In  order  to  prove  locking,  it  is  necessary  to  demonstrate  that 
both  laser  stripes  oscillate  at  the  same  frequency.  This  is  discussed  in  more 
detail  in  Appendix  E.  For  these  reasons,  we  carefully  investigated  the  laser 
frequencies  to  measure  locking  and  studied  devices  L10B  and  L10C. 

The  spectrum  of  laser  stripe  arrays  made  from  L10B  is  shown  in  Fig.  4- 
28.  The  top  trace  is  that  of  stripe  A  and  the  bottom  trace  is  that  of  stripe 
B.  Since  these  stripes  were  pumped  in  parallel,  it  was  not  possible  to  look 
at  each  spectrum  without  the  effects  of  the  adjacent  stripe.  However,  it  is 
clear  that  stripe  A  has  a  spectrum  that  peaks  about  0.8625  pm,  whereas 
strip  B  has  a  spectrum  that  peaks  about  0.8607  pm.  In  addition  to  the  peak 


N.E.  Basov,  E.M.  Belenov,  and  V.S,  Letokhov,  "Diffraction  Synchonization  of 
Lasers,"  Sov.  Phys.-Tech.  Phys.,  845  (1965). 
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Fig.  4-28.  Spectra  for  Unlocked  Coupled  Lasers.  (a)  Output  from 
stripe  A.  (b)  Output  from  stripe  B. 
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wavelength.  It  can  be  seen  that  each  stripe  contains  some  light  at  the 
wavelength  of  the  adjacent  stripe,  causing  the  spectrum  of  stripe  A  to  have  a 
subsidiary  maximum  at  0.8610  pm  and  stripe  B  to  have  a  subsidiary  maximum  c 
0.8623  pm.  These  subsidiary  maxima  are  caused  by  optical  coupling  between  the 
adjacent  laser  stripes,  which,  in  this  case,  were  separated  by  25  Pm. 

When  the  laser  stripes  are  closer  together,  they  are  locked.  This  is 
demonstrated  by  the  spectra  in  Fig.  4-25.  These  data  were  obtained  on  sample 
L10D  whose  stripes  are  separated  by  13  Pm.  In  this  case,  both  lasers  have 
very  similar  spectra,  proof  that  they  are  locked  together.  The  spectra  of 
locked  lasers  are  particularly  simple.  The  numerical  analysis  of  this 
physical  situation  is  discussed  in  more  detail  in  Appendix  E. 

4.4  BEAM  DEFLECTION  AND  EXPANSION 

Beam  expansion  can  be  achieved  with  the  use  of  a  distributed  Bragg  de¬ 
flector  (DBD).  This  device  consists  of  a  corrugated  waveguide,  with  a  grating 
slanted  at  an  angle  0  with  respect  to  the  incident  beam.  The  geometry  of  the 
DBD  is  shown  in  Fig.  4-29.  The  incident  beam  enters  from  the  left  into  a 
confined  slab  guide  with  a  corrugated  surface.  Bragg  reflection  occurs  from 
the  grating,  reflecting  the  beam  at  an  angle  equal  to  the  incidence  angle. 
However,  the  reflections  from  the  corrugations  add  up  coherently  only  If  the 
Bragg  condition  Is  satisfied,  that  is. 


A 


_ Xm _ 

2n  sin  (0/2) 


(4-11) 


In  addition  to  beam  expansion,  applications  of  the  DBD  include  beam 
splitters  (power  division),  deflectors  (as,  for  example,  to  create  ring 
geometries),  polarizers  and  analyzers  (with  90  deg  deflection),  and  multi- 
plexers-demultiplexers  (with  chirped  gratings). 

The  theoretical  analysis  of  the  device  was  completed  during  the  initial 
phases  of  the  contract  and  is  presented  in  Appendix  F.  For  this  reason,  only 
the  results  of  that  theoretical  analysis  are  presented  in  the  section. 
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It  was  found  that  the  incident-light-beam  amplitude  is  reduced  while 
traveling  through  the  DBD  region  by  an  exponential  factor  with  distance,  that 
is,  exp(-Yz),  because  of  deflection.  The  results  of  the  calculation  are: 


and 


where 


%  "  V  COS0B 


(4-12) 
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log  I—  -J  (cos  6  -  1)  -  cos  0. 
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ln2  -  nLJ  1 


(4-13) 
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(4-14) 


The  quantities  are  as  defined  in  Fig.  4-29,  with  Xq  the  wavelength  of  incident 

light,  l  the  Bragg  scattering  order,  and  te  the  effective  guide  thickness, 

(t  =  t  +  1/5,  where  5  is  the  exponential  decay  constant  of  the  mode  in  the 
e 

substrate). 

Consider  the  deflection  of  a  TM  polarized  beam  at  90  deg  which  is  the 
configuration  of  greatest  interest  for  beam  expander  applications: 

(4-15) 


We  have  made  use  of  the  fact  that  nj  ■  air  “l. 

It  is  convenient  to  relate  these  results  to  those  of  retroreflection  in 
a  corrugated  waveguide.  By  means  of  the  coupled  mode  theory,  it  is  predicted 
that  a  corrugated  waveguide  will  couple  forward  and  backward  waveguide  modes 
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with  a  coupling  coefficient  K.  Light  originally  travelling  forward  transfers 
to  the  backward  wave  exponentially  with  length,  with  a  coefficient  given  by 
K.  Thus,  the  effectiveness  of  a  grating  can  be  described  by  describing  its 
coupling  coefficient.  For  square  corrugations  such  as  those  used  to  derive 
Eq.  (4-13),  the  coupling  coefficient  is 


2  2 

2v  An  (d_\ 

“  3  (21  -  1)  Xn  lt  > 
g  e 


(4-16) 


A  related  parameter  is  the  characteristic  distance  at  which  the  power  in  the 
initial  mode  decreases  by  a  factor  of  1/e.  This  characteristic  length  is 
given  by  lg  =  (2K)-^.  We  then  write  the  DBD  exponential  coefficient  as 
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Since  the  deflected  beam  has  an  amplitude  that  decreases  with  distance 
along  the  DBD  according  to  an  exponential  with  a  coefficient  Y,  the  intensity 
width  of  the  expanded  beam  will  be  approximately  1/2  Y. 

Defining  the  expansion  factor  as  the  width  of  the  deflected  beam 
divided  by  the  incident  beam  width,  we  write  the  expansion  factor  as 
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where  we  have  defined  the  grating  by  its  characteristic  length 


It  is  easy  to  obtain  large  expansion  factors,  since  the  characteristic 
length  of  gratings  can  be  typically  larger  than  W.  The  research  program  was 
designed  first  to  demonstrate  beam  expansion  in  high  optical  quality  vapor- 
phase  epitaxy  (VPE)  GaAs  waveguide  material  and  then  to  apply  the  concept  to 
liquid-phase  epitaxy  (LPE)  GaAs  double  heterostructures,  so  that  the  laser  and 
beam  expander  could  be  combined  into  the  same  substrate. 

Demonstration  of  beam  expansion  in  VPE  GaAs  epilayers  was  successful. 

A  perspective  view  of  the  Bragg-effect  polarizer-analyzer  that  we  fabricated 

is  shown  in  Fig.  4-30.  The  asymmetric  slab  waveguide  structure  consisted 

K.  16  -3 

of  a  1-pm  thick,  n-type,  ~10  cm  GaAs  layer  epitaxially  grown  over  a 

18  -3 

heavily  doped  n-type,  ~5  x  10  cm  GaAs  substrate;  it  supported  only  the 
lowest  order  one-dimensional,  transverse  electric  (TE)  and  transverse  magnetic 
(TM)  polarized  modes.  Ninety-degree  Bragg  scattering  of  the  incident 

1.15-pm  radiation  was  accomplished  by  means  of  a  second-order  periodic  surface 
corrugation  (A  =  0.478  pm)  that  was  ion  milled  into  the  waveguide  surface  and 
overcoated  with  a  0.1-pm-thick  layer  of  aluminum  in  order  to  increase  its 

deflection  efficiency.  As  revealed  by  SEM  examination  of  a  cross  section 

through  the  device,  the  corrugations  bore  a  nearly  triangular  tooth  profile. 
The  total  corrugation  length  L  and  input  beam  width  W  were  3.2  mm  and  50  pm, 
respectively. 

A  semilogarithmic  plot  of  the  near-field  deflected-beam  intensity 
profile  of  one  of  our  90-deg  Bragg  deflectors  is  shown  in  Fig.  4-31.  The  data 
indicate  that  the  deflected  beam  decreases  by  a  factor  of  1/e  in  1.3  mm 
length.  With  an  input  beam  of  50  pm,  this  corresponds  to  an  expansion  factor 
of  38.  From  these  data,  it  can  be  inferred  that  the  effective  coupling  length 
of  the  grating  was  Lg  =  385  pm,  with  a  coupling  coefficient  of  K  =  0.0013. 
The  reason  for  the  weak  coupling  coefficient  is  the  use  of  VPE  GaAs.  The 
index  change  caused  by  the  reduction  of  the  number  of  free  carriers  in  the 

epi  layer  is  very  small.  As  a  result,  the  guided  mode  extends  into  the 
substrate,  and  the  effective  guide  width  is  larger  than  the  1-pm  physical 
thickness  of  the  epilayer.  For  the  physical  parameters  used  in  the 
experiment,  the  effective  guide  thickness  was  several  microns,  causing  the 
coupling  coefficient  to  be  small. 
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Fig.  A— 3 1  -  Semilogarithmic  Plot  of  Near-Field  Bragg  Scattered  Beam  Intensity 
Profile  of  90-deg  Deflector 


Ji. 


I 


V 


4-50 


Because  of  the  weakness  of  guiding  in  the  VPE  material,  the  effective 
waveguide  thickness  was  large  and  the  grating  coupling  coefficient  was 
small.  We  demonstrated  that  even  a  grating  with  a  low-coupling  coefficient 
can  be  used  successfully  as  a  polarizer-analyzer.  A  study  of  this  effect  was 
undertaken;  the  results  of  this  study  have  been  submitted  for  publication®  and 
are  summarized  in  this  report. 

Operation  of  the  DBD  as  a  polarizer-analyzer  can  be  seen  by  observing 
that  0  f°r  a  deflection  angle  of  90  deg.  Polarizer-analyzer  performance 
parameters  of  interest  include  the  extinction  ratio  p,  which  we  define  as 
being  equal  to  the  ratio  of  TM  to  TE  mode  power  deflected  by  the  device,  or 

yTM  6TE  1  "  exp(_2  6tmL) 

P  -  yTE  d™  1  -  exp(-2  ^L)  (4"19> 


and  the  TM  deflection  efficiency  n»  which  we  define  as  being  equal  to  the 
ratio  of  deflected  to  incident  TM  mode  power,  or 


n  =  [i  ~  exp(-2  6tmL)  ]  (4-20) 

In  Eqs.  (4-19)  and  (4-20),  <$TE  and  6  are  the  total  (including  background)  TE 
and  TM  mode  attenuation  coefficients,  respectively,  and  L  is  the  device 
length. 

Experiments  were  performed  on  the  polarization  dependence  of  beam  de¬ 
flection  using  the  VPE  GaAs  waveguide  beam  expander  shown  in  Fig.  4-30. 
By  measuring  the  near-field  intensity  profile  of  the  beam  which  was  Bragg 


g 

H.M.  Stoll  and  W.E.  Soady,  "Bragg  Effect  Polarizer/Analyzers  and  Lensless 
Waveguide  Beam  Expanders,"  presented  at  Integrated  Optics  Conference,  Lake 
Tahoe,  1980. 
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deflected  by  the  polarizer-analyzer,  6p.  was  h^ermined  to  be  4.19  cm  x. 
Knowledge  of  6™,  together  ”ith  measurements  of  the  undeflected  throughput 


power 


P  (0)  =  P.  [cos  0  exp(-2  6  L) 
out  in  TM 


+  sin  U  exp(-2  6TgL)] 


(4-21) 


0  =  0  and  tt/2  were  then  used  to  yield  6^,^  =  3.11  cm  y  =  0.154  cm  ^ ,  and 

n  =  0.034.  Direct  measurement  yielded  an  extinction  ratio  of  343  from  which 

-4  -1 

it  was  inferred  that  =  3.6  x  10  cm  .  In  Fig-  4-32,  the  Pout(0) 
and  ^-^5(0)  (the  total  TM  power  deflected)  are  plotted,  and  the  experimental 
data  summarized. 

The  passive  beam  expander  in  VPE  GaAs  demonstrated  that  beam  expansion 
was  indeed  possible  with  the  DBD.  However,  it  was  also  demonstrated  that  the 
efficiency  of  deflection  goes  to  zero  for  the  TE  mode  at  a  deflection  angle  of 
90  deg. 

In  Figure  4-33,  a  DBD  at  45  deg,  which  was  fabricated  in  VPE  GaAs,  is 
shown.  The  beam  expansion  is  shown  in  Fig.  4-33  (bottom)  as  a  streak  to  the 
left  of  the  beam.  The  difficulty  in  observing  this  device  is  the  necessity  of 
polishing  the  end  face  to  obtain  the  output  at  45-deg.  The  advantage  is  that 
the  grating  can  be  used  in  fundamental  order.  The  deflection  of  the  DBD  in 
the  VPE  waveguide  was  very  weak,  because  the  waveguide  was  so  thick.  Experi¬ 
ments  on  45-deg  deflection  in  thinner  LPE  GaAs  waveguides  are  planned  in  the 
future. 


To  determine  the  required  grating  spacing  for  a  DBD,  it  is  necessary  to 

measure  0.,  the  propagation  constant  in  the  guide.  For  the  VPE  grown  GaAs 

1  -1 
waveguides,  the  measured  value  of  B  was  18.887  ±  0.002  pm  ,  determined  for 

the  output  angle  of  43.23  deg  for  a  grating  spacing  of  0.4217  pm.  The 

experimental  setup  used  to  make  this  measurement  is  shown  in  Fig.  4-34. 


POWER  arb  units 


INPUT  POLARIZATION  ANGLE  deg 


Pig.  4-32.  Plots  of  Throughput  Power  P(  (0)  and  Bragg-Scattered  Power 
P-p,^  (0)n  as  Functions  of  Input-Bean  Polarization  Angle  0. 
Also  shown  art*  experimental  data  A  and  0. 
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Fig.  4-34.  Output  Grating  Coupler  Used  to  Measure  Waveguide  Propagation 
Constant  B.  (a)  Top  view.  (b)  Side  view. 


Knowledge  of  the  output  coupler  period  to  one  part 
mination  of  the  guided-wave  propagation  constant 
according  to: 


„  2 it  2iT  cos  6 

e  “  J  +  - x — 

oc 

where  A  is  the  output  coupler  period,  A  is  the 

wavelength,  and  0  is  the  output  extraction  angle. 

requisite  DBD  grating  period,  A  ,  for  second  order 

DBD 

given  by: 

A  =  — 

DBD  28 


in  10^  permitted  deter- 
8,  to  similar  accuracy 

(4-22) 

free-space  measurement 
Once  8  is  known,  the 
90  deg  deflection  can  be 
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5.  CONCLUSIONS  AND  SUMMARY 


Although  the  goal  of  building  a  fractional  watt  laser  array  with  low-beam 
civergence  was  not  met,  substantial  progress  was  made.  Distributed  Bragg 
reflector  (DBR)  lasers  were  fabricated  and  tested.  Distributed  Bragg 

deflectors  (DBD)  were  fabricated  and  proved  to  be  useful  for  beam  expanders 
and  polarizer-analyzer  elements  in  optical  circuits.  Output  couplers  and  DBDs 
were  integrated  on  the  same  chip,  and  experiments  were  performed  on  coupled- 
stripe  geometry  lasers  to  study  the  basic  physics  of  injection  locking. 

Device  integration,  in  the  form  of  a  high-brightness  laser,  requires  the 
ability  to  routinely  fabricate  DBR  lasers,  DBD  gratings,  and  output  couplers, 
and  to  be  able  to  place  such  devices  on  a  single  crystal  with  close 
tolerances.  Current  laboratory  efforts  are  directed  toward  increasing  the 
vield  of  our  devices  by  refining  our  fabrication  processes. 

In  addition,  device  integration  requires  suitable  material  free  of 
irregularities  with  good  optical  and  electrical  qualities  over  large  areas. 
Toward  this  end,  in  addition  to  refining  our  liquid  phase  epitaxy  (LPE) 
capability,  we  are  constructing  a  metal  organic  chemical  vapor  deposition 
(MOCVD)  system  with  which  to  grow  material. 

We  are  hopeful  that  our  current  efforts  will  result  in  a  high-brightness 
GaAlAs  laser. 
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LPE  CRYSTAL  GROWTH  PROCEDURE 


A.  1 
A.  1.  1 


A.  1.2 


MELT  PREPARATION 

Source  Material 

1.  Sn.  5-9s  pure  1/8-in. -diameter  shot. 

2.  Ge.  High-purity  ingot,  small  broken  pieces. 

3.  A1  5-9s  pure  wire,  0.020-in.  diameter. 

1C  O 

4.  GaAs.  (100)  orientation,  7.5  x  10  /cm  concentration, 

EPD  «  1000-2000/cm^;  0.018-in.  thick.  Cleaved  into 

small  pieces. 

5.  Ga.  7-9s  pure.  Poured  as  a  liquid  from  a  plastic  squeeze 
hot  tie. 

Weighing 

1.  Cut  or  cleave  dopants  into  pieces  weighed  on  analytical 
balance  to  ±0.0005  g.  A1  is  weighed  to  ±0.0001  g. 

2.  Ga  is  weighed  in  an  HF-cleaned  glass  beaker  as  a  liquid.  With 
the  use  of  an  evedropper,  it  can  be  weighed  to  ±0.0010  g. 

Cleaning 

1.  Dopants  Sn  and  Ge  and  GaAs  are  placed  in  a  clean  beaker  with 
tot  isopropyl  alcohol  for  approximately  1/2  hr.  Dopants  are 
i hen  Ar-dried  and  placed  in  a  covered  petri  dish. 

2.  A1  is  not  put  into  isopropyl  alcohol  because  it  can  react  with 
any  trace  water  to  form  a  hydroxide. 

3.  Ga,  weighed  as  a  liquid  in  a  glass  beaker,  has  H2O  plus  5  tc 
10  drops  of  HC1  added.  It  is  heated  for  about  10  min.  Then, 
the  beaker  is  placed  in  an  ultrasonic  cleaner  for  cleaning  in 
xylene  and  rinsing  in  isopropyl  alcohol.  After  an  Ar-dry,  the 
thickness  is  measured. 
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4.  Both  the  substrate  and  lead  crystals  are  then  cleaned  in  hot 
isopropyl  alcohol  held  in  a  vertical  position  by  a  glass 
holder  submerged  in  isopropyl  alcohol.  They  soak  for  about 
1  hr  while  the  melt  material  is  being  weighed.  They  are 
rubbed  with  Q-tlps  wet  with  alcohol  and  rinsed  in  the  hot 
isopropyl  alcohol  several  times  to  rid  the  surface  of  any 
dirt.  After  argon  drying,  the  samples  are  then  inspected  by 
microscope.  Sometimes  they  are  cleaned  in  isopropyl  alcohol 
in  the  ultrasonic  cleaner,  but  this  is  avoided  if  possible 
because  the  ultrasonic  sometimes  causes  breakage.  The 
crystals  are  Ar-dried  and  kept  covered  until  loaded  into  the 
boat . 


A. 2  BOAT 

1.  A  six-well  graphite  boat  has  been  used  with  a  sliding  graphite 
plug  top.  This  plug  top  can  be  pushed  forward  to  force  the 
Al,  wnich  is  resting  on  the  top  edge  of  the  melt  hole,  to  fall 
into  the  molten  melts  at  elevated  temperatures  (600°C). 

2.  The  boat  is  shown  in  Figs.  A-l  through  A-4,  and  drawings  are 
included  at  the  end  of  this  appendix. 

3.  The  boat  is  kept  under  in  the  quartz  furnace  tube  until  it 

is  loaded,  to  help  keep  excess  O2  from  absorbing  into  the 
porous  graphite. 

4.  The  Sn  and  Ge  are  placed  on  the  bottom  of  the  melt  hole,  then 
the  solid  Ga  is  put  on  top,  with  the  small  cleaved  pieces  of 
GaAs  source  material  placed  on  top  of  the  solid  Ga.  The  Al  is 
put  on  the  top  edge  of  the  graphite,  so  it  is  not  in  contact 
with  melt,  which  has  oxides  that  react  with  Al.  To  rid  the 
boat  of  any  water  vapor,  the  melts  are  taken  to  600°C  for  10 
to  15  min  to  dissociate  Ga  oxides.  At  this  temperature,  Al 
reacts  with  graphite  to  form  Al^C^,  which  is  insoluble  in  Ca 
and  also  is  a  source  of  contamination  to  the  substrate  sur¬ 
face.  Therefore,  the  soaking  time  is  kept  at  15  min  or  less 
before  the  Al  is  pushed  into  the  oxide-free  melts.  The 
furnace  temperature  is  then  taken  to  860 °C  and  held  for  about 
4  hr. 

5.  Cleaved  small  pieces  of  GaAs  are  used  to  float  on  top  of  Ga, 
because  one  large  plate  on  top  would  prevent  the  Al  from  fall¬ 
ing  into  the  melt.  On  the  other  hand,  very  small  powder-like 
material  would  allow  pieces  to  be  trapped  on  the  substrate 
surface  while  the  crystal  is  growing  and  would  cause  imperfec¬ 
tions. 

6.  During  the  loading  procedure,  cotton  gloves  are  worn  to  keep 
the  graphite  boat  and  quartz  furnace  from  being  contaminated. 


A- 2 


FURNACE 


A.  3 


1.  The  Thermco  electric  furnace,  tube  type,  Kanthal  element,  200 
to  1200°C  control,  three-zone  control  has  good  temperature 
stability.  It  has  a  flat  temperature  zone  of  approximately 
10-in.  (center).  An  induced  temperature  increase  of  l^C  from 
the  first  melt  to  the  last  melt  keeps  etch-back  from 
occurring.  Little  or  no  temperature  gradient  can  be  measured 
from  top  to  bottom  of  the  melt  area. 

2.  A  quartz  tube  approximately  2-1/2  in.  in  diameter  and  5  ft 
long  is  used  to  hold  the  boat  under  a  pure  H2  atmosphere  while 
the  melt  is  in  the  boat. 

3.  H2  is  put  through  a  Pd  purifier  before  entering  the  furnace. 
It  has  an  oxygen-free  purity  of  parts  per  billion. 

4.  The  temperature  programmer  is  a  Eurotherm  Type  125,  which  is 
solid  state,  permitting  a  very  smooth  temperature  drop  of  0.1 
to  l°C/min. 

A. 4  PROGRAM 

1.  The  furnace  is  evacuated  after  the  boat  is  placed  in  its 
center  and  back-filled  with  N.  It  is  evacuated  the  second 
time,  then  filled  with  pure  hydrogen.  The  use  of  this  vacuum 
permits  leaks  to  be  detected. 

2.  The  furnace  is  turned  on  to  400°C  with  H2  bubbling  through  at 
about  1/2  1/min.  It  is  soaked  overnight  at  400°C  to  bake  out 
any  water  vapor. 

3.  The  next  morning,  the  furnace  is  raised  to  600°C  for  1  hr. 
Then  the  A1  is  pushed  into  the  melts,  which  are  fairly  free  of 
oxides  at  this  temperature.  This  greatly  helps  in  limiting 
the  very  active  A1  from  oxidizing,  which,  in  turn,  keeps  the 
AI2C3  from  falling  on  the  growing  surface. 

4.  The  furnace  temperature  ■  is  then  raised  to  800  to  850°C  and 
soaked  for  3  to  4  hr. 

5.  The  lead  crystal,  which  precedes  the  substrate  by  one  melt,  is 
pulled  under  the  first  melt  and  left  soaking  for  10  to  30  min 
to  ensure  against  melt  back. 

6.  The  programmer  is  turned  on  to  a  cooling  rate  of  about  0.0004 
mV /min,  which  is  approximately  0. 3°C/min.  A  100%  Pt-13%  Rh-Pt 
thermocouple  is  used.  The  temperature  is  lowered  2  to  3°C  to 
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supercool  the  melts,  which  may  help  smoothness  and  flatness  at 
the  substrate  surface  because  of  the  rapid  growth  of  GaAs  on 
the  surface. 

7.  The  substrate  is  then  pulled  to  the  first  melt.  The  tempera¬ 
ture  drop  is  from  0.004  to  0.002  mV/min  or  0.3  to  0.15°C/min, 
depending  on  desired  growth  rate  for  each  layer.  The  active 
laser  layer  usually  grows  xor  5  to  20  sec  at  a  rate  of 
0.3°C/min.  This  produces  a  layer  thickness  of  600  to  3000  A. 

8.  After  the  last  layer  is  grown,  the  substrate  is  pulled  to  the 
front  of  the  boat,  where  it  remains  under  the  graphite  melt 
holder  until  it  is  removed  from  the  furnace. 

9.  After  the  furnace  has  reached  room  temperature,  N2  is  used  to 
purge  it  of  H?.  The  boat  is  removed  from  the  furnace,  and  the 
substrate  and  lead  crystals  are  placed  in  hot  isopropyl 
alcohol  to  rid  their  surfaces  of  excess  Ga  melt.  A  Q-tip  is 
used  to  rub  off  Ga. 

10.  The  total  time  required  for  one  LPE  growth  and  a  preliminary 
crystal  analysis  is  3  days. 

SUMMARY 

1.  The  substrate  is  not  pre-etched  in  H2SO4  before  growth. 

2.  Trichlorethylene  is  not  used  to  clean  the  substrate. 

3.  A1  is  pushed  into  melts  at  600°C  to  reduce  oxidation. 

4.  Melts  are  supercooled  from  0.5  to  j°C  to  help  eliminate  etch- 
back  and  to  help  grow  smooth  surface  layers  without  meniscus 
lines. 

5.  Slow  growing  rate  with  smooth  cooling  curve  helps  eliminate 
dislocations. 

6.  Cleaning  Ga  in  +  HCl  in  the  ultrasonic  cleaner  rids  excess 
oxidation.  Freezing  Ga  in  1^0  +  1  to  2  drops  of  HCl  reduces 
the  oxidation  rate  of  Ga  while  loading  it  in  the  boat  in  air. 

7.  Sellows  valves  and  leak-tight  fittings  eliminate  09  leaks. 
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OPTIMALLY  COUPLED  GaAs  DBR  LASERS 

B.l  INTRODUCTION 

Optimal  designs  are  developed  for  large  optical  cavity  (LOC)  GaAs 
DBR  lasers.  Cavity  and  reflector  lengths  were  determined  that  permit  single¬ 
mode  single-ended  operation  over  a  continuous-wave  optical  output  power  range 
of  5  to  50  mW.  With  an  internal  quantum  efficiency  of  90%  assumed,  external 
quantum  efficiencies  for  the  optimized  structures  are  found  to  range  from 

14.1%  (for  5-raM  output)  to  57.8%  (for  50-mW  output).  Threshold  current  densi- 
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ties  are  found  to  remain  less  than  1000  amp/cm  over  the  entire  optical  output 
range  considered. 

Two  types  of  periodic  feedback  lasers  have  been  proposed  for  use  as 

sources  in  monolithic  semiconductor  integrated  optical  circuits:  the  distrib- 
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uted  feedback  (DFB)  laser  and  the  DBR  laser.  While  both  conf igurations 
offer  obvious  advantages  over  the  conventional  Fabry-Perot  semiconductor  laser 
(e.g.,  freedom  of  orientation  within  the  integrated  optical  circuit  and  lack 
of  facet  pitting  at  high  continuous-wave  output  levels),  the  DBR  laser 
appears  to  be  superior  to  the  DFB  laser  for  a  number  of  operational  and  fabri- 
cative  reasons.  First,  because  of  the  spectral  selectivity  of  its  mirrors, 
the  DBR  laser  may  be  designed  to  lase  in  a  single  longitudinal  mode  instead  of 
in  the  longitudinal  mode  pair  characteristic  of  the  index-periodic  DFB  laser. ^ 
Second,  because  the  laser  reflectors  can  be  located  away  from  the  active 
(pumped)  region  of  the  device,  nonradiative  recombination  centers  introduced 
during  fabrication  of  the  DBR  laser  reflector  gratings  can  be  disregarded 
(the  feedback  grating  of  the  DFB  laser  must  be  contiguous  to  or  very  near  the 
active,  p-n  junction).  Third,  because  the  DBR  laser  reflectors  are  easily 
accessible,  they  can  be  individually  tailored  in  order  to  achieve  single-ended 
single-longitudinal-mode  operation.  Fourth,  unlike  the  DFB  laser,  the  semi¬ 
conductor  DBR  laser  can  be  fabricated  by  means  of  a  single  planar  epitaxial 
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growth  sequence. 
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Of  the  optically^  and  electrically^-^ ’ pumped  DBR  lasers 

fabricated  to  date,  we  choose  to  analyze  the  large  optical  cavity  (LOC)  design 
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first  proposed  and  demonstrated  by  Wang  et  al.  •  This  configuration  is  com¬ 
pletely  amenable  to  batch  planar  processing  and  may,  as  will  be  shown,  be 
designed  to  operate  at  pump  current  thresholds  of  less  than  1000  amp/cm  . 
Within  certain  parameter  constraints  (e.g.,  minimum  practical  or  "comfortable" 
active  and  waveguide  epilayer  thicknesses),  an  LOC  DBR  laser  design  is  devel¬ 
oped  that  simultaneously  maximizes  the  external  quantum  efficiency  of  the 
device  and  provides  for  single  transverse  and  longitudinal  mode  operation  over 
a  continuous-wave  optical  output  power  range  of  5  to  50  mW. 

B. 2  BASIC  DESIGN  CONSIDERATIONS 

A  cross  section  through  the  DBR  laser  cavity  to  be  analyzed  is 
shown  in  Fig.  B-l(a).  The  central  portion  of  the  device,  the  active  region, 
is  electrically  pumped  and  experiences  gain.  The  two  end  regions  of  the 
device,  the  passive  waveguiding  regions,  accommodate  the  distributed  Bragg 
reflectors  (DBRs )  and  are  designed  to  be  transparent  to  the  radiation  emitted 
by  the  active  region.  In  order  to  achieve  efficient,  single-longitudinal-mode 
operation,  one  DBR  must  bear  low  reflectivity  for  optimal  output  coupling  and 
the  other  DBR  must  bear  narrow  bandwidth  and  high  reflectivity  for  single¬ 
mode,  single-ended  emission.  Hybrid  Fabry-Perot/DBR  laser  configurations 
wherein  one  reflector  is  a  crystal  facet  are  not  considered,  either  because  of 
the  difficulty  involved  in  optimally  adjusting  the  reflectivity  of  a  cleaved 
facet  (assuming  the  DBR  served  as  the  high-reflectivity  narrow  bandwidth 
reflector)  or  because  of  the  lack  of  spectral  selectivity  associated  with 
low-reflectivity  DBRs  (assuming  the  cleaved  facet  were  coated  for  high 
reflectivity). 

The  DBR  laser  configuration  shown  in  Fig.  B-l(a)  may  be  analyzed  by 
considering  the  formally  equivalent  structure  shown  in  Fig.  B-l(b).^  Here, 
the  DBRs  are  replaced  by  discrete  narrow-bandwidth  mirrors  bearing  frequency- 
dependent  complex  reflectivities 
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(B-l) 


R1  2  ^  =  ri  2  ^  exp  l1*!  2 

With  Eq.  (B-l)  and  accounting  for  active-to-passive  region  coupling  losses  by 
inclusion  of  the  coupling  factor  f(<l),  the  DBR  laser  mode  spectra  and  thresh¬ 
old  gains  are  found  in  the  usual  manner  by  setting  the  net  round-trip  cavity 
gain  equal  to  unity  or 

f4  r1r2  exp  [1(4^  +  <p 2  -  2BLa)]  exp  [2  (gth  -  aa>  La  ]  =  1  (B-2) 

where  gt^  and  are  the  distributed  amplitude  gain  and  loss  coefficients, 
respectively,  of  the  active  region,  is  the  length  of  the  active  region, 
and  B  is  the  propagation  constant  of  an  oscillation  mode.  The  fourth  power 
exponent  of  f  in  Eq.  (B-2)  is  justified  by  noting  that  each  encounter  with  a 
DBR  by  a  wave  circulating  within  the  laser  cavity  results  in  a  net  coupling 
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efficiency  of  f  :  one  factor  of  f  for  active-to-passive  traversal  and  one 

factor  of  f,  by  reciprocity,  for  passive-to-active  reentry.  The  meaning  of 
Leff,  the  effective  cavity  length  in  Fig.  B-l (a),  will  be  explained  in  Sec¬ 
tion  B.5.  From  Eq.  (B-2),  it  is  evident  that  basic  DBR  laser  design  consid¬ 
erations  involve  both  adjustment  of  the  active-passive  confinement  structure 
in  order  to  maximize  f  and  minimize  gtj1  and  adjustment  of  Rj  2  to  achieve 
optimally  coupled  single-mode  emission. 

B.  3  DBR 

The  DBR  consists  of  a  periodically  perturbed  section  of  waveguide 

designed  to  effect,  as  a  result  of  backward  Bragg  scattering,  retroref lection 
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of  an  incident  wave.  In  our  analysis,  we  assumed  the  geometry  of  Fig.  B- 
1(a),  wherein  the  DBR  consists  of  a  periodically  deformed  (or  corrugated) 
section  of  waveguiding  medium.  The  amplitude  transmittivity  T  and  amplitude 

reflectivity  R  of  such  a  structure  may  be  found  by  means  of  the  coupled  mode 
1  3 

theory.  Accordingly, 

_ Y  exp  [~i  (B  ~  AB)  L] _ 

(a  +  iAB)  sinh  (yL)  +  ycosh  (yL) 

P 

B-A 


T 
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and 


R  = 


-  i  icsinh  (yL) 


(a  +  iAB)  sinh  (yL)  +  Ycosh  (yL) 
P 
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where  is  the  distributed  loss  coefficient  in  the  waveguiding  (passive) 

region,  L  is  the  DBR  length,  k  is  the  corrugation  coupling  coefficient,  12,14 
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YZ  =  <  +  +  iAB)  , 


and  AB  is  the  measure  of  deviation  from  exact  satisfac¬ 


tion  of  the  Bragg  condition  given  by 


0)  -  to 
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where  v  is  the  phase  velocity  of  light  within  the  DBR,  m  =  1,2,3...  is  the 
Bragg  scattering  order,  and  A  is  the  period  of  the  surface  corrugation. 


The  magnitude  r  and  the  phase  tp  of  R  are  plotted  as  functions 
of  ABL  in  Figs.  B-2  and  B-3  for  several  combinations  of  the  attenuation  and 
coupling  parameters  a^L  and  kL.  The  principal  effects  of  attenuation  within 
the  DBR  are  a  reduction  in  the  line  center  (ABL  =  0)  reflectivity  magnitude 
and  less  distinct  reflectivity  sidelobes.  In  addition,  note  that  in  all 
cases,  the  reflectivity  phase  varies  almost  linearly  with  ABL  over  the  central 
reflectivity  lobe.  For  small  ABL,  the  latter  observation  permits  us  to 
approximate  <p  by 
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a  +  («  +  ap  ) 
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for  IyL|  >  2  or  by 


for  I  yL  I  S.  0.5. 

Design  optimization  (Section  B.5)  requires  knowledge  of  the 
DBR  amplitude  reflectivity  bandwidth  B  =  FWHM  of  r  (ABL),  line  center  power 
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Fig.  B-2.  DBR.  (a)  Amplitude  reflectivity  r  as 
function  of  deviation  (A0L)  from  Bragg 
scattering  condition.  (b)  Reflectivity 
phase  <J>  plotted  as  function  of  deviation 
from  Bragg  condition. 


.  DBR.  (a)  Amplitude  reflectivity  r  plotted 
as  function  of  deviation  (A&L)  from  the 
Bragg  scattering  condition.  (b)  Reflec¬ 
tivity  phase  <f>  plotted  as  function  of 
deviation  from  the  Bragg  condition. 


reflectivity  =  R^(0),  and  line  center  power  transmittivity  S  =  T^(0). 

Line  center  values  for  .•/}  and  .7  are  chosen  for  convenience;  justification  for 

their  use  is  made  in  Section  B.5.  Figure  B-4(a)  is  a  parametric  plot  of  B  as 

a  function  of  both  a^L  and  kL.  Figure  B-4(b)  is  a  composite  plot  of  both 

■*’  and  .7  is  functions  of  a  L  and  kL.  For  arbitrary  values  of  a  L  and  kL.  the 

P  P  ’ 

latter  plot  can  be  used  to  calculate  the  total  distributed  loss  within  a  given 

DBR  =  1  - -It  -  .7. 

B.4  ACTIVE-PASSIVF,  COUPLING  FACTOR 

A  cross  section  through  the  interface  that  separates  the  active 
region  of  a  DBR  laser  from  one  of  its  DBRs  is  shown  in  Fig.  B-5.  The  active 
region  is  comprised  of  five  epitaxially  grown  GaAlAs  layers,  which  together 
form  a  localized  gain  region  (LGR)  or  large  optical  cavity  (LOC)  laser  struc¬ 
ture.  * ^  The  transverse  field  distribution  of  the  active  region  mode  is 
assumed  to  be  defined  by  the  upper  cavity,  active,  waveguide,  and  isolation 
layers  and  to  be  concentrated  primarily  within  the  active  and  waveguiding  epi- 
lavers.  The  passive  region,  within  which  the  distributed  Bragg  reflectors 
reside,  is  formed  by  removing  the  cap,  upper  cavity,  and  (highly  absorptive) 
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active  layers  from  the  active  region.  The  transverse  field  distribution 
within  the  passive  region  is  determined  by  the  superstrate  (assumed  to  be  air) 
and  by  the  waveguide  and  isolation  epilayers. 

As  accounted  for  by  the  f^  factor  in  Eq.  (B-l),  different  active 
and  passive  transverse  field  distributions  result  in  incomplete  optical 
coupling  between  active  and  passive  regions  of  the  DBR  laser.  The  active 
passive  coupling  factor  f  can  be  found  in  a  straightforward  way  by  considering 
the  field  radiated  into  the  passive  region  by  the  lowest  order  active  region 
mode  incident  at  the  active  passive  interface  (it  is  assumed  that  the  active 
and  passive  regions  each  support  only  one  confined  mode).  The  lowest  order 
active  region  mode  e^°^(x)  is  first  expressed  as  an  expansion  over  the  ortho- 
normal  passive  region  mode  sot  as  follows: 
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where  c^°\x)  is  the  lowest-order  passive  region  confined  mode.  e^V^(x)  repre- 
P  P 

sents  the  continuum  of  passive  region  radiation  modes,  and  f  and  are  expan¬ 
sion  coefficients.  With  the  use  of  the  orthonormality  of  the  passive  region 
eigenmode  set,  Eq.  (B-8)  is  then  manipulated  to  yield  for  the  coupling  factor 


f  =  f  (x)  (x)  dx 

J- »  a  P 


(B-9) 


In  deriving  Eq.  (B-9),  it  has  been  assumed  that  all  modes  are  transverse  elec¬ 
tric  (TE)  polarized  (electric  field  vector  aligned  parallel  to  the  y- 
direction);  extension  to  the  transverse  magnetic  (TM)  case  is  obvious  and  will 
not  be  considered  here. 

The  (TE)  field  distributions  in  the  active  and  passive  regions  are 
found  by  solving  the  scalar  wave  equation  and  requiring  continuity  of  the 
tangential  electric  and  magnetic  fields  at  the  various  epilayer  interfaces 
(Fig.  B-5) .  Taking  the  origin  of  the  x-axis  co  be  at  the  upper  cavity  active 
layer  interface,  the  active  region  confined  mode  field  distribution  is  given 
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where  m  «  0,1,2,...  Is  the  mode  index  number  and  the  transverse  q,  h,  l,  p, 
longitudinal  6,  and  free  space  kQ  »  2tt/ Xq  propagation  constants  are  related  to 
one  another  by 
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Similarly,  for  the  three-layer  passive  region: 
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Finally,  the  normalization  constants,  C|  and  C2  are  chosen  so  that 
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B.  5  DESIGN  OPTIMIZATION 

The  DBR  laser  operational  characteristics  are  governed  by  Eq. 

(B-2).  Solution  of  Eq.  (B-2)  yields  for  the  propagation  constants  8^  and 
threshold  gains  g*1*1  of  the  DBR  laser 
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where  i  =  0,  ±1,  ±2,  ...  is  the  longitudinal  mode  index  and  AB^  the  deviation 
from  exact  satisfaction  of  the  Bragg  condition  for  the  ith  longitudinal 
mode.  Moreover,  from  Eqs.  (B-6),  (B-7),  and  (B-10),  the  longitudinal  mode 
spacing  is  given  to  good  approximation  by 
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where  the  effective  laser  cavity  length  is  given  by 
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In  Eq.  (B-18),  and  are  the  lengths  of  the  left-hand  and  right-hand, 
respectively,  DBRs,  and  power  extraction  is  assumed  to  occur  from  the  right- 
hand  low-reflectivity  DBR. 

Design  optimization  requires  single  transverse  and  longitudinal 
mode  operation.  Operation  in  a  single  transverse  mode  in  the  y-direction 
[refer  to  Fig.  (B— 5 )  ]  may  be  achieved  by  appropriate  adjustment  of  the  DBR 
laser  lateral  (stripe)  dimension;  lateral  mode  control  in  the  x-direction 
(perpendicular  to  the  junction  plane)  will  be  seen  to  result  automatically 
from  our  attempts  to  maximize  the  external  quantum  efficiency  of  the  device. 
Restriction  to  oscillation  in  a  single  longitudinal  mode  requires  that  the 
intermode  spacing,  defined  by  Eq.  (B-12),  be  equal  to  or  greater  than  the  net 
reflectivity  bandwidth  of  the  laser  cavity.  While  the  latter  quantity  is 
accurately  defined  as  being  equal  to  the  FWHM  of  r^*^,  out  of  convenience,  we 
shall  use  the  FWHM  of  r^  instead.  Our  assumption  that  the  cavity  bandwidth  is 
governed  by  the  high-reflectivity  DBR  is  justified  by  the  fact  that  in  all 
cases  of  interest  the  FWHM  of  ^  is  considerably  greater  than  the  FWHM  of  rj. 

The  DBR  design  parameters  required  for  single  longitudinal  mode 
operation  can  be  found  graphically  [Figs.  B-4(a)  and  B-4(b)].  We  assume  that 
the  Bragg  condition  has  been  exactly  satisfied  for  the  longitudinal  mode  of 
interest  (this  may  be  accomplished,  for  example,  by  thermally  varying  the 
laser  cavity  length  La  in  such  a  way  that  A0  =  0).  Certain  line-center  DBR 
reflectivities  rj  =  [->^(0)]  and  r£  =  [-^(O)]  are  then  assumed  and 

their  corresponding  parametric  curves  in  Fig.  B-4(a)  noted  as  representing  the 
loci  of  all  possible  design  solutions.  Specific  design  solutions  are  found  by 
examining  the  intersections  of  these  curves  with  straight  lines  corresponding 
to  fixed  values  of  and  *c.  The  intersection  points  yield  values  of  and 
1-2  that  can  then  be  used,  with  the  aid  of  Fig.  B-4(b)  and  Eqs.  (B-17)  and 
(B-18),  to  compute  the  longitudinal  mode  spacing  A  and  the  FWHM  of  r^. 
Once  A  has  been  found  to  be  greater  than  or  equal  to  the  FWHM  of  r^,  the 
design  problem  will  have  been  solved. 
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As  initial  assumptions,  we  take  for  the  high-reflectivity  DBK,  rj  = 
0.95  (or  'Jft  j  =  0.9),  and  for  the  low-reflectivity  DBR,  r  2  --  0.45  (or  -J/i  9  = 
0.2).  The  (lower  bound)  value  chosen  for  r^  reflects  our  wish  to  induce 

efficient  single-ended  laser  emission.  The  (upper  bound)  value  chosen  for  Ty 
anticipates  the  effects  of  a  high  active  region  gain  on  later  attempts  to 
maximize  the  external  quantum  efficiency  of  the  laser.  Parametric  curves 

for-^j  =  0.9  and  .^9  =  0.2  are  shown,  together  with  representative  DBR  band¬ 
width  curves,  in  Fig.  B-b  as  functions  of  a^L  and  kL.  Solutions  to  the  design 

problem  are  found  by  examining  the  common  intersections  of  these  two  sets  of 
curves  with  the  straight  lines  that  correspond  to  fixed  values  of  and  *c. 

The  straight  lines,  also  shown  in  Fig.  B-6,  pass  through  the  origin  and  to 

bear  slopes  of  a^/<.  Throughout  our  analysis,  we  shall  assume 

that  a  =  1  cm  '  and  that  <  is  the  same  for  both  OBRs.  ""  9  value  chosen 

P 

for  a  represents  a  lower  bound  for  GaAlAs  waveguide  att  auation  coeffi- 

1  1  h 

cients;  the  assumption  that  <  remains  the  same  for  both  DBRs  simplifies  con¬ 
siderably  the  actual  DBR  laser  fabrication  process. 

With  a  =  1  cn  *,  L  =  300  pm,  and,  as  an  initial  estimate,  <  =  20 

p  a 

cm  l,  we  find  from  Fig.  B-b  and  Eq.  ( B— 1 8 )  that  Lj  =  1.47  mm,  L2  =  240  pm,  and 
Leff  =  660  pm,  thereby  yielding  for  the  longitudinal  node  spacing,  48  cm 
and,  for  the  high  reflectivity  DBR  bandwidth,  54  cm  '.  Since  in  this  case, 
the  mode  spacing  is  less  than  the  high-reflectivity  DBR  bandwidth,  either  L_ 
or  <  must  be  adjusted  in  order  to  achieve  single  longitudinal  mode 

operation.  Increasing  <  worsens  the  situation,  whereas  a  decrease 

in  k  requires  progressively  longer  (and  more  lossy)  reflector  lengths  in  order 
to  achieve  the  same  amplitude  reflectivity  (0.95).  A  satisfactory  solution  to 
the  problem  results  from  setting  k  =  20  cm  *  and  Lg  =  275  pm,  in  which  case  Lj 
=  1.47  mm,  ho  =  240  pm,  L  =  501  pm,  and  the  iniermode  spacing  and  DBK 

bandwidth,  both  equal  54  cm  Note  that,  in  order  to  achieve  single-mode 

single-ended  laser  operation,  it  has  been  necessary  to  adjust  all  previously 
free-cavity  parameters. 

Having  determined  the  DBR  coupling  coefficient  and  cavity  length 
parameters  required  in  order  to  achieve  efficient  single  longitudinal  mode 
operation,  we  now  concentrate  on  maximizing  the  DBR  external  quantum 


efficiency  of  the  laser  by  appropriate  adjustment  of  the  epilayer  thicknesses 
and  compositions  of  the  device.  The  external  quantum  efficiency  is  maximized 
by  minimizing  the  pump  current  density  required  in  order  to  achieve  a  given 
continuous-wave  optical  output  power. 

Thompson  et  al.*^  derived  an  empirical  expression  (since  derived 
1  8 

theoretically  by  Casey  )  for  the  threshold  gain  of  lasers  bearing  active 
regions  similar  to  chose  under  consideration.  Their  expression  is 
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where  J  h  is  the  threshold  current  density  (amp/cm  ),  d  is  the  thickness  of 
the  active  region  (pm),  and  v  is  a  filling  factor  equal  to  the  fraction  of 
laser  mode  energy  that  occupies  the  active  layer, 
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The  threshold  current  density  is  found,  by  simultaneous  solution  of  Eqs. 
(B-16)  and  (B-19),  to  be 
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The  pump  current  density  Jp  required  in  order  to  achieve  a  given  optical  power 


output  PQ  is  related  to  by 


Jp  Jth  +  AE  r  n. 
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where  A  (cm  )  is  the  area  pumped,  is  the  internal  quantum  efficiency, 
is  the  active  region  bandgap  and  r  the  ratio  of  useful  to  total  loss  suffered 
per  round  trip  cavity  pass,  is 


In  deriving  Eq.  (3-22),  small  ohmic  losses  occasioned  by  electrical  contacts 
and  bulk  resistivity  have  been  ignored.  In  Eq.  (B-23),  the  first  term  on  the 
right-hand  side  accounts  for  useful  power  coupled  outside  the  active  region; 
the  second  term  accounts  for  imperfect  (i.e.,  <100%)  reflectivity  by  the  high- 
reflectivity  DBR  and  for  distributed  losses  suffered  as  a  result  of  passage 
through  both  DBRs. 

For  a  given  optical  output  power  PQ,  the  pump  current  density  Jp 

was  minimized  by  varying  the  following  parameters:  the  active  layer  thickness 

d,  the  waveguide  thickness  t,  the  refractive  index  of  the  upper  cavity  layer 

nj,  the  refractive  index  of  the  active  layer  n2»  the  refractive  index  of  the 

waveguide  layer  03,  the  refractive  index  of  the  isolation  layer  n^,  and  the 

amplitude  reflectivity  of  the  output  DBR  r2*  Throughout  the  analysis,  it  was 

assumed  that  (the  oscillation  wavelength)  =  0.9  pm,  w  (the  laser  stripe 

width)  =  10  pm  E^  =  1.38  V,  aa  =  9  cm  *^®),  ancj  ^  =  o.9.^®  The  refractive 

indices  of  various  layers  were  obtained  using  the  data  of  Casey  et  al.^  In 

addition,  the  following  constraints  were  placed  on  the  parameters  varied: 

(1)  that  the  active  layer  thickness  be  no  less  than  0.05  pm,  (2)  that  the 

active  layer  contain  0.05  mole  fraction  AlAs,  (3)  that  the  difference  in  AlAs 

mole  fraction  between  the  active  and  waveguiding  layers  be  no  less  than  0.06, 

and  (4)  that  the  AlAs  mole  fractions  of  the  upper  cavity  and  isolation  be  no 

more  than  0.5.  In  the  first  constraint  the  limitations  of  state-of-the-art 

20 

liquid-phase  epitaxy  were  recognized,  in  the  second  constraint,  the  need  for 

2  1 

strain  relief  between  the  active  and  contiguous  epilayers  was  recognized, 

the  third  constraint  ensured  reasonable  transparency  of  the  waveguide  layer  at 

0.9  pm,  and  the  fourth  constraint  indirectly  prevented  the  waveguide  layer 

2  3 

from  becoming  inconveniently  thin. 

The  epilayer  parameters  referred  to  above  were  varied  sequentially 
until  Jp  reached  a  minimum,  either  as  a  result  of  parameter  tradeoffs  or  as  a 
result  of  a  particular  parameter  reaching  one  of  the  enumerated  constraints. 
Remarkable  design  consistency  was  evidenced  over  the  entire  optical  output 
range  considered  (PQ  =  5  to  50  mW).  In  all  cases,  Jp  reached  a  local  minimum 
for  a  waveguide  thickness  of  0.225  pm  and,  simultaneously ,  was  minimized  as 
each  of  the  other  epilayer  parameters  readied  their  limiting  values. 
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In  Fig.  B-7,  both  Jp  and  are  plotted  as  functions  of  t  for  PQ 

=  50  nW.  The  limiting  (optimal)  values  of  the  remaining  epilayer  parameters, 
which  were  held  constant  as  t  was  varied,  are  also  listed.  In  addition,  the 
values  of  t2  required  in  order  to  minimize  Jp  for  a  given  waveguide  thickness 
are  shown  above  the  curve.  The  minimum  value  of  Jp,  obtained  just  prior 

to  waveguide  cutoff  (at  t  =  0.15  ins),  is  understood  to  occur  as  a  result  of 
tradeoffs  between  the  filling  factor  v  and  the  coupling  factor  f.  The  latter 
situation  is  illustrated  in  Fig.  B-8  wherein  v  is  seen  to  decrease  monotoni- 
callv  with  t,  while  f  is  seen  to  increase  monotonically  with  t.  Optimal 
values  of  v  and  f  are  noted,  from  Figs.  B-7  and  B-8,  to  be  0.113  and  0.814, 
respectively. 

The  fact  that  Jp  is  simultaneously  minimized  as  the  remaining  epi¬ 
layer  parameters  reach  their  limiting  values  is  seen  to  result  both  from  the 
strong  dependence  of  Jj,  on  the  active  layer  thickness  d,  Eq.  (B-16),  and  from 
enhancement  of  the  active-to-passive  coupling'  factor  f.  The  latter  effect 
occurs  both  because  of  improved  mode  confinement  within  the  active  and  passive 
regions  (induced  by  higher  AlAs  mole  fractions  in  the  upper  cavity  and  isola¬ 
tion  confinement  layers)  and  because  of  improved  mode  matching  between  the 
active  and  passive  regions  (induced  as  a  result  of  the  waveguide  AlAs  mole 
fraction  approaching  that  of  the  active  layer). 

In  Fig.  B-9,  plots  V  Jth>  and  r->  are  plotted  as  functions  of  PQ, 
and  the  results  of  our  design  analysis  are  summarized.  A  list  of  the  optimal 
epilayer  parameters  is  given  in  Fig.  B— 9 ;  the  coupling  coefficient  and  cavity 
length  parameters  used  were  those  derived  in  Section  B.3.  If  bulk  and  elec¬ 
trical  contact  losses  are  ignored  as  before,  the  external  quantum  efficiency 
ranges  from  14.1%  for  PQ  =  5  mW  to  57.8 Z  for  PQ  -  50  mW.  Also,  note  that 
Jtq  <1000  amp/cm  over  the  entire  optical  output  range  considered. 

B.b  CONCLUSION 

Within  certain  practical  design  constraints,  optimal  DBR  laser 
designs  have  been  developed  that  permit  single-mode  single-ended  operation 
over  a  5  to  50  mW  continuous-wave  optical  output  power  range.  It  has  been 
determined  that  only  the  reflectivity  (or  length)  of  the  output  DBR  need  be 
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Fig.  B-7 .  Threshold  and  Pump-Current  Densities  Plotted 
as  Function  of  Waveguiding  Layer  Thickness 
for  Optical  Output  Power  of  50  mW 
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Fig,  B-8.  Active  Region  Filling  Factor  v  and  Active- 
Passive  Coupling  Factor  f  Plotted  as  Func¬ 
tions  of  Waveguiding  Layer  Thickness  for 
Optical  Output  Power  of  50  nnW 
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varied  in  order  to  achieve  optimal  output  coupling  over  the  optical  power 
range  of  interest.  With  an  internal  quantum  efficiency  of  90%  assumed, 
external  quantum  efficiencies  have  been  found  to  range  from  14.1%  (5-mW  out¬ 
put)  to  57.8%  (50-mW  output),  with  threshold  current  densities  never  exceeding 
1000  amp/cm^. 
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C. 1  SAMPLE  EXAMINATION 

1.  Etch  sample  In  1:1:8  A-B  etch  (AiBtl^O  mp)  for  5  sec. 

2.  Photograph  all  four  corners  of  sample  to  determine  cap 
thickness  variations. 

C. 2  CAP  LAYER  REMOVAL 

1.  Soak  clean  sample  in  warm  (>60°C)  trichlorethylene  for  2  to  5 
min. 

2.  Rinse  sample  in  dry  isopropyl  alcohol  for  1  to  2  min. 

3.  Blow  sample  dry  with  nitrogen. 

4.  Dehydration  bake  sample  at  170°C  for  30  min. 

5.  Mount  sample  on  glass  microscope  slide  cover  with  Apiezon  wax. 

6.  Cool  and  then  whirl  clean  sample  with  a  one  to  one  mixture  of 

HMDS  and  AZ  thiner. 

7.  Apply  unthinned  filtered  AZ  1373  to  sample. 

8.  Spin  sample  at  4000  rpm  for  30  sec  (resultant  photoresist 

thickness  should  be  approximately  3  pm). 

9.  Air  dry  sample  for  15  min. 

10.  Oven  bake  sample  at  90°C  for  20  min.  (pre-exposure  bake). 

11.  Use  mask  GE-2  with  mask  alignment  instrument  and  expose 

photoresist  for  30  sec. 

12.  Develop  photoresist  for  1  min  using  equal  parts  Shipley  MF  312 
developer  and  deionized  water.  Rinse  sample  in  deionized  1^0 
for  60  sec,  and  blow  dry  with  N. 


13.  Oven  bake  sample  at  90’C  for  20  min  (postexposure  bake). 


14.  Remove  sample  from  glass  microscope  slide  cover  by  heating 
on  hot  plate.  Rinse  sample  quickly  in  xylene  at  room 
temperature. 

15.  Ion  beam  etch  sample  to  remove  cap  layer.  Mount  sample  on 
rotary  stage  with  heat  conductive  zinc  oxide  paste.  Approxi¬ 
mate  system  settings:.  0^500-mA  beam  current,  500-V 
acceleration  voltage,  7  x  10  Torr  Ar,  (resultant  milling 
rate  is  approximately  455  A/min  or  22.2  min/um).  Sample 
should  be  rotated  at  30  deg  to  the  incident  ion  beam. 

16.  Check  cap  layer  removal  by  observing  the  photoluminescence  of 
the  sample. 

C.  3  TOP  BARRIER  LAYER  REMOVAL 

1.  Chemically  etch  sample  with  A-B  etchant,  1:1:4  or  1:1:8  ratio 
of  A,  B,  and  1^0  (millipore),  for  approximately  one  min. 

2.  Check  depth  of  etching  with  photoluminescence  measurements. 

Repeat  Step  1  until  cap  layer  is  removed.  Check  integrity  of 
photoresist  after  etching. 

C.4  ACTIVE  LAYER  REMOVAL 

1.  Ion  mill  sample  to  remove  active  layer.  Use  same  parameters 
as  in  Section  C.2,  Step  15.  Or,  perform  Step  2. 

2.  Chemicallv  etch  sample  with  A-B  etchant  as  in  Section  C.3, 

Step  1. 

3.  Check  depth  of  etching  with  photoluminescence  measurements. 

Repeat  Step  1  or  Step  2  until  active  layer  is  removed. 

4.  Remove  photoresist  from  mesas  with  the  plasma  line  etcher. 
Use  0t  at  maximum  power  for  20  min. 

C.  5  DBR  GRATING  APPLICATION* 

1.  Soak  clean  sample  in  warm  (>60°C)  trie lorethylene  for  2  to  5 
min. 

2.  Rinse  sample  in  dry  isopropyl  alcohol  for  l  to  2  min. 


Omit  this  section  if  active-passive  devices  are  desired. 
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3.  Blow  dry  sample  with  nitrogen 


4.  Dehydration  hake  sample  at  170°C  for  30  min. 

5.  Spray  deposit  a  1500-A-thick  layer  of  photoresist  (one  part  AZ 
1350B  diluted  with  four  parts  of  AZ  thinner)  over  entire 
sample.  Typical  system  parameters  are  a  spray  duration  of 
0.125  sec  and  a  spray  distance  of  11  cm  from  the  sample. 

6.  Oven  bake  sample  at  65°C  for  25  min. 

7.  Expose  the  sample  with  interfering  beams  of  the  4579  A  line  of 
an  argon  laser.  The  typical  exposure  time  is  about  30  sec  for 
a  combined  beam  intensity  of  10  mW/cm^.  10  mW/cm  corresponds 
to  an  indicated  voltage  of  0.370  from  the  solar  cell  mounted 
on  the  stage.  The  sample  is  aligned  by  reflecting  a  He-Ne 
beam  off  of  a  cleaved  face  back  to  the  HeNe  laser.  Each  of 
the  interfering  beams  should  be  of  equal  intensity  within 
10%.  A  standing-wave  pattern  can  be  observed  on  the  beam¬ 
splitting  prism.  By  slightly  adjusting  both  beam  angles,  the 
standing-wave  pattern  can  be  made  symmetrical,  indicating  both 
beam  angles  are  symmetric  with  respect  to  the  normal  to  the 
sample. 

8.  Develop  photoresist  for  1  sec  using  equal  parts  Shipley  MF  312 
developer  and  deionized  water.  Rinse  sample  in  dionized  1^0 
for  60  sec,  and  blow  dry. 

9.  Measure  the  grating  periodicity  using  the  ellipsometer .  Both 
first-order  Littrow  angles  are  measured  to  within  20  sec  of 
arc  and  averaged.  The  grating  orientation  with  respect  to  the 
crystal  cleave  is  simultaneously  measured  using  the  rotatable 
stage  mounted  on  the  ellipsometer.  The  latter  measurement  is 
accurate  to  within  1  min  of  arc. 

If  either  the  grating  period  or  grating  orientation  is  wrong, 
the  photoresist  can  be  removed  from  the  sample  with  acetone 
and  steps  1  through  9  repeated. 

10.  Ion  beam  etch  grating  into  the  GaAs.  Mount  sample  on  rotary 
stage  with  heat  conductive  Si02  paste.  Approximate  settings 
are  described  in  Section  C.3,  Step  15. 

OXIDATION  AND  STRIPE  CONTACT  DEFINITION 

1.  Soak  clean  sample  in  warm  (<60°C)  trichloroethylene  for  2  to  5 
min . 

2.  Rinse  sample  in  dry  isopropyl  alcohol  for  1  to  2  min. 
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3.  Blow  dry  sample  with  nitrogen. 

4.  Dehydration  bake  sample  at  170°C  for  30  min. 

5.  Deposit  3000  A  thick  Si02  layer  over  entire  sample  surface. 

6.  Repeat  Steps  1  through  4. 

7.  Mount  sample  on  glass  slide  with  Apiezon  wax. 

8.  Cool  and  then  whirl  clean  sample  with  a  1:1  mixture  of  HMDS 
and  AZ  thinner. 

9.  Apply  unthinned  filtered  AZ  1350J  to  sample. 

10.  Spincoat  sample  at  5000  rpm  for  30  sec. 

11.  Air  dry  sample  for  15  min. 

12.  Oven  bake  sample  at  90°C  far  25  min. 

13.  Use  mask  0E-1  with  mask  alignment  instrument  and  expose  sample 
for  16  sec.  Mask  GE-1  consists  of  a  five-row  matrix  of  25-  by 
450-um  windows. 

14.  Develop  photoresist  for  1  min  using  equal  parts  Shipley  MF  312 
developer  and  deionized  water.  Rinse  sample  in  deionized 
water  for  60  sec  and  blow  dry  with  nitrogen. 

15.  Oven  bake  sample  at  90°C  for  20  min. 

16.  Etch  windows  (25  by  450  pm)  through  the  SiCK  using  the 

plasmaline  etcher  at  maximum  power  with  CF^  for  about  20  min. 

C.7.  METALIZAT10N  AND  MOUNTING  OF  INDIVIDUAL  DEVICES 

1.  Clean  sample  in  a  boiling  mixture  of  acetone  and  xylene  (1:1) 

for  3  min.  Rinse  sample  in  three  changes  of  hot  xylene,  and 

blow  dry  with  nitrogen. 

2.  Mount  sample  face  down  on  plating  fixture  with  conductive 

(Dotite)  paint.  Coat  all  conducting  areas  of  the  fixture  and 
the  sides  of  the  sample  with  Apiezon  solution  (Apiezon  wax  in 
xylene).  Let  the  Apiezon  solution  dry. 

3.  Immerse  sample  in  HC1:H20  (1:1)  for  3  min.  Do  not  rinse. 

Transfer  sample  directly  to  the  Pt  plating  solution.  The  Pt 
plating  solution  is  a  commercial  halite  solution. 
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A  current  density  of  200  pA/mil  has  been  established  as  being 
an  optinum  current  density  for  platinum  plating  in  this 
procedure.  The  desired  thickness  of  platinum  is  between  0.2 
and  0.3  pm,  which  takes  approximately  two  min. 

4.  Overplate  the  Pt  layer  with  Au  using  a  ^standard  cyanide 

solution.  A  current  density  of  2  pamp/mil  typically  gives 
a  0. 5-pm-thick  layer  in  about  1  min. 

5.  Repeat  Steps  2  with  the  face  of  the  sample  up. 

6.  Repeat  Step  3. 

7.  Use  a  Pt  plating  current  of  166  pA  per  window  for 

approximately  1  min  to  obtain  a  0.2  to  0.  3-pm-thick  layer  of 
Pt. 

8.  Overplate  the  Pt  layer  with  a  0.5-pm  layer  of  Au.  A  plating 
current  of  16  pA  per  window  for  2  min  is  required. 

9.  Demount  sample,  clean  with  hot  xylene,  and  blow  dry  with 

nitrogen. 

10.  Cleave  and  scribe  individual  lasers  from  the  sample. 

11.  Scrub  botton  (n-side)  of  the  laser  to  a  10  to  18  transistor 
header  using  a  AuSi  eutectic  preform  at  3b5°C. 

12.  A  1-mil  Au  wire  is  thermal  compression  ball  bonded  to 

the  25  by  450-pm  Pt-Au  window.  The  other  end  of  the  gold  wire 
is  thermal  compression  ball  bonded  to  the  header  post. 

TEST  LASERS 


1.  Check  current-voltage  characteristics  on  the  Tektronix  by  573 
transistor  curve  tracer. 

2.  Check  luminescence  from  diode. 

3.  Check  spontaneous  and  lasing  spectra. 
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APPENDIX  D 


ELECTROCHEMICAL  DEPOSITION  OF  OHMIC  CONTACTS 
FOR  DIODE  LASERS 


The  typical  method  used  for  making  ohmic  contacts  on  GaAs  high- 
power  devices  was  to  resistively  evaporate  AuGe  for  the  back  (n)  contact,  and 
Cr  and  An  through  a  physical  mask  for  the  front  (p)  contact,  followed  by 
sintering  (alloying)  the  contacts  at  an  elevated  temperature.  These  contacts 
were  not  reproducible  and  had  a  fairly  high  impedance.  In  addition,  the 
metallization  for  the  n-contact  must  be  wet  by  the  AuSi  used  to  attach  the 
device  to  the  package,  and  the  front  p-contact  must  be  suitable  for 
thermocompression  wire  bonding. 

The  resistively  evaporated  contacts  proved  to  be  of  limited 
usefulness  for  small  geometries  evaporated  through  the  physical  mask  (for  the 
p-contact)  because  of  difficulties  encountered  in  the  alignment  of  the  mask  to 
the  wafer  and  also  the  edge  definition.  Photolithography  and  chemical  etching 
of  the  front  p-contact  was  attempted  but  was  found  to  be  impractical  because 
of  the  surface  attachment  of  the  GaAs  by  the  etches. 

At  the  suggestion  of  A.  B.  Chase,  electroplated  metallization  of  a 
Pt  and  Au  ohmic  contact  was  tried. 

The  electroplating  procedure  developed  involves  the  application  of 
a  given  contac  geometry  in  photoresist  (Shipley  AZ1350)  for  the  p-side 
contact.  This  p.  u-oresist  mask  has  a  two-fold  purpose:  (1)  to  provide  a  mask 
for  the  confinement  of  the  surface  area  to  be  plated  and  (2)  to  provide 
sidewall  confinement  of  the  plated  material.  At  this  point,  the  wafer  is 
attached  to  the  A1  strip  with  conductive  paint  (Dotite),  and  all  surfaces 
except  those  to  be  plated  are  insulated  with  Apiezon  wax  and  xylene  solution. 

A  very  critical  consideration  of  successful  electrochemical 
deposition  is  that  the  surface  be  free  of  any  foreign  or  native  oxides,  which 
appears  to  be  satisfied  by  a  predeposition  immersion  of  the  wafer  in  a  1:2 


solution  of  HCLtDit^O  for  3  min  at  ambient  temperature.  The  wafer  is  then 
transferred  directly  to  the  Pt  plating  solution  (the  Pt  plating  solution  is  a 
commercial  Halite  solution)  without  rinsing  to  avoid  any  possibility  of 
oxidation. 

2 

A  current  density  of  200  pA/mil  has  been  established  as  being  an 
optimum  current  density  for  Pt  plating  in  this  procedure.  Voltage  should  be 
applied  prior  to  insertion  of  the  wafer  into  the  flowing  plating  solution,  to 
avoid  an  alteration  of  the  GaAs  surface.  A  ph  of  3.5  was  chosen  for  the  Pt 
plating  solution  to  prevent  both  surface  damage  and  undercutting  of  the 
masking  material.  Although  Pt  plating  is  typically  performed  at  ambient 
temperatures,  an  elevated  temperature  plating  solution  produces  superior- 
quality  depositions,  as  well  as  an  accelerated  deposition  rate.  Routinely 
successful  plating  was  accomplished  at  a  constant  temperature  of  ”60  C  for 
approximately  2  min.  The  thickness  of  the  deposited  film  was  between  0.2  and 
0.3  pm.  Films  thicker  than  this  lifted  and  pulled  GaAs  from  the  bulk  crystal 
because  of  differences  in  the  thermoexpansions  of  the  materials. 

At  this  time,  the  wafer  is  ready  for  the  Au  overplating  (plating 

solution  being  a  standard  cyanide  solution),  which  is  conducted  the  same  as 

2 

the  Pt  deposition  except  for  the  current  density  of  2  pA/mil  ,  ambient 
temperature  plating  solution,  a  film  thickness  of  0.6  to  0.7  pm,  and  an 
exclusion  of  the  HC1  surface  preparation. 

The  n-side  plating  is  changed  only  by  a  reduction  in  the  Au 
thickness  to  0.2  to  0.3  pm.  The  mounting,  surface  preparation,  current 
density,  and  temperatures  are  the  same  as  the  p-side. 

The  results  we  have  obtained  both  physically  and  electrically  from 
this  plating  procedure  have  been  excellent.  Physical  adhesion  to  both  p-  and 
n-sides  has  been  far  superior  to  any  other  contact  previously  used,  thus 
greatly  increasing  contact  integrity  between  p-contact  and  the 
thermocompression  bonded  wire  contact,  as  well  as  the  n-contact  and  the 
package. 

The  results  obtained  electrically  indicate  that  the  contacts  are 
able  to  withstand  repeated  current  pulses  greater  than  100,000  amp/cm4  without 
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noticeable  degradation.  The  contacts  are  also  of  very  low  resistance,  on  the 

“8  2 

order  of  approximately  6  x  10  Q/cm  . 

The  areas  of  further  study  should  include  the  establishment  of 
current-time  values  and  thickness  parameters  to  obtain  a  given  thickness  of 
film  on  a  given  geometry  and  also  to  increase  our  knowledge  of  the  ohmic 
contact  itself. 
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APPENDIX  E 


COUPLED  LASER  THEORY 


The  interaction  of  two  lasers  was  investigated  for  optical 

coupling  and  phase  locking.  Previously,  Adler^  examined  phase  locking  for  two 

2 

oscillators  in  tandem  and  Basov  et  al.,  examined  phase  locking  for  two 

oscillators  in  parallel  or  side  by  side.  In  the  analysis  described  here,  the 
two  cases  are  combined  into  a  more  general  theory  and  the  work  of  Basov  et  al. 
is  extended.  For  the  case  when  the  lasers  are  placed  in  parallel,  experi¬ 

mental  results  are  given  and  compared  to  the  theory. 

E.l.  DESCRIPTION  OF  PHASE  LOCKING 

Lasers  have  a  characteristic  or  natural  frequency  when  they  are 

free,  not  coupled  to  any  other  lasers.  When  two  lasers  are  coupled,  their 
output  changes.  Some  possible  ways  of  coupling  are  by  end  coupling  (when  two 
lasers  are  placed  In  tandem),  diffraction  or  directional  coupling  in  lasers 
placed  stde-by-si ie ,  and  backscattering  in  ring  lasers  (Fig.  E-l).  The  case 
for  lasers  placed  side  by  side  in  parallel  is  treated  here. 

Very  weak  coupling  between  lasers  introduces  no  change  in  the 
fields  from  that  of  free  lasers.  In  this  case,  the  lasers  can  be  considered 
free.  For  stronger  coupling,  changes  are  seen  in  the  output  spectrum  (Fig. 
E-2).  A  new  frequency  appears  along  with  the  characteristic  frequency.  This 
new  frequency  is  the  same  as  that  of  the  coupled  laser.  We  shall  call  this 
regime  of  weak  coupling  "optical  coupling."  The  fields  of  the  lasers  are 
incoherent  with  respect  to  each  other  because  the  relative  phase  shift  varies 
with  time.  Thus,  the  fields  will  not  produce  interference  fringes  in  a 
diffraction  pattern. 


*R.  Adler,  Proc.  IEEE  61,  1380  (1973). 

~N.  G.  Basov,  E.  M.  Belenov,  and  V.  S.  Letokhov,  "Diffraction  Synchronization 
of  Lasers,”  Sov.  Phys.  -  Tech.  Phys.  10,  845  (1965). 
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Fig.  E-I.  Geometries  for  Coupled  Laser.  (a)  Tandem, 
end  coupling.  (b)  Parallel,  dif fractional 
or  directional  coupling.  (c)  Ring, 
scattering. 
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Lg.  E-2.  Laser  Spectra.  (a)  Free  uncoupled, 
(b)  Weakly  coupled. 


I>t 


\ 


\ 


The  intermediate  coupling  region  covers  the  range  of  coupling 
strengths  between  the  weak  coupling  regime  and  the  threhold  for  locking. 
Sideband  frequencies  appear,  with  the  intensity  and  spacing  determined  by  the 
coupled  fields  and  the  initial  frequency  shift.  In  general,  two  free  lasers 
will  not  oscillate  exactly  at  the  same  frequency  so  there  is  an  initial 
frequency  difference.  If  this  frequency  shift  is  small  and  the  coupled  field 
is  large,  many  sideband  frequencies  appear  with  a  small  spacing  in  frequency 
(Fig.  E-3).  Ultimately,  for  large  coupling,  the  frequency  spacing  is  zero, 
and  the  lasers  are  locked.  The  opposite  is  true  for  large  frequency  shifts 
and  small  coupled  fields.  The  overall  effect  is  a  pulling  of  the 
characteristic  frequency  toward  a  frequency  that  is  between  the  characteristic 
frequencies  of  the  lasers  and  is  called  the  locking  frequency.  The  locking 
frequency  is  determined  by  the  relative  strengths  of  the  laser  fields.  As  the 
coupling  strength  approaches  the  threshold  for  locking,  the  coherence  length 
of  the  lasers  becomes  greater.  Interference  patterns  appear  in  diffraction 
patterns  and  become  more  distinct  as  the  locking  threshold  is  approached. 

The  strong  coupling  region  covers  locked  lasers.  As  in  the  case 
for  intermediate  coupling,  the  degree  of  locking  is  also  determined  by  the 
initial  frequency  difference.  Locked  lasers  oscillate  at  a  common  frequency, 
the  locking  frequency  (Fig.  E-4).  Because  they  have  a  common  frequency,  the 
relative  phase  shift  between  the  lasers  is  constant  in  time.  The  lasers  are 
thus  coherent  with  each  other  and  exhibit  interference  fringes  in  their 
diffraction  pattern.  Note  that,  although  the  relative  phase  shift  between  the 
lasers  is  constant,  it  is  not  zero. 

These  coupling  regions  are  not  rigidly  defined,  and  the  transitions 
between  the  regions  are  gradual.  However,  from  this  discussion,  the 
determination  of  whether  or  not  lasers  are  locked  can  be  made  by  looking  at 
the  spectrum  of  the  lasers,  or  at  their  diffraction  pattern. 

E-2.  THEORY  OF  OPTICAL  COUPLING 

Optical  coupling  is  the  weak  coupling  limit  of  the  more  general 
problem  of  phase  locking.  Optical  coupling  occurs  when  the  field  from  laser  1 
is  coupled,  by  some  means,  to  Laser  2.  The  output  frequencies  of  Laser  2  is 
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Fig.  E-3.  Intermediately  Coupled  Lasers  with 
Stronger  Coupling 


E-5 


v 


LASER  1 


LASER  2 


1 


\ 


i 


that  of  the  coupled  field  of  Laser  1  and  the  characteristic  frequencies  of  a 
free  and  uncoupled  Laser  2.  In  general,  optical  coupling  is  a  mutual  effect 
so  that  the  outputs  of  both  lasers  are  different  from  the  outputs  of  a  free, 
uncoupled  laser.  However,  the  field  of  one  laser  may  be  much  greater  than  the 
field  of  the  other  laser,  in  which  case  only  one-way  coupling  need  be 
considered. 


Optical  coupling  can  be  described  by  the  regenerative  amplifier 
model.  A  laser  below  the  lasing  threshold  is  a  regenerative  amplifier.  In 
this  case,  although  a  laser  may  be  lasing  at  its  characteristic  frequency,  it 
does  not  lase,  but  only  amplifies  the  coupled  field  frequency. 

In  developing  the  regenerative  amplifier  model,  the  equations  for 
an  amplifier  are  derived.  It  can  be  shown  that  a  simple  amplifier  model  is 
not  sufficient  to  describe  the  experimental  data.  However,  when  feedback  is 
included  in  the  model  to  form  a  regenerative  amplifier,  the  model  compares 
favorably  with  the  data. 


The  two-laser  geometry  is  shown  in  Fig.  E-5.  For  clarity,  one-way 
coupling  only  is  shown.  It  is  assumed  that  the  coupled  frequence  is  within 
the  gain  envelope  of  the  laser. 


If  the  net  gain  of  Laser  2  at  frequency  w^  is  small  and  the 
feedback  minimal,  it  is  an  amplifier  at  that  frequency.  The  differential 
equation  for  an  amplifier  with  distributed  coupling  and  gain  is 


dE 

^r  =  gE2+KE1  (E-l) 

where 

E2  =  field  of  the  amplifier  (Laser  2) 

E|  =  coupled  field  of  the  laser  (Laser  1) 
g  =  gain-loss  coefficient 

K  =  field  coupling  coefficient  between  the  laser  and  the 
amplifier 
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The  change  in  the  field  of  the  amplifier  is  the  result  of  the  net  gain  of  the 
amplifying  material  plus  the  coupled  field  from  Laser  1. 


If  no  light  was  injected  initially  (by  end  coupling),  then 


E2(x  =  0)  =  0 


The  field  in  the  amplifier  as  a  function  of  propagation  distance  is  then 


(E-2) 


K 0 ( x )  =  k  [(exp  (gL)  -  ig)]  E 


(E-3) 


(E-4) 


where  L  =  amplifier  length. 

During  one  pass  through  an  amplifier  with  no  end  coupling,  the 
coupled  field  is  amplified  bv  a  factor  of  y.  For  an  initially  injected  field 

°f  E.n, 


E?(x)  =  E^exp  (gL  +  i9)  +  yE  j 


where  0  is  the  phase  shift  incurred  while  propagating  through  the  amplifier. 

An  amplifier  with  feedback  is  a  regenerative  amplifier.  From  Eq . 
(E-5)  and  Fig.  K-6 


'out  =  KineXp  (rL  +  i9)  +  YE1 


(K-6) 


The  feedback  equation  is 


E in  =  [r  EQ(it  exp  (gL  +  i9)  +  yE^r 


(E-7  ) 


The  output  of  the  amplifier,  K7,  is 


E  =  1  -  r  E 
l  out 


0>8) 


► 


Fig.  E-6 .  Regenerative  Amplifier  < 
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where  r  is  the  field  reflection  coefficient.  By  solving  Eqs.  (E-6),  (E-7), 

and  (E-8),  the  equation  for  the  intensity  ratio  is  obtained. 
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gL  , 


(E-9) 


where 


Ij  =  intensity  of  the  laser 


19  =  intensity  of  the  amplifier 
E  =  length  of  the  amplifier 
R  =  field  reflection  coefficient 
g  =  gain-loss  per  unit  length 
K  =  coupling  coefficient  per  unit  length 

9  =  phase  shift  resulting  from  propagation  through  the  amplifier 

The  intensity  ratio  is  proportional  to  the  coupling  coefficient  squared.  The 
frequency  dependence  of  the  intensity  ratio  is  contained  in  the  gain  g.  The 
gain  distribution  is  given  in  the  following  section.  Equation  (E-9)  is  valid 
for  gains  less  than  the  threshold  for  lasing.  As  the  threshold  is  approached, 
the  gain  starts  to  saturate  and  at  threshold  the  intensity  ratio  is  actually 
finite  instead  of  blowing  up  as  predicted  by  Eq .  (E-9).  For  the  data 
presented  in  the  following  section,  all  gains  are  below  threshold. 

The  threshold  gain  can  be  found  from  the  condition  for  lasing 


1  =  r  e 


gE 


(E-10) 


Hence,  the  threshold  gain  is 
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This  gain  is  the  net  gain,  the  gain  produced  by  pumping  minus  any  losses 

resulting  from  scattering  or  absorption.  A  typical  threshold  gain  is  gt  =  14 

-1 

cm 

E.2.1  Gain 

The  gain  for  a  given  pumping  level  j  can  be  approximated  by  a 
Boltzmann-Lorentzian  distribution 


g  .(v)  =  g  exp(-ghv) 
1  o 


(E-12) 


where 

v^=  frequency  of  the  peak  gain 

6v  =  half-width  of  the  gain  profile  at  half  maximum 
g  =  1/kT  =  1/38  eV  at  room  temperature 

The  gain  coefficient  gQ  can  be  determined  by  finding  a  longitudinal  mode  at 
threshold  and  comparing  it  with  the  mode  exhibiting  the  peak  gain. 
Rewriting  g.(v)  in  terms  of  wavelength, 

1  +  (AXt/<5A  .) 

1  +  (aX/SXj)2 


g.(X)  =  gt  exp 


124003 


AA 

A  2 
o 


(E-13) 


where 

A  =  threshold  wavelength 
Aq  =  wavelength  of  peak  gain 
AA  =  A  -  A 


E-12 


<5A  =  half-width  of  gain  profile 
gt  =  threshold  gain 

Experimental  data  are  given  in  Section  E.3  and  compared  with  this  theory. 

F, .  3  OPTICAL  COUPLING  DATA 

Two  double-heterostructure  CaAlAs  lasers  were  placed  side  by  side 
as  shown  in  Fig.  E-5.  The  lasers  were  defined  by  SiO^  stripes  laid  on  the 
wafer.  A  cross  section  of  the  laser  geometry  is  shown  in  Fig.  E-7.  The  top 
GaAs  layer  provided  .,n  ohmic  contact  between  the  p-type  GaAlAs  and  the  metal 
contact.  Because  the  metal  contact  was  common  to  all  lasers,  all  of  the 
lasers  were  simultaneously  pumped  in  parallel. 

The  injected  current  pulse  amplitude  was  slowly  increased  until  one 
of  the  diodes  started  lasing  (Diode  1).  Measurements  of  the  light  outputs 
from  the  lasing  diode  and  the  nonlasing  diode  (Diode  2)  were  made  at  various 
currents  between  the  thresholds  of  Diode  1  and  Diode  2.  Spectra  of  each  of 
the  diodes  were  monitored  with  a  resolution  of  0.3  A.  In  addition,  the  image 
of  the  diodes  could  be  translated  across  the  monocrometer  slit  to  obtain 
spatial  near  or  far-field  profiles  at  a  given  frequency. 

Typical  spectra  are  shown  in  Fig.  E-8.  The  periodicity  was  40  pm, 
with  a  stripe  width  of  15  pm.  The  longitudinal  mode  envelope  of  Diode  1 
ranged  from  0.8700  pm  to  beyond  ('.8725  pm  peaking  at  A  =  0.8723  pm.  Diode  2 
had  two  sets  of  mode  envelopes  oscillating.  The  shorter  wavelength  modes 
peaked  at  0.8707  pm  and  died  out  at  0.8711  pm.  These  modes  were  the  natural 
modes  of  Diode  2.  The  other  set  of  modes  at  the  longer  wavelength 
(A  ~  0.8720  pm)  was  beyond  the  natural  mode  envelope  of  Diode  2.  These  modes 
were  the  result  of  amplified  optical  coupling  from  Diode  1.  Matching  the 
spectrum  of  Diode  1  with  the  spectrum  of  Diode  2  indicates  that  the  coupled 
modes  of  Diode  2  are  at  the  same  frequency  as  that  of  Diode  1.  Note  that  the 
optical  coupling  efficiency  without  amplification  is  about  1%.  This  would 
make  coupled  light  unobservable.  But  with  amplification  by  a  gain  medium, 
optical  coupling  can  be  observed. 
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Fig.  E-7.  Coupled  Laser  Stripes  Cross-Sectional  Geometry 


These  experimental  results  for  optical  coupling  are  compared  with 
the  theory  in  Section  E.2.  The  gain  of  the  amplifier  is  given  by  Eq.  (E-13). 
With  n  =  3.6,  R  =  (n  -  1/n  +  l)2,  1  =  400  pm,  6  =  1/38  eV,  X^  =  8611  A, 
X  -  8607  A,  AX^  =  4  A,  and  6X  =  6  A,  Eq.  (E-13)  becomes 


,  ,s  1.44 

g  (X)  =  g  - 2 - 

J  C  1  +  (AX) /  36 


(E-14) 


At  the  coupled  wavelength  X  =  8618  A,  the  gain  threshold  for  lasing  [from 
Eq.  (E-l l ) ]  is  gt  =  14.3  cm~*.  In  this  case,  Eq.  (E-14)  is 


g  .(X)  =  4.7  cm 


(E-15) 


Working  backwards  from  the  measured  intensity  ratio,  I 2 / I j  =  0.2, 
to  predict  the  coupling  coefficient  for  a  gain  of  4.7  cm  ' ,  and  Eq.  (E-9) 
predicts 


|  KL| 2  =  0.023 


(E— 16) 


Two  and  a  half  percent  of  the  intensity  of  Diode  1  was  coupled  to  the 
amplifier. 

A  theoretical  expression  for  the  coupling  coefficient  of  the  fields 
in  the  strong  coupling  limit  is 


(Ac  w)2k3 

|K|  = - lH1  exp 


j^-(Ae  w)k2sj 


(E-l 7 ) 
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where 


Ac  -  dielectric  constant  step  between  the  diodes  and  t  tie 
intervening  areas 

w  =  width  of  the  laser  stripes 

s  =  separation  between  the  near  edges  of  the  diode  stripes 
Fig.  K-7 

k  =  2/A 

o  o 

n  =  refractive  index  in  the  diodes 

K  =  coupling  coefficient  between  fields  in  the  two  diodes 
Maximizing  K  with  respect  to  (Ac  w),  the  maximum  value  of  K  is 

(K)max  =  —  — --nr  -  ( E—  1 8 ) 
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For  a  tvpical  diode  laser,  X  =  0.862  um,  s  =  25  pm,  n  =  3.6,  L  =  400  pm,  and 
(KL)  =  0.33  cm'*.  These  values  permit  the  calculation  of  the  fraction  of 

laser  light  that  is  coupled  per  pass: 

(KL)2  =  1.3% 


which  compares  with  the  value  determined  from  the  measurements  of  2.5%. 

Optical  coupling  has  been  observed,  and  the  measured  coupling 
coefficient  compares  favorably  with  the  theoretical  value  when  the 
regenerative  amplifier  model  is  used.  In  this  model,  it  is  assumed  that  the 
'■unli  1  field  is  within  the  gain  envelope  of  the  diode  to  which  it  is  coupled, 


and  that  secondary  effects  such  as  coupling  of  the  coupled  field  bad:  into  the 
original  diode  are  negligible. 


R.4  THEORY  OF  1NJFXTI0N  LOCKING 

The  theory  of  injection  locking  has  been  discussed  in  E.2  for  the 
case  when  one  oscillator  is  weakly  coupled  to  the  other.  The  behavior  of 
coupled  oscillators  near  and  above  the  threshold  for  phase-locking  are 
described  here.  While  the  model  of  the  regenerative  amplifier  works  well  for 
weakly  coupled  oscillators,  it  does  not  explain  the  phenomena  of  strongly 
coupled  oscillators.  The  regenerative  amplifier  model  does  not  permit  the 
pulling  of  frequencies,  the  locking  of  lasers  at  a  common  frequency,  or  any  of 
the  other  locking  characteristics  such  as  locking  times. 

Three  models  have  been  proposed  to  describe  phase  locking.  The 
first  was  the  phasor  model  of  Adler. ^  Adler's  model  dealt  with  the  relative 
phase  difference  between  two  oscillators.  His  model  assumes  a  master-slave 
arrangement,  that  coupling  occurs  in  only  one  direction.  This  model  is  good 
when  there  is  an  optical  isolator  between  lasers  or  when  one  laser  is  much 
more  intense  than  the  other.  However,  when  the  laser  fields  are  of  comparahle 
magnitude,  Adler's  model  breaks  down,  since  the  model  does  not  allow  for 
mutual  coupling. 

q 

The  second  model  is  that  of  Statz  and  Tang.  Their  model,  again  a 
master-slave,  is  based  on  oscillator  equations  at  two  frequencies;  the  free, 
uncoupled  frequency  of  the  slave  and  the  coupled  frequency  of  the  master. 
This  model  breaks  down  when  the  master-slave  relationship  does  not  hold,  and 
does  not  demonstrate  frequency  pulling  and  the  dynamic  effects  of  locking. 

2 

The  third  model  was  developed  by  Basov,  Belenov,  and  Letokhov  and 
is  the  one  used  here.  This  model  is  based  on  an  oscillator  equation  for  each 
oscillator,  with  additional  terms  and  equations  to  describe  the  coupling. 


3 

H.  Statz  and  C.  L.  Tang,  Phase-Locking  of  Laser  Oscillators  by  Injected 
Signal,"  J.  Appl.  Phys.  38,  323  (1967). 
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In  Basov's  node l ,  diffraction  coupling  is  assumed,  instead  of  the  end  coupling 
of  Adler  and  Statz  and  Tang.  '  The  described  theory  here  includes  both  forms 
of  coupling.  Whereas  Basov  onlv  gave  the  conditions  for  locking,  other 
locking  phenomena  such  as  frequency  pulling,  locking  frequency,  and  locking 
tines  are  derived. 

K .-'*.1  Locking  Conditions 

The  oscillator  equations  for  the  lasers  are 


u>  i  #  0 

C1  +  IT  ei  "  woa  l2  +  “1  £1  =  0 

( E-20  ) 

..  u>^  #  0 

e2  +  if  ^2  '  “oa  f  +  u2  e2  =  ° 

(E-21) 

where 

€  .  = 

.1 

fie  lil  of  osc  i  1  lat  or  j 

Q  = 

cold  cavity  Q  of  the  oscillators 

CO  .  = 

J 

natural  frequency  of  j 

II 

c 

3 

locking  frequence  of  the  oscillators 

a  = 

coupling  efficiency 

It  is  assumed  that  the  cavitv  Os  for  both  oscillators  are  equal.  The  second 
tern  represents  the  losses  in  I  he  oscillator,  and  the  third  tern  represents 
the  coupled  radiation.  This  nav  he  provided  by  end  coupling  or  side  coupling 
(  lire,  t  i  >nal  coupling  and  diflrntion.il  coupling!.  It  these  oscillators  were 
unronp ' ,  s' ,  ''lev  would  oscillate  at  their  natural  frequeuev  «  and  o»., ,  respec- 

’  '  •  h  ■  ■  ‘  Mirth  tern'. 
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Since,  in  general,  the  fields  have  different  phase  shifts,  can 


be  rewritten. 


e  .  =  E  .exp  ito  (t  +  6  .) 
3  j  o  j 


(E-22) 


where  Ej  is  slowly  varying  over  an  optical  cycle  time  period,  l/2TrwQ.  The 
change  in  the  phase  shift  4  .  is  also  small  over  a  time  period.  Cold  cavity  Qs 

1  j 

of  lasers,  such  as  GaAs,  are  of  the  order  of  lCr  so  that  1/Q  «  1.  The 
coupling  efficiency  a  is  of  the  order  of  a  few  percent.  By  substituting 
Eqs.(E-22)  and  (E-23)  in^o  F.qs.  (E-20)  and  (E-21),  and  regrouping  terms 
according  to  their  order  of  magnitude,  an  equation  for  the  change  in  the 
electric  field  phase  shift  is  obtained. 
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where  j  *  i. 


For  the  frequencies  of  interest, 
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The  conditions  for  locking  are  obtained  by  examining  Eq.  (E-25)  in 


the  steadv  state  when  4  .  =  0 
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(E-26) 
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Equation  (E-26)  provides  the  frequency  range  over  which  locking  can  occur  for 
a  given  coupling  efficiency  and  oscillator  fields. 

The  coupling  can  be  provided  by  side  coupling  or  end  coupling.  For 
side  coupling,  the  coupling  process  can  be  directional  or  diffractional 

9 

coupling.  Basov-  assumed  diffractional  coupling  to  derive  his  results,  which 
agree  with  Eq.  (E-26).  For  directional  coupling  K  =  Xo/2un.  For  end 
coupling,  the  coupling  is  dependent  upon  the  reflectivity  of  the  end  mirrors 
or  the  cavity  Q.  Thus,  F,q.  (E-26),  in  terns  of  the  cavity  Q,  is 


( E— 27  ) 


whe-c 


X  a  —  1/1  InR 
o  o 


n  =  refractive  index 

I,  =  oscillator  length 

R  =  mirror  reflectivity 

a  =  distributive  loss  coefficient 
o 

"nuitiei  ( K-1'7)  is  identical  to  the  results  of  Adler'  and  Statz  and  Tang. 
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PHASE-LOCKING  DATA 


K.3 

Phase  locking  has  been  observed  in  a  pair  of  lasers  with  a  25-pm 
separation  distance  and  verified  hy  two  different  methods:  (1)  comparing  the 
longitudinal  node  spectrum  of  each  laser  and  (2)  observing  interference 
fringes  in  the  far-field  diffraction  pattern.  The  spectra  of  two  locked 
lasers  are  shown  in  Fig.  F,-9.  As  can  be  seen,  the  longitudinal  nodes  of  both 
lasers  peak  at  a  common  wavelength. 

Two  lasers,  when  phase  locked,  are  coherent  with  respect  to  each 
other  and,  hence,  exihibit  interference  fringes  in  the  far-field  diffraction 
pattern.  The  far-field  diffraction  pattern  for  phase-locked  lasers  is  shown 
in  Fig.  E-10.  Interference  between  the  two  lasers  is  clearly  evident.  From 
the  interference  pattern,  it  was  inferred  that  the  lasers  were  10  pm  apart, 
using  the  double  slit  model  for  plane  waves,  which  is  in  reasonable  agreement 
with  the  stripe  separation  of  13  pn. 

The  individual  free  laser  characteristics  were  not  measured  because 
the  lasers  were  not  individually  pulsed.  However,  the  range  of  wavelengths 
over  which  phase  locking  occurs  can  be  calculated.  The  coupling  coefficient 
can  he  determined  from  the  equation  for  the  directional  coupling  coefficient 
per  unit  Length.  For  X  =  0.86  pm,  L  =  400  pn,  and  n  =  3.6; 


K 


o 


2.82  x  10 


From  Hq.  (E-26),  the  locking  range  for  E^/E,  =  1  is 


AX  =  0. 12  A 


E-22 
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Thus,  for  a  pair  of  mutually  coupled  lasers,  twice  this  range,  i.e.,  0.24  A, 
is  the  maximum  wavelength  difference  between  the  free  and  uncoupled 
wavelengths  for  which  phase  locking  can  occur. 

Experiments  are  currently  under  way  in  which  individually  pulsed 
lasers  are  being  used.  This  should  provide  more  information  about  the 
individual  lasers  both  above  and  below  the  threshold  for  locking  since  the 
power  of  each  laser  can  be  independently  controlled. 
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APPENDIX  F 


DBR:  A  MULTIFUNCTION  INTEGRATED 

OPTICAL  DEVICE 


F. 1  INTRODUCTION 

The  properties  of  wide-angle,  integrated  optical  Bragg  deflectors 
in  which  slab-coupled  optical  waveguides  are  used  are  analyzed.  Specifically 
considered  is  the  interaction  which  occurs,  through  the  internediary  of  a 
periodic  waveguide  perturbation,  between  an  incident  waveguide  within  the  core 
region  of  the  structure  and  a  Bragg-deflected  beam  guided  within  the  slab 
region  of  the  structure.  The  deflection  efficiencies  and  far-field  deflected- 
beam  intensity  patterns  characteristic  of  this  distributed  Bragg  deflector 
(DBD),  are  derived  for  both  transverse  electric  (TE)  and  transverse  magnetic 
(TM)-polarized  incident  waves  and  for  deflection  angles  between  ir/4  and  3ir/4 
rad. 

A  large,  number  of  passive  and  active  integrated  optical  devices  may 
he  developed  by  exploiting  tire  phenomenon  of  wide-angle  intrawaveguide  Bragg 
scattering  off  fixed  or  tunable  periodic  waveguide  perturbations.  The  prin¬ 
cipal  advantages  of  such  devices,  compared  to  their  shallow-scattering-angle 
counterpart s ,  »*•* J >  include  compactness,  good  incident-  and  deflected-bean 
separation  and  accessibility,  and  considerable  device  versatility  and  design 
flexibility.  Several  wide-angle  intrawaveguide  deflection  devices  have  been 
demonstrated,  which  have  relied  upon  optical  waveguides  overmoded  in  either  or 
both  dimensions  transverse  to  the  direction  of  propagation.  Multimode  devices 
that  tiave  depended  vipon  fixed  periodic  waveguide  perturbations  as  the  scatter¬ 
ing  internediary  have  included  thin-film  modulators , ^ ^  deflectors, ^  and  bean 
splitters. ^  Tunable  periodic  waveguide  deflectors  formed  by  magnetic  stripe 
domains  have  been  used  as  nulti-mode  waveguide  modulators  and  suggested  as 
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promising  wide-angle  intrawaveguide  bean  scanners.  It  is  the  purpose  here  to 
describe  the  operational  characteristics  and  potential  applications  of  wide- 
angle  f nt rawavoguide  Bragg  deflection  devices  that  utilize  single-node 
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two-dimensional  optical  waveguides  of  the  type  one  would  expect  to  be  used  in 
future  integrated  optical  circuits. 

The  single-mode  waveguiding  structure  has  been  generically  termed  a 

q 

slab-coupled  waveguide  by  Marcatili.  In  the  present  context,  the  usefulness 
is  that  it  may  support  both  one  and  two  dimensional  guided  waves,  the  former 
within  the  slab  region  of  the  structure,  and  the  latter  within  its  core.  By 
locating  a  periodic  waveguide  perturbation  of  appropriate  period  and  orienta¬ 
tion  within  or  adjacent  to  the  core  region,  a  two  dimensionally  confined  wave 
initially  propagating  within  the  core  region  may  be  gradually  converted  by 
means  of  coherent  Bragg  scattering  into  a  one  dimensionally  confined  mode  that 
propagates  within  the  slab  region.  The  resultant  periodically  perturbed  slab- 
coupled  waveguiding  structure  is  a  distributed  Bragg  deflector  (DBD). 

Passive  integrated  optical  applications  of  the  DBD  structure 
include  guided-wave  beam  deflection,  division,  and,  for  90  deg  Bragg  scatter¬ 
ing,  polarization  analysis.  In  addition,  because  of  the  gradual  or  distri¬ 
buted  nature  of  the  Bragg  scattering  process,  a  properly  configured  DBD  may  be 
used  to  effect  lensless,  intrawaveguide  beam  expansion.  When  combined  with  a 
suitably  oriented  output  grating  coupler,  the  latter  application  indicates 
that  the  utility  of  the  DBD  as  a  device  for  reducing  the  beam  divergence  of 
fully  integrated  lasing  sources  (provided  with  DBR  or  distributed  feedback 
cavities  and  mated  with  a  passive  waveguiding  region  in  which  the  beam 
expansion  could  be  accomplished),  or  as  a  means  by  which  two-dimensional  slab- 
coupled  waveguides  could  be  accessed  by  large-diameter  beams  externally 
incident  on  the  waveguiding  structure.  With  regard  to  active  integrated 
optical  devices,  the  relative  insensitivity  of  the  Bragg  scattering  condition 
within  the  DBD,  appears  to  limit  its  active  applications  to  those  that  may  be 
realized  through  the  direct  manipulation  of  its  spatial  (i.e.,  period  or 
perturbation  orientation)  properties.  As  a  consequence,  wide-angle  DBD  beam 
scanners  or  switches  based  upon,  for  example,  tunable  magnetic  stripe  domains 
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appear  feasible,  whereas  wide-angle  DBD  devices  dependent  upon  a  modification 
or  violation  of  the  Bragg  condition  do  not."7 
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F. 2  DEVICE  GEOMETRY  AND  ANALYSIS  OVERVIEW 

Top,  side  cross-sectional,  and  end  cross-sectional  views  of  the 
prototypical  DBD  structure  to  be  analyzed  are  shown  in  Fig.  F-l.  With 
reference  to  the  core  region,  in  which  two-dimensional  guiding  takes  place, 
mode  confinement  in  the  z-direction  is  assumed  to  occur  as  a  result  of  step 
discontinuities  in  refractive  index  at  the  superstrate  guide  and  guide 
substrate  interfaces.  Mode  confinement  within  the  core  region  in  the  y- 
direction  is  assumed  to  result  from  the  dispersive  effects  associated  with 
either  a  height  variation  in  the  vicinity  of  the  core  region  or  superstrate 
loading  over  the  core  region  (Fig.  F-2).  Figure  F-2(a)  is  a  cross-sectional 
view  through  an  imbedded  ridge^  DBD  structure  (lateral  confinement  caused  by 
height  variation),  and  Fig.  F-2(b)  is  a  cross-sectional  view  through  a  strip- 
loaded^  DBD  structure  (lateral  confinement  caused  by  superstrate  loading). 
For  analysis,  the  lateral  confinement  effect,  induced  either  by  a  height 
variation  or  by  superstrate  loading,  is  accounted  for  by  the  geometry  of  Fig. 
F-l(c)  in  which  the  core  region  has  been  assigned  a.  effective  refractive 
index,  n2  and  a  height  (in  the  z-direction)  equal  to  that  of  the  slab.  n->  is, 
of  course,  taken  to  be  greater  than  n,.,  the  effective  refractive  index  of  a 
one-dimensi onally  guided  mode  propagating  in  the  slab  region  of  the  structure. 

Coupling  between  the  core-propagating  and  slab-propagating  nodes  is 
assumed  to  occur  as  a  result  of  a  suitably  located  and  oriented  periodic  wave¬ 
guide  perturbation.  Although  this  analysis  is  applicable  to  any  type  of  per¬ 
turbation  geometry,  only  one  specific  case  will  be  considered:  the  square- 

wave  type  surface  corrugation  shown  in  Figs.  F-l(b)  and  F-l(c).  The  extent  of 
this  perturbation  in  the  z-direction  will  be  governed  by  the  depth  of  the 
surface  corrugations;  its  extent  in  the  y-direction  will,  for  convenience,  be 
taken  equal  to  the  equivalent  width1'1  (in  the  y-direction)  of  the  core- 
propagating  mode.  This  approximation,  illustrated  by  the  mode  profile  shown 
in  Fig.  F-l(c),  m,i”  be  exported  to  change  the  results  of  the  analysis  only 
slight lv. 

The  Brag:',  dr1].-.  ■■  ;.r.».  ess  to  be  analyzed  is  schematically  por- 
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t  raved  i  n  F  i  g.  F-l  ( ,.  > 


A  core-propagating  wave 


n,  Isuperstrate) 


Fig.  F-2.  DBD.  (a)  F,nd  cross-sectional  view  of 

imbedded  ridge.  (b)  End  cross-sectional 
view  of  strip-loaded  DBD. 


incident  on  the  distributed  Bragg  deflector  from  the  negative  x-direction  is 
converted,  by  a  periodic  waveguide  perturbation,  into  a  slab-propagating  wave 
that  is  deflected  through  an  angle  9„  with  respect  to  the  x-axis.  The  problem 

a 

addressed  is  that  of  finding  the  field  distribution  of  the  deflected  wave  . 
Coherent  waveguide  scattering  or  mode  conversion  problems  such  as  this  may  be 
treated  with  either  the  coupled-mode  formalism*'*’***  or  a  Green's  function 
technique.*"’  Given  the  interaction  geometry  shown  in  Fig.  F-l,  wherein  an 
insignificant  amount  of  the  deflected  wave  is  expected  to  be  fed  back  into  the 
incident  wave  and  for  which  the  incident  wave  field  distribution  within  the 
interaction  region  may  be  well  approximated,  the  Green's  function  technique  is 
preferred  and  will,  as  a  consequence,  be  the  one  used.  Accordingly,  the 
analytical  approach  to  be  followed  will  be  to:  (1)  define  an  appropriate 

Green's  function,  (2)  assume  a  form  for  the  core-propagating  wave  that 
accounts  for  the  Bragg  scattering  process  through  the  inclusion  of  an  attenua¬ 
tion  factor  exp(-ax),  (3)  calculate  the  deflected  wave  field  distribution  by 
convolving  the  incident  wave  field  distribution  with  the  Green's  function,  (4) 
determine  the  power  carried  away  by  the  deflected  wave,  and  (5)  derive  an 
expression  for  a  by  invoking  conservation  of  power  for  the  incident,  trans¬ 
mitted  (within  the  core,  past  the  DBD),  and  deflected  waves.  Once  a  is  known, 
the  deflector  efficiency  and  deflected  wave  far-field  patterns  can  be  found. 
The  analysis  will  be  carried  out  for  the  both  TE  (electric  fields  polarized 
parallel  to  the  x-y  plane)  modes  and  TM  (magnetic  fields  polarized  parallel  to 
the  x-y  plane).  The  angular  deflection  range  to  be  considered  will  be  that 
for  which  rr / 4  <  0  <  3w /4 . 
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F. 3  DEFLECTED  WAVE -FIELD  CALCULATIONS 

The  TE  and  TM  vector  Helmholtz  equations  describing  wave  propaga¬ 
tion  within  the  DBD  structure  are  for  TE 

-d  <  z  <  o 

otherwise  (F-l) 


('2+"?k  o=)- ! 


TE 


and  for  TM 
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(,2  +  \  ko)  5  ■ 


-d  <  z  <  o 


otherwise 


(F-2) 


where 


^TE  +  "koAn2(?>  E  -  7  |e  7  log  [eQn2  +  eQAn2(?)j  j  (F-3) 


=  -ko^n2 (r )H  -  7  log  [^n2  +  EQAn2(r)  x  7  x 


(F-4) 


In  Eqs.  ( F — 1 )  through  (F-4),  kQ  is  the  magnitude  of  the  free  space  wave- 

vector,  e  is  the  permittivity  of  free  space,  n^  refers  to  the  refractive 

°  2  ■+■ 
index  of  medium  "i"  (i  =  1,  2,  3,  4),  "7"  is  the  del  operator,  and  An  (r)  the 

periodic  waveguide  perturbation.  For  the  square-wave  surface  corrugation 

under  consideration, 


An2(r) 


(CO 

\"1  n2j)2  it  2-J  (24  -  i: 


cos  [(24-l)i<  •  r]. 
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where  r  is  the  position  vector  ( I r I  2  =  x2  +  y2  +  z2),  and  K  is  the  fundamental 
mechanical  phase  vector  of  the  perturbation  ( I K 1  =  2m/A;  A  =  perturbation 
period). 

The  electric  and  magnetic  fields  that  satisfy  Eqs.  (F-l)  and  (F-2) 
are  composed  of  core-  and  slab-propagating  modes  that,  in  accordance  with  the 
procedure  outlined  in  Section  F.2,  can  be  written 


*  +C  c 

FJ=E  +  E  =  E  (y,z)  exp 


.  is 
a  +  E 
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(F-6a) 
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and 


H  =  HC  +  HS  =  Hc  (y,z)  exp  [-i  x] 


a  +  HS 


(F-6b) 
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where  a  is  a  unit  vector  in  the  y-direction,  E  (or  H  )  is  the  assumed  field 
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distribution  of  the  core-propagating  mode,  E  (or  H  )  is  the  unknown  field 

distribution  of  the  Bragg  deflected  slab-propagating  mode,  and  the  core- 

~c  ~c 

propagating  phase-vectors  k^^  and  k^  have  been  taken  to  be  complex  in  order 
to  account  for  the  attenuation  that  results  from  the  Bragg  scattering  process. 
The  core-propagating  modes  Ec  and  Hc  (following  Marcatili's  analysis)  16  are 
assumed  to  be  well  enough  confined  so  that,  to  a  good  approximation,  their 
field  distributions  can  be  separated  in  the  y  and  z  coordinates  as  follows: 


E  (y,z)  = 


W(y)  £ ( z ) 


.  2 1€ 

Hc  (y , z )  =  W(y).*f(z) 

\ /  kTM 


(F-7a) 
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In  Eqs.  (F-7a)  and  (F-7b),  Ec  and  Hc  have  been  normalized  to  reflect  the 
assumption  that,  prior  to  entering  the  DBD  structure,  each  mode  carried  1  W  of 
power  in  the  positive  x-direction.  Also,  k^  and  k^  |=  Re  |  k^)  and  Re^k!^), 
respect ivelyj  have  been  taken  to  be  the  magnitudes  of  the  phase  vectors  of 
core-propagating  modes  in  the  absence  of  any  mode  conversion. 

W(y ) ,  the  transverse  field  variation  in  the  y-direction  for  either 
TE-  or  TM-polarized  core-propagating  modes,  is  approximated  as 


W(y)  + 


/I  cos  M 
w  \w  / 


(F-8) 


where  w  equals  both  the  equivalent  width1^  of  the  actual  field  distribution 
and  the  physical  extent  of  the  periodic  DBD  perturbation  in  the  y-direction. 
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(Errors  in  later  field  calculations  arising  from  this  approximation  are 
expected  to  be  small  for  modes  that  are  reasonably  well  confined  in  the  y- 
direction.)  g(z)  and  .#(2),  the  TE-  and  TM-mode  field  variations,  respec¬ 
tively,  in  the  z-direction,  are  found  by  applying  the  appropriate  boundary 
conditions  to  the  solutions  of^ 


(F-9a) 


^4^)2-(*™)2]: 


(z)  =  0;  i  =  1,2,3 


(F-9b) 


which  result  from  substituting  Eqs.  (F-7a),  (F-7b),  and  (F-8)  into  the  scalar 
Helmholtz  equation.  Because  of  the  form  chosen  for  W(y)  and  the  normalization 
imposed  by  (F-7a)  and  (F-7b),  e(z)  and  .'J{{ z)  are  themselves  normalized  as 
follows 


e  (z)dz  = 


•X2(z) 


dz  =  1;  i  =  1,2,3 


(F-10) 


With  Eqs.  (F-5),  (F-6a)  and  (F-6b)  used,  the  TE  and  TM  driving 

terms  of  Eqs.  (F— 1 )  and  (F— 2 )  are  now  calculated  to  be 
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where 


\  =  2(  -1-  -  log  /— )  [(2)1-1)  K  K  +  k^  K  1 

x  X  6  \n2  /  *-  x  y  TE  yJ 


A  -<HL 

y  * 


r\2  /  2  2  \ 

ko(nl  "  n2  ' 


(2,-d  +  2(2t-l)  Ky  log 


£) 


B  kc  lo6(2>) 

x  7T  y  TM  yn2  f 


(F-12a) 


(F-l 2b) 


(F-l 2c ) 


and 


B  =  - 

y 


„  r,  2  ,  2  2\ 

(-1)*  |ko(nl  -  °2 ) 


(2*-l) 


+  2 


k™  Kx  (si) 


(F-12d) 


In  deriving  Eqs.  (F-lla)  through  (F-12d),  it  has  been  assumed 

that  1  kTE ,TM 1  "  kTF,TM  and  that’ 
and  deflection  angles  of  interest 


that  |k^„  _ .  |  »  _w  and  that,  for  those  wavelengths,  waveguide  dimensions 
TE,TM  TF,TM 


|_  EC(y,z)  |_  EC(y,z)  «  Kx>y 


(F-l 3) 


The  condition  for  Eq .  ( F—  1 3 )  is  most  seriously  violated  in  the  case  of  long 

wavelengths,  very  small  transverse  waveguide  dimensions  (comparable  to  or 
smaller  than  a  wavelength),  and  either  very  shallow  or  very  large  deflection 
angles. 
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Given  analytical  expressions  for  QTE  and  Q^,  we  now  seek  suitable 
TE-  and  TM-mode  dyadic  Green's  functions  that  are  solutions  of  the  following 
inhomogeneous  wave  equations: 


(v2  +  n2k2j  Gte  (r i r " )  =  I<  (r-r' ) ;  i  =  1,2, 3, 4  (F-14a) 

and 

(v2  +  n^l-TM  (r  I  ?" )  =  I  6  (r-r)  ;  i  =  1,2, 3, 4  (F-14b) 


where  I  is  the  idem  factor  and  6(r  -  r')  is  the  three-dimensional  delta 
function. 

The  dyadic  Green's  function,  (1,  is  a  solution  of 

V2G  +  n2k2  C  =  I  {  (r  -  r' ) ;  i  =  1,  2,  3  (F-15) 

Given  the  simple  boundary  conditions  that  £  must  satisfy  (for  either  TE-  or 
TM-polarized  waves),  £  can  be  written 
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where  e  is 
m 

the  scalar 
cylindrical 


any  set  of  curvilinear  coordinate  system  unit  vectors,  and  G  is 
Green's  function.  Expressing  Eq .  (F-15)  in  terns  of  circular 
coordinates  and  making  the  following  coordinate  transformation 
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the  wave  equation  satisfied  by  G  becomes 
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Since  G  is  now  recognized  as  being  a  function  only  of  R  and  z, 
solutions  of  Eq.  C F— 17)  within  the  source-free  region  nay  be  found  by  writing 
G(R,z)  =  f(R)h(z)  and  using  separation  of  variables,  which  yields 


I  ±_ 

R  dR 


(F-18  ) 


and 


~  +  (n2  k2  -  k2)  h  =  0  (F—  1 9  ) 
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where  k“  is  the  separation  constant  and  the  solutions  of  Eq,  (F— 19)  are  recog¬ 
nized  as  being  either  £  (z)  (for  the  TE-polarized  waves)  or  ,,(z)(for 

o,v  r  o,v 

the  TM-polarized  waves). 

The  admixture  of  h  and  f  required  to  construct  G  can  be  found  by 
integrating  Eq.  ( F  —  1 7 )  over  the  cylindrical  volume  depicted  in  Fig.  F-3.  By 
ignoring  contributions  from  the  cylinder  ends  (valid  for  long,  slender 
cylinders),  this  integration  yields 


z=z  - 


(F-20) 


where  use  has  been  made  of  the  divergence  theorem,  and  e  is  the  radius  of  the 
cylinder.  From  Eq.  (F-20),  it  is  implied  that 

2*e  ||  I  =  <5 ( z  -  z')  (F-21) 
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or,  expressing  6(z  -  z')  in  terms  of  the  complete,  orthonormal  sets, 
eQ  v(z)  and  ,>4^  y(z),  that 
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for  TE-polarized  modes  and 
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for  TM-polarized  modes. 

Finally,  integrating  Eqs.  (F-22)  and  (F-23)  over  R  and  imposing  the 
necessary  discontinuities  on  GTE  and  as  £  +  R  +  0  yields  GTE  and  G^  as 
given  in  Eqs.  (F~24a)  and  (F-24b). 
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where  R  =  I  r— r"  | ,  H^^(x)  is  the  zeroth  order  Hankel  function  of  the  second 

kind;  e  (z)  and*>r(z)  are  the  lowest-order  solutions  of  Eqs.  (F-9a)  and  (F-9b) 
o  o 

(we  shall  be  concerned  with  guides  that  confine  only  the  lowest  order  core- 
propagating  modes);  e^(z)  andJS^(z)  are  the  continuum  mode  solutions  of  Eqs. 
(F-9a)  and  (F-9b)*®;  k^,  and  k^  are  the  magnitudes  of  the  corresponding  con¬ 
tinuum  mode  phase  vectors;  and  j  =  1,  2,  3.  G__  and  G__,  (which  lend  them- 
selves  naturally  to  the  interaction  geometry  of  Fig.  F-l  are,  as  noted, 
approximate  in  the  sense  that  Bragg-scattered  wave  components  reflected  a\  the 
core-slab  interface  have  been  ignored.  This  approximation,  which  requires 

that  n.  "  n„  and,  hence,  that  e  (z),  e  (z),  CK  (z),  and  -X  (z)  retain 
4  2  o  v  o  v 

essentially  the  same  shape  throughout  the  slab-coupled  structure,  remains 
valid  for  all  but  very  shallow  or  very  wide  deflection  angles. 

Knowledge  of  the  dyadic  Green's  functions  Gj-g  and  GjM  and  of  the 
driving  terms  ($TE  and  finally  permits  calculation  of  the  various  field 

quantities  by  evaluation  of 


-*■  -*C 

E  =  E 


+  ES 


-//A, 


(rlr* 


)  Oj-g  (EC  +  ES)  dx'  dy'  dz'  (F-23a) 


and 


+ 
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&C  +  HS 


-///$* 


(rl?- 


(»C  +»S) 


dx'  dy'  dz' 


(F-25b) 


Equations  (F-25a)  and  (F-25b)  are  simplified  considerably  because  the  DBD  per¬ 
turbation  induces  coupling  only  between  confined  modes  (core-  and  slab- 
propagating  modes)  and  rescattering  from  the  slab-propagating  modes  back  to 
the  core-propagating  mode  can  be  ignored.  Accordingly,  Eqs.  (F-25a)  and 
(F-25b)  become 


<^HTe  <EC)  dx'  dy'  dz' 


(F-26a ) 
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and 


HS 


■Iff* 


P 

-TM 


(r|r') 


Qtm  (Hc>  dx' 


dy'  dz' 


(F-26b) 


Where  the  only  terms  retained  from  the  right-hand  sides  of  Eqs.  (F-25a)  and 
(F-25b)  have  been  those  which  contribute  synchronously  (or  nearly  synchro¬ 
nously)  to  the  slab-propagating  fields,  Es  and  Hs,  and  where  Gp  and  GP  rep- 

— Ik  — TM 

resent  partial  dyadic  Green's  functions  obtained  by  dropping  all  but  the  first 
(confined  mode)  terms  on  the  right-hand  sides  of  Eqs.  (F-24a)  and  (F-24b). 

Evaluation  of  the  right-hand  sides  of  Eqs.  (F-26a)  and  (F-26b) 
yields  for  the  slab-propagating  field  distributions 


(Oil 


s  ,  i  /  —  e*P  + 
2k  V 


H  -  ° !)  [fd  «20(«-w-  e0(z)  exp  ( jtx)  f(e> 


[u^  (A^  cos  0  +  Ay  sin  9)  +  Ug  (-A^  sin  ®  +  \  cos  0)1>  j  *  1,2,3 


(F-27a) 


and 


♦.  i  /!%  /■  l \/-r°  liz"> . : 

2k  J  *  exp  (  4  2/  (/d  n2  (Z-)  Z  , 


(.)  f(0) 


*  (u^(Bv  cos  0  +  Bu  sin  0)  +  ufl(-Bw  sin  0  +  Bw  cos  8)1;  j  -  1,2,3 

(F-27b) 
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where  E  and  H  have  been  expressed  in  terms  of  the  circular  cylindrical 

A  A 

coordinate  system  indicated  in  Fig.  F-l  (u  ^and  u  q  are  the  radial  and  polar 
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unit  vectors  for  this  coordinate  system);  the  large-argument,  asymptotic  form 

(2) 

for  H  (x)  has  been  used 
o 


-i(x-n) 

A 


(F-28) 


the  complex,  core-propagating  phase  vector  k  has  been  taken  to 
be  k  =  k  -  ia,  and 


f  C  0 )  =  ff  W(y")  exp  [-iu(0)y"]  exp  [-ax"  -  iv(0)x"]  dx"  dy"  (F-29) 


For  clarity,  all  TE  and  TM  subscripts  have  been  omitted  in  Eqs.  (F-27a) 
through  (F-29)  (the  identities  of  all  TE-  and  TM-related  quantities  are 
evident  from  context)  and  the  angle-dependent  factors,  u(0)  and  v(0),  which 
appear  in  Eq.  (F-29),  are  given  by 


and 


u(0)  =  (2A  -  1)  Ky  -  k  sin  0 


(F-30) 


v(0)  =  k  (1  -  cos  0)  +  (2*  -  1)  Kx  (F-31 ) 

where  and  Ky  are  the  Cartesian  components  of  the  DBD  corrugation  mechanical 
phase  vector. 

In  order  to  calculate  the  core-propagating  attenuation 

coefficients  and  of  Eqs.  (F-27a)  and  (F-27b),  respectively,  we  now 

invoke  the  principle  of  power  conservation  among  the  incident  and  transmitted, 
core-propagating  modes  and  the  spectrum  of  Bragg-deflected  slab-propagating 
modes.  Specifically,  we  first  integrate  Poynting's  vector  for  the  deflected 
field  distribution  over  a  cylindrical  surface  whose  axis  of  symmetry  corre¬ 
sponds  to  that  of  the  cylindrical  coordinate  system  indicated  in  Fig.  F-l  and 
whose  radius  is  large  compared  with  L,  the  length  of  the  DBD.  We  then  set  the 
result  of  this  integration  equal  to  one  watt  (the  assumed  power  content  of  the 
incident  core-propagating  mode)  less  the  power  transmitted  down  the  core  past 
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the  DBD  structure.  The  resultant  equations  yield  explicit  expressions  for 
both  aTE  and  o^.  Accordingly,  the  power  conservation  equation, 


Y  J'J' Re  [E  *  H  ]  rd6  z  =  1  -  exp  (- 


2aL) 


(F-32) 


becomes,  for  the  TE-  and  TM-polarized,  deflected  fields: 


r2 

_TE_  fr 

8ltkTE  J 


■A  sin  9  +  A  cos  0]2  |f(0)|2  d9 
x  y 


1  -  exp  (-2nTE  L) 


(F-33a) 
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(F-35) 


and 
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cos  [u  (9)  w/2]  2 
l  -  lu  (0)  w/it]2 


(F-36) 
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As  before  the  specific  definition  of  a  (i.e.,  either  a  or  a  )  in  Eqs. 

TE  In 

(F-32)  and  Eqs.  (F-24)  through  (F-36),  is  understood  from  context. 

The  integrals  on  the  left-hand  sides  of  Eqs.  (F-33a)  and  (F-33b) 

are  easily  evaluated  once  the  relative  variations  with  respect  to  6  of  the 

2 

Integrand  terms  g  -(9),  g  >(6),  (-A  sin  9  +  A  cos  9)  ,  and  (-B  sin  9  +  B 
2  x  y  x  y  x  y 

cos  9)  have  been  examined.  The  full  widths  at  half  maxima  of  gx-(9)  and 

g  -(9)  (which  are  plotted  in  Fig.  F-4  for  the  case  of  exact  satisfaction  of 

the  Bragg  condition)  are  given  by 


A0  „  ■  - 

x  k  s<n  (6g) 


2* 

k  sin  (9„ 


(F-37) 


40^  -  <2 


9^  near  x/2 


w  k  cos  (£ 


otherwise 


(F-38) 


where  9  ,the  Bragg  scattering  angle,  is  recognized  as  being  the  deflection 
angle  for  which  u(9g)  =  v(9^)  =  0  and,  hence,  the  angle  for  which  g^,  and 
gy,  are  simultaneously  maxima  (the  latter  conclusion  requires  that  k  »  a, 
a  condition  that  will  always  be  met  for  those  geometries  under  consideration). 
From  Eqs.  (F-37)  and  (F-38),  we  now  note  that,  for  9  within  the  angular  range 

D 

of  interest  (i.e.,  3x/4  >  9  >  x/4)  and  for  weak  deflectors  (i.e.,  k  »  «), 

both  (-A  sin  9  +  A  cos)  and  (-B  sin  0  +  B  cos  9)4  vary  slowly  with 

x  y  x  y 

respect  to  both  g  „  (9)  and  g  ,  (0)  and  can,  as  a  consequence,  be  drawn  out- 
x  y 


side  the  integrals  of  Eqs.  (F-33a)  and  (F-33b). 


Similarly, 


for  a  »  w.  40  ^  «  40  *  and  e  -  (9„)  can  be  withdrawn  from  the  integrals, 
*  x  y  °y  B 

leaving  g  ^  (9)  as  the  only  term  that  must  be  integrated  in  Eqs.  (F-33a)  and 
(F-33b)  and  the  term  which,  because  of  its  origin  Eq.  (F-32),  determines  the 
far-field  pattern  of  the  deflected  wave. 


i 


Before  performing  the  integrations  in  Eqs.  (F-33a)  and  (F-33b)  in 
order  to  obtain  a_„  and  a_,,  a  comment  should  be  made  regarding  the  angular 
behavior  of  (0)  which,  as  noted  earlier,  determines  the  far-field  radia¬ 
tion  pattern  of  the  DBD.  Accordingly,  note  that  g  ..  (0)  [as  well  as  g  ,  (0)] 

x  y 

is  double  valued  over  the  angular  range  for  which  2tt  >  0  >  0  (Fig.  A-3),  re¬ 
flecting  the  possibility  that  more  than  one  deflected  beam  may  emerge  from  the 
DBD.  By  referring  to  Eq.  (F-36),  this  possibility  is  eliminated  by  requiring 
that 


u  (-0B)  5 
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3ir 

2 


or  that 


w  > 


_ 3tt 

2k  sin  (0_) 

D 


(F-39) 


The  physical  basis  for  Eq.  (F-39)  is  the  requirement  that  the  DBD  be  made  wide 

enough  to  permit  the  core-propagating  mode  to  "sense"  the  transverse  component 

of  the  perturbation  phase  vector  Ky,  and  to  prevent,  thereby,  energy  from 

being  deflected  away  from  the  structure  at  an  angle  -0D.  The  condition  in  Eq. 

-1  “ 

(F-39)  together  with  the  requirement  that  w  <<  a  ,  are  seen  to  constrain  w  to 
lie  within  a  well  defined,  but  large  and  useful  range.  For 
example,  0O  =  x/4,  a  *  =  1  mm,  and  k  =  20  pm  *  (for  ~  0.8  pm  light  propagating 

D 

in  a  GaAs  waveguide)  requires  that  1000  pm  »  w  >  0. 3  pm. 

With  the  use  of  the  techniques  and  approximations  discussed  earlier 
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and  of  the  well-confined  mode  approximation,  *  which  implies  that  k  B  ^k 
and  that 


r 


TE 


TM 


Equations  (F-33a)  and  (F-33b)  finally  yield  for  the  TE  and  TM  core-propagating 
attenuation  coefficients 
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which,  normalized  by  the  factor 


(2A  -  1) 

v  =  - - — 
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(^)6 


(F-A2) 


(where  t  and  d  are  the  waveguide  thickness  and  perturbation  groove  depth,  re¬ 
spectively)  have  been  plotted  in  Fig.  A-5  for  a  GaAlAs  DBD  with  n^  =  1, 
x\2  -  3.6,  and  =  0.85  pm.  Equations  (F-AO)  and  (F-Al)  and,  hence,  the 
methods  used  in  their  derivations  have  been  checked  by  allowing  the  DBD  width 
to  approach  infinity  and  considering  the  limiting  case  of  180-deg 
deflection.  Under  these  circumstances,  A8^  «  A0x>  >  gx^(*)  Can  ^  w^t'1<^rawn 


from 


the 


integrals 


of 


Eqs. 


(F-33a)  and  (F-33b);  and  the 


resultant  (0_  =  it)  and  <*_,  (0D  =  it)  reduce  to  the  familiar  contradirec- 
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tional  coupling  coefficients  of  the  DBR. 

As  discussed  in  Section  F.2,  knowledge  of  and  permits 
calculation  of  the  DBD  deflection  efficiencies,  far-field  deflected  beam 
patterns  and  far-field  beam  widths.  Accordingly,  the  DBD  deflection  effi¬ 
ciencies,  for  the  TE-  or  TM-polarized  waves,  are  given  in  a  straightforward 
manner  by 


nTE,TM  1  “  eXp  ^“2  LaTE,TfP 


(F-A3) 
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Moreover,  for  aL  ^  2  and,  as  discussed  earlier,  A0 
mode  far-field  patterns  and  far-  field 
by  (0)  and  A0^,,  respectively,  or 


.  «  A0  ,,  the  TE-  and  TM- 
x  y  ’ 

beam  widths  are  given 


and 
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(0tTE, 
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_  Afl  =  ™ 
TM  ~  x'  2k  sin  (0D) 
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where  Iq-£  ^  has  a  Lorentzian  shape  with  a  full-width  at  half 
maximum  A®TE  ^  approximately  equal  to  a  wavelength  within  the  guiding  medium 
divided  by  the  e-folding  length  of  the  DBD  (~  1  / The  latter  observa¬ 
tion  is  not  surprising  since  the  Fourier  transform  (Fraunhofer  diffraction 
pattern)  of  a  decaying  exponential  (the  near-field  pattern  of  the  deflected 
wave  must  bear  the  same  exponential  dependence  in  the  x-direction  as  the  inci¬ 
dent  wave  by  conservation  of  power)  is  a  Lorentzian. 

We  conclude  by  considering  the  effects  on  a  of  small  departures 
from  exact  satisfaction  of  the  Bragg  condition.  Specifically,  given  that 

K  =  2k  sin  =  2k  cos  (<j>)  (F-46) 

for  exact  satisfaction  of  the  Bragg  condition,  we  consider  the  reductions 

in  aTE  and  that  occur  as  a  result  of  small  errors  in:  (1)  the  magnitude 

of  the  DBD  grating  phase  vector  Ak,  (2)  the  angular  alignment  of  the  DBD 

grating  phase  vector  A<J> ,  and  (3)  the  magnitude  of  the  core-propagating  phase 

vector  Ak.  With  it  assumed,  as  before,  that  A0  „  «  A0  „,  it  is  easily  shown 

x  y 

that  the  reduced  attenuation  coefficients  olL  and  olL  are  given  by 


aTE,  TM  =  aTE,  TM  8y'  (U  W) 


(F-47 ) 
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where 
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AK  -  k  A<j> 
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cos  (efi)  -  1 
sin  (8g) 
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(F-48) 


The  reduction  factor,  g  „  (u  w),  which  equals  one  when  condition  (F-46)  is 
satisfied,  is  plotted  in  Fig.  F-6. 

F.4  DEVICE  APPLICATIONS 

The  operational  characteristics  of  the  basic  DBD  structure  suggest 
that  a  number  of  active  and  passive  integrated  optical  device  applications  are 
possible.  Passive  applications  include  intrawaveguide  beam  expansion,  divi¬ 
sion,  deflection,  and,  for  90-deg  deflection,  polarization  analysis.  Pos¬ 
sible  active  applications  include  modulation  and  wide-angle,  intrawaveguide 
beam  steering. 

Used  as  a  lens-less  intrawaveguide  beam  expander-collimator,  the 

DBD  appears  particularly  well  suited  for  reducing  the  beam  divergence  of  fully 

integrated  DFB  or  DBR  lasing  sources  (wherein  the  expansion-collimation  would 

take  place  in  a  passive  waveguiding  region  lying  adjacent  to  the  active, 
2 1 

lasing  region).  Integrated  optical  device  applications  of  such  a  laser-DBD 

22 

configuration  might  include  real-time  spectrum  analyzers,  which  require 

high-resolution  acousto-optic  beam  scanning  and  high-radiance  (or  low  beam 

divergence)  semiconductor  laser  sources.  As  an  example,  taking  the  case  of  a 

GaAlAs  DBD  for  which  n,  =  I,  n5  =  3.6,  A  =  0.85  um,  8  =  tt/2,  l  =  1  (A  = 

0.167  pm),  w  =  10  pm,  t  =  2.0  pm,  d  =  0.136  pm,  L  =  2  mm,  and  assuming  the 

well-confined  mode  approximation  is  valid,  Fig.  F-4  and  Eqs.  (F-43)  and  (F-45) 

— 

yield  a  deflected  wave  beam  width  of  A9  ^  =  1.5  x  10  rad  or  an 

x 

intrawaveguide  beam  expansion  factor  of  M  =  157,  where  M  is  defined  by 

(X  /n  _  w) 

M  s  A6  :  - 
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As  an  efficient  intrawaveguide  polarizer-analyzer,  the  DBD's  potential  utility 

is  immediately  suggested  by  the  plots  of  (0g)  and  (0g)  shown  in  Fig. 

F-4.  Operation  of  a  DBD  near  0„  *  ir/2  rad  ensures  both  good  polarization  dis- 

0 

crimination  (a^  is  roughly  a  factor  of  5  smaller  than  near  0^  -  n/2  and, 

indeed,  equals  zero  for  9  -  n/2)  and  considerable  insensitivity  to  small 

B 

errors  in  0  (a™  goes  through  a  broad  null  near  0n  *  n/2).  The  DBD  may  be 

used  as  a  simple  intrawaveguide  beam  deflector  or  beam  divider.  For  example, 
with  L  =  0. 2  mm  and  all  other  device  parameters  as  before  (for  the  beam  ex¬ 
pander)  yields  a  TM  deflection  or  beam  division  efficiency  of  n^  =  0.55. 

With  regard  to  use  of  the  DBD  as  a  modulator,  Fig.  F-6  and  Eqs. 
(F-47)  and  (F-48)  indicate  that  (with  modulation  assumed  to  be  accomplished 
through  an  external  manipulation  of  either  a^E  or  a^j)  only  relatively  large 
changes  in  AK  and  A$  are  effective.  From  this  observation,  it  is  implied  that 

O 

a  DBD  modulator  based  upon  electronically  alterable  magnetic  stripe  domains, 
as  an  example,  is  feasible  whereas  one  thct  relies  upon  a  relatively  small 
electro-optic  variation  in  Ak  is  not.  Similar  reasoning,  of  course,  applies 
to  the  case  of  wide-angle,  intrawaveguide  beam  scanners  that  require  large 
simultaneous  changes  in  both  AK  and  A$. 

F. 5  CONCLUSIONS 

A  Green's  function  technique  has  been  used  to  find  the  deflection 
efficiencies  and  far-field  deflected-beam  patterns  of  an  integrated  optical 
wide-angle  deflection  device, a  distributed  Bragg  deflector  (DBD).  Intrawave¬ 
guide  interactions  involving  deflection  angles  between  x/4  and  3n/4  rad  and 
both  TE-  and  TM-polarized  waves  have  been  considered. 

Because  of  the  gradual  or  distributed  nature  of  the  Bragg  deflec¬ 
tion  process  for  those  geometries  of  interest,  the  DBD  may  be  employed  as  an 
efficient  compact  intrawaveguide  beam  expander.  Moreover,  because  of  its 
polarization-dependent  deflection  efficiencies,  the  90-deg  DBD  is  promising  as 
a  simple  integrated  optical  polarizer-analyzer.  Other  potential  passive  inte¬ 
grated  optical  embodiments  of  the  basic  DBD  structure  include  intrawaveguide 
beam  deflectors  and  beam  splitters.  With  regard  to  active  device  applica¬ 
tions,  the  relative  insensitivity  of  the  DBD's  deflection  efficiency  to  small 
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changes  in  phase  vector  magnitude  appears  to  preclude  use  of  the  DBD  as  an 
electro-optic  Bragg  modulator.  Integrated-optical  modulators  and  intrawave¬ 
guide  beam  scanners  based  upon  an  external  manipulation  of  either  the  magni¬ 
tude  or  the  orientation  of  the  DBD  k-vector  do,  however,  appear  feasible. 
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LABORATORY  OPERATIONS 


Th*  Laboratory  Operation*  of  Th*  Aaroapaca  Corporation  la  conducting 
aaparlaantal  and  thaoratlcal  InvaatlRatlona  nacaaaary  for  the  evaluation  and 
application  of  aclantlflc  advancaa  to  new  allltary  concept*  and  ayataae.  Var- 
aatlltty  and  flaalblllty  have  bean  developed  to  a  high  dagraa  by  the  laborato¬ 
ry  paraonnal  In  dealing  ulth  the  eany  problea*  ancountarad  In  th*  Nation'* 
rapidly  developing  apac*  ayataaa.  Exparti**  in  th*  latent  aclantlflc  davelop- 
nanta  la  vital  to  th*  acconpl lehaant  of  taaka  related  to  th***  problan*.  Th* 
laboratorla*  that  contribute  to  tht*  raaearch  era: 

Aaropbyalca  Laboratory:  Aarodynaalce;  fluid  dynaalcai  plaanadynaalce; 
chanlcal  klnactca;  engineering  aachanlca;  flight  dynaaica;  heat  tranafar; 
high-powar  gaa  laaara,  contlnuoua  and  pulaed,  IR,  vlalbla,  UV;  laaer  phyalca, 
laaar  raaonator  optica;  laaer  affect*  and  eountaraeaaurae. 

Chanlatry  and  Phyalca  Laboratory;  Ataoapharlc  react lona  and  optical  baefc- 
g round*;  radiative  tranafar  and  ataoapharlc  tranaalaalon;  tharnal  and  atata- 
apaclflc  reaction  rate*  In  rocket  plunaa;  chanlcal  thernodynaalca  and  propul - 
alon  chanlatry;  laaar  laotope  aaparation;  chanlatry  and  phyalca  of  partlclaa; 
apaca  aavlronnaaral  and  contanlnattoa  affacta  on  apacacraft  aatarlala;  lubrica¬ 
tion;  aurfaca  chanlatry  of  lnaulatora  and  conductora;  cathode  natarlala;  aan- 
aor  notarial*  and  aanaor  optica;  applied  laaar  apactroacopy;  atonic  frequency 
atandardai  pollution  and  toalc  aatarlala  aonltorlng. 

Electronic*  Raaearch  Laboratory;  Rlactrnaagnatlc  theory  and  propagation 
phanonana  ;  alcrouav*  and  aanlconductor  device*  and  Integrated  circuit*;  quan- 
Cun  electronic*,  laaara,  and  elect ro-opt lea;  coanunieatloo  aclancaa,  applied 
alactronlca,  aupar conduce  lag  and  electronic  device  phyalca i  al  1 llnater-vava 
and  far- Infrared  technology, 

Hatatlala  Science*  Laboratory:  Developaant  of  nan  aatarlala;  coapoalt* 
aatarlala;  graphite  and  caranlca;  polynarlc  natarlala;  waapoaa  affect*  and 
hardened  natarlala;  natarlala  for  electronic  device*;  dlaanalonally  atabl* 
aatarlala;  chanlcal  and  atructural  analyaaa;  atraaa  corroaloa;  fatigue  of 
aatala. 

Space  Science*  laboratory;  Ataoapharlc  and  toaoapherlc  phyalca,  radia¬ 
tion  tran  the  afanaphara,  danaTty  and  conpoaltlon  of  the  atnoaphara,  aurora* 
and  alrglou;  aagnatoapharlc  phyalca,  coanlc  ray*,  fan* radon  and  propagation 
of  plaan*  uavaa  la  th*  angnatonphara;  aolar  phyalca,  i-ray  aatronony;  th*  affacta 
of  nuclear  a^laalona,  angnatlc  atone,  and  aolar  activity  on  th*  earth'* 
atnoaphara,  lenaaphara,  and  angnatoaphora;  the  affacta  of  optical,  electronag- 
aatlc,  and  particulate  radiation*  la  apnea  on  apac*  ayatana. 
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